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Abstract 

The application of the Correlated basis function theory and of the Fermi hypernetted 
chain technique, to the description of the ground state of medium-heavy nuclei is reviewed. 
We discuss how the formalism, originally developed for symmetric nuclear matter, should 
be changed in order to describe finite nuclear systems, with different number of protons and 
neutrons. This approach allows us to describe doubly closed shell nuclei by using microscopic 
nuclcon-nucleon interactions. We presents results of numerical calculations done with two- 
nucleon interactions of Argonne type, implemented with three-body forces of Urbana type. 
Our results regard ground-state energies, matter, charge and momentum distributions, nat- 
ural orbits, occupation numbers, quasi-hole wave functions and spectroscopic factors of ^^C 
, , 40Ca , ^^Ca and ^ospb nuclei. 



PACS: 21.60.-n; 21.10.Dr; 21.10.Ft; 21.10.Gv; 21.10.Jx; 21.90-Ff 



1 Introduction 



Aim of the many-body theories is the description of composite systems in terms of their ele- 
mentary components. In our present understanding of nature, only leptons, quarks, and gauge 
bosons are considered to be elementary. In principle, any composite system should be described 
in terms of these entities. In reality, pragmatical difficulties hinder the accomplishment of such 
an ambitious program. For example, the difficulties in dealing with the non perturbative features 
of Quantum Chromodinamics (QCD) complicates the description even of the individual hadrons. 
The description of even more complex systems in terms of quarks and gluons is evidently not 
practicable, at least at present. 

A more pragmatical, and fruitful, approach for the study of composite systems abandons 
quarks and leptons and uses effective degrees of freedom and interactions. For nuclear systems, 
the most convenient choice is to consider the nucleon as the basic degree of freedom. The nucleon- 
nucleon (NN) interaction, even if modeled in various manners, is always fixed to reproduce the 
NN elastic scattering data and the deuteron properties [H El O lU [5] . 

After choosing the basic degrees of freedom and their interaction it is necessary to define the 
theoretical framework to use. In the case of nuclear systems, if one is interested in the ground 
state properties and in low energy phenomenology (we mean here energies well below a GeV) 
the role of antiparticles can be neglected. Therefore a good description of these systems can be 
provided by the Schrodinger equation. 

In the last decade the Schrodinger equation has been solved without approximations by using 
Green function Montecarlo (GFMC) techniques, for several light nuclei up to A=12 [El 13 EJ [9] . 
In these calculations the two-body interaction has been implemented with a three-body force 
chosen to reproduce the triton binding energy. The great success in describing binding energies 
and low-lying spectrum of these nuclei, is the pragmatical demonstration of the validity of the 
non relativistic many-body approach. 

Various techniques to solve the many-body Schrodinger equation without approximations 
have been developed. Some of them are tailored to describe only few-body systems [10]. Others, 
formulated to handle any number of nucleons, like the GFMC, are limited for computational 
reasons to deal with light nuclei. Recently, a new Montecarlo approach, called Auxiliary Field 
Diffusion Montecarlo (AFDMC) [H], has been developed and it shows potentialities to be applied 
to the description of medium-heavy nuclear systems. 

The many-body theories try to overcome the difficulties in solving the many-body Schrodinger 
equation by using approximations which simplify the problem, but still provide a proper descrip- 
tion of the relevant physics effects. Because of the strong repulsion in the scalar channel of the 
NN interaction at short internucleonic distances the use of the most traditional, and simplest, 
approximations, such as the mean-field approximation, fails badly. For example, the use of mi- 
croscopic interactions in Hartee-Fock (HF) calculations leads to unbound nuclei [12]. Dealing 
with the strong repulsion at short distances is the major issue of the nuclear many-body theories. 

Loosely speaking the various nuclear many-body theories can be classified in two categories 
depending upon how they treat the short-range repulsion problem. There are theories acting 
on the interaction, and others working on the trial many-body wave function. The Brueckner 
theory, and all the theories constructing effective interactions from microscopic ones, belong 
to the first category. One of the theories belonging to the second category is the Correlated 
Basis Function (CBF) theory whose most recent extensions and applications to medium-heavy 
nuclei will be presented in this report. Strictly speaking, our approach is the lowest order 
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approximation of the CBF theory |13| . 

The starting point of the CBF approach is the solution of the many-body Schrodinger equa- 
tion by means of the variational principle: 



= 0. 



(1) 



< > 



The search for the minimum is done by using trial wave functions of the form: 



=F(1,..,A)$(1,.... 



A) . 



(2) 



In the above equation <I> describes the system as a set of A particles moving independently from 
each other. We call Independent Particles Model (IPM) this picture, which, in our calculations, 
is modified by the correlation function F. In its easiest form, we use for F the expression |14] : 



where f{rij) is a scalar function of the distance between two particles of the system. 

The peculiarity of our approach consists in the technique used to calculate the expectation 
values of Eq. ([1]) . This technique is inspired by the cluster expansion method used in statistical 
mechanics to describe liquids [15]. The particles, correlated by the function /, form clusters. A 
topological study of the various clusters shows that it is possible to construct a set of integral 
equations which allows one to sum in a closed form the contributions of all the clusters with 
some specific topological properties. This set of integral equations, called HyperNetted Chain 
(HNC) equations, can be used to describe both classical and bosonic systems jTB] . 

In the mid 1970's, the cluster expansion techniques were extended to include also the Pauli 
exclusion principle, and the Fermi HyperNetted Chain (FHNC) equations were formulated [171 
[T8| [19] . The complexity of the NN interactions requires the use of correlation functions that are 
more complex than those of Eq. ([3|). These new correlations contain operator dependent terms 
which commute neither with the nuclear hamiltonian, nor among them. Also for this reason it 
became necessary to extend the FHNC equations to deal with this new type of correlations [20] . 
The computational difficulties require the use of an approximation called Single Operator Chain 
(SOC). The resulting set of equations is called FHNC/SOC [21], and it has been successfully 
applied to describe infinite systems [22 | [23 | [24 l [25]. 

In this review, we are concerned about the application of the FHNC/SOC computational 
scheme to medium-heavy nuclei. The extension of the FHNC theory to finite Fermi systems 
was introduced by Fantoni and Rosati in the late 1970s [26]. In their works they have shown 
that a cluster expansion with an infinite numbers of terms can be formulated even for finite 
systems. Consequently, the basic set of FHNC equations, can be used also for finite systems. 
However, one has to consider that the loss of translational invariance in these systems produces 
the so-called vertex corrections. We shall refer to the new set of equations as Renormalized 
Fermi HyperNetted Chain (RFHNC) equations. 

The results of the first numerical application of the RFHNC equations to finite nuclear 
systems were presented in Ref. [27] . In that article, model nuclei were described. Protons and 
neutrons wave functions were produced by a unique mean field potential, and in a Is coupling 
scheme. The NN interactions considered had only central terms, and the correlations were scalar 
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F(i,....,^)= n f^n,) 



(3) 



j>i=l 
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functions. This simplified situation was used to test the theoretical, and numerical, feasibility 
of the approach. Results for binding energies of ^^O and ^"Ca model nuclei were presented in 
[27] while the momentum distributions where shown in a following article |28j . 

A more realistic description of doubly closed shell nuclei was given in [29] , where proton and 
neutrons were separately treated, and the single particle wave functions were expressed in a jj 
coupling scheme. The RFHNC equations required a non trivial reformulation. Binding energies, 
matter densities and momentum distributions, have been calculated for various doubly magic 
nuclei up to ^^^Pb . However, also in this case, simple central interactions and scalar correlations 
were used. 

In a following step, the RFHNC equations were extended to treat the correlation terms 
commuting neither with the hamiltonian, nor among themselves. This involved the extension of 
the SOC approximation. Because of the technical difficulties the RFHNC/SOC equations have 
been first formulated to deal with spin and isospin saturated nuclei, and with single particle 
wave functions in a Is coupling scheme. Again only ^^O and ^"Ca nuclei could be treated. The 
results of these calculations have been presented in Refs. [30 } [311 [32] . 

A formulation of the RFHNC/SOC equations general enough to handle separately protons 
and neutrons in the more realistic jj coupling scheme was finally done. Binding energies and 
density distributions have been shown in Ref. [33] for the ^^C , ^^O , ^^Ca , '^^Ca and also ^osp^j 
nuclei. Here fully realistic microscopic interactions, with tensor and spin-orbit terms were used. 
The hamiltonian included also three-body interactions. 

In the literature, there are various reviews regarding the FHNC/SOC formalism applied to 
infinite nuclear systems [21^ [HH [351 [36] 137] . but there is a void regarding finite nuclei. For the 
sake of brevity in writing journal articles, the formalism presented in the papers quoted above is 
incomplete. The aim of the present article is to provide a complete, coherent, and self-contained 
presentation of the FHNC/SOC formalism for finite nuclear systems, and to review the most 
recent results. 

We recall in Sect. [2] the HNC, FHNC and FHNC/SOC equations for infinite systems. They 
are important, not only because we want to give a self-contained presentation, but especially 
because the RFHNC/SOC formulation for the finite systems is constructed by modifying that 
of infinite systems. The RFHNC/SOC set of equations will be presented in Sect. [3l and it will 
be applied in Sect. [H to evaluate the energy of the system. A selected, but significant, set of 
recent numerical results will be presented, and discussed, in Sect. [H In Sect. [6] we provide a 
short overview of the possible extensions of the formalism. Conclusions are presented in Sect. 

13 

To improve the readability of the paper we present many technical details of the derivation 
of the various expressions in the Appendices. Furthermore, because of the large use of acronyms 
and symbols, we list them in the Appendices iGl and [Hi respectively. 
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2 Infinite systems 



In this section we present the HNC and FHNC equations for infinite systems. This presentation 
does not aim to substitute, or update, the excellent review articles describing in detail the 
derivations of the various expressions, see for example Refs. [211 IMl ESJ [36l Elj. Our purpose 
is to recall the main ideas and to emphasize those details which should be reconsidered in the 
description of finite systems. 



2.1 Bosons 

We start to present the CBF approach by describing a system composed of A bosons contained in 
a volume, V. We are interested in getting an infinite system by using the thermodynamic limit, 
i.e. A and V go to infinity keeping the density, p = A/V, constant. We consider a homogeneous 
and translationally invariant system, with a constant density, p. The wave function describing 
the system when the interaction between the particles is switched off is: 

^{xi,X2, . . . ,xa) = S(^(j)i{xi) ■ ■ ■ (pAixA)^ , (4) 

where we have indicated with S the symmetrization operator, with (j)i{xi) the single particle 
wave functions, and with Xi the generalized coordinate of the i-th particle. 

In this IPM description of the ground state of the system, all the bosons occupy the lowest 
single particles state. We consider spin zero bosons, and because of the translational invariance 
of the system, the single particle wave functions are eigenfunctions of the momentum k, and 
they can be expressed as: 

,^,(x,) = -^e^^^-. . (5) 

In this case, the generalized coordinate x corresponds to r. 

The density of the system can be obtained by using Eqs. dH) and ([5]), 

po{x)=AcP*ix)<f>{x) = ^=p , (6) 

which is constant, as expected. 

As we have already discussed in the introduction, the idea is to solve the Schrodinger equation 
by means of the variational principle by using a trial wave function of the form: 

^'(xi, j;^) = ...,xa) , (7) 

where, in this case, the expression of $ is that of Eq. 

For this specific bosonic case we describe the many-body correlation function F(xi, ■■■xa) by 
using the so called Jastrow ansatz [38| [Tl]: 

A 

F{xu....,xa) = n /(^^i) ' (8) 
j>i=i 

where the two-body correlation function (TBCF), f{rij), is a scalar function of the distance 
between the i-th and j-th. particles. 



4 



In the calculation of the energy functional 

it is very useful to employ the two-body distribution function (TBDF) defined as: 

A{A-1) / dx3 . . . dxA'i'*{xi, . . . ,xa)'^{xi, . . . ,xa) 

g{xi,X2) = J — . (10) 

/ dxidx2 . . . dxA'^*{xi, . . . ,XA)'if{xi, . . . ,xa) 

The expectation value of any two-body operator, such as the two-body interaction, is obtained 
by integrating the TBDF on the two coordinates xi and X2- 



< O >=]^p^ j dxidx2gixi,X2)0{xi,X2) . (11) 



The evaluation of the TBDF allows the calculation of the many-body effects independently from 
the explicit expression of the operator. 

By using Eqs. dl]), ([5]), ([6]) and the expressions ([7]) and the numerator and the denomi- 
nator of Eq. (fTO]) can be written respectively as: 



„A-2 r 

M={A-l)^j^ / dxsdxi...dxA Yl finj) , (12 

i<j 



and 



P = ^ y dxidx2...dxA ll f{'r^J) ■ (13) 

i<j 

The cluster expansion is done by defining a new function h{rij) such as: 

f{ri,) = l + h{ri,). (14) 

The product of factors can be rewritten by collecting all the terms with the same number 
of /i-functions. Let's first consider the denominator T>, Eq. (|13p . which can be written as: 

P = / dxidx2...dxA 1^1 + ^ /i(rij) + 3 ^ h{rik)h{rkj) 

i<j i<j<k 

+ Yl HrijMrki) + ■ . .] . (15) 

i<j<k<l 



A convenient way of investigating the structure of the various terms of Eq. (llSp is to 
use the graphical representation introduced by Yvon and Mayer [15]. In this formalism the 
integrated points Xj = rj, which are called internal points, are represented by solid circles, and 
the /i-functions by dashed lines. The expression of Eq. (I15p is obtained by associating to each 
integrated point the contribution of the density. In the present case the density is constant, 
therefore its contributions can be factorized out of the integral. This will not be the case for 
finite systems. 
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I I J I J 

A B C 



Figure 1: Graphical representation of some terms contributing to Ea. (fT5)l . The dashed hne represents 
the correlation function h. The black dots represents the integration points. 

The first sum of Eq. (jl5p is represented by the diagram A of Fig. [TJ The second sum of 
Eq. (jlSp is represented by the diagram B. In this case, the point k is in common with the two 
/i-functions of the sum. The total contribution of this type of term is: 

1(A — 1^ — 2) ^3 J ^^.^^.^^^^^.^^^^^.-^ ^ i^iQ^ 

where the {A — 1){A — 2) factor is due to the fact that the sums on the i,j and k indexes are 
limited to i < j < k. 

In the third sum the two /i-functions involves four different points. Its contribution is repre- 
sented by the diagram C of Fig. [T], and it is given by: 



A\ p4 



4\{A-A)\ A^ 



dxidxjdxkdxih{rij)h{rki) , (17) 



J I 
m •- 



/ 2 1 2 1 2 

ABC 



Figure 2: Graphical representation of some terms contributing to Eq. 
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The use of Eq. (jl4p in the numerator of Eq. (jlOp . allows us to obtain the expression: 



(A-l) 
A 



l + 2^(A-2)^ I dx,h{n,) 



+ 



04-2)04-3) p2 
2 ^2 



E 

j>i>2 



dxidxjh{rij) + 



(18) 



A new symbol is required for the graphical representation of the numerator, since it is necessary 
to indicate the two coordinates which are not integrated. These coordinates are called external 
points, and we have labeled them 1 and 2. The external points are indicated by white circles as 
it is shown in Fig. [21 where we represent the lowest order terms of Eq. (|18p . The uncorrelated 
term is represented by the A diagram. The B diagram represents the terms of the first sum of 
Eq. (jlSp . where the /i-function connects an external and an internal point. Also the second sum 
of Eq. (jlSp contains only a single /i-function but it connects, in this case, only internal points. 
The contribution of this sum is represented by the diagram C of Fig. [2l 

The numerator and the denominator of the TBDF (jlOp are expressed by Eqs. (jlSp and (jlSp 
as sums of terms characterized by the number of the /i-functions, and by that of the external, 
and internal points. Each term of these sums forms a cluster of particles, and can be described 
by a diagram. We proceed now by doing a topological classification of the various diagrams. 

The C diagrams of Fig. [T]can be written by factorizing the non connected terms: 



4!{A-4)! 



J / dridrjdrkdrih{rij)h{rki) 



11(1 



-A + 
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6)p2 



dridrjh{rij) ■ p drkdrih{rki) 



Any diagram that can be factorized in two or more independent pieces is called unlinked. Also 
the C diagram of Fig. [2] is unlinked. The diagrams that cannot be expressed as a product of 
independent parts as the diagram B of Fig. [H are called linked. 



0. 



a •- 



o o 



1 2 1 2 1 2 

ABC 



Figure 3: Examples of reducible diagrams in the numerator of the TBDF, Eq. ([18 



The linked diagrams of Fig. [Sj have the properties of being reducible. From the graphical 
point of view the reducible diagrams are characterized by the presence of, at least, one point 
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linking a part of the diagram containing the external point, and another part containing internal 
points only. Because of the translational invariance of the system, the contributions of these two 
parts can be factorized. In general, every linked diagram whose contribution to the TBDF can 
be expressed as a product of independent integrals is called reducible. In bosonic systems both 
reducible and unlinked diagrams are factorizable. These factorizable diagrams of the numerator 
simplify, up to the 1/A order, all the diagrams of the denominator. The rigorous proof of this 
property is given in Ref. [18j. 

So, in the expression (jlOp of the TBDF, the denominator diagrams compensate the contri- 
bution of the unlinked and of the reducible diagrams of the numerator. Therefore, the TBDF 
can be expressed as the sum of all the irreducible linked diagrams containing the two external 
points 1 and 2: 

giru) = firu) Yl ^^rriru) = f{ri2) (1 + ^(ris) + C{n2)) ■ (19) 

all orders 

The translational invariance of the infinite system makes the TBDF dependent only on the 
relative distance between the external points, ri2. A further topological classification of these 
irreducible diagrams, divides them into simple and composite, and in the above equation we 
have called S{ri2) and C(ri2) the corresponding contributions to the TBDF ()19p . 
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Figure 4: Example of irreducible diagrams, classified as composite and simple. This latter type of 
diagrams is sub-divided into nodal and elementary ones. 

The composite diagrams are those composed by parts that are connected only through the 
two external points 1 and 2, see Fig. HI They can be expressed in terms of simple diagrams. 
Since there is no integration on the external points, the contribution of a composite diagram 
is given by the product of the simple diagrams connected to the external points. For example, 
the contribution of all the composite diagrams which are formed by two simple diagrams only. 
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is S'^{r 12). Since the exchange of all the particles of one subdiagrams with those of the other 
one, produces the same composite diagram, we have to divide by 2 to avoid double counting. 
Repeating the same procedure we find that the contribution of the composite diagrams formed 
by three simple diagrams is 5^(ri2)/3!, and so on. The total sum of composite diagrams can be 
written as: 



The TBDF, ()19p . can be rewritten as: 



1 + S{ri2) + — — + — — + 



g{ri2) = f{ri2) 

= firu) exp[5(ri2)] , (21) 



where the last equality appears because our system has an infinite number of particles and is 
called hypernetted connection. 

The above equation expresses the TBDF in terms of simple diagrams only, which are further 
classified as nodal and elementary ones. In a nodal diagram there is at least one point where all 
the paths going from one external point to the other one have to pass. This point is called a 
node. In the literature, the diagrams without nodes are called elementary or bridge diagrams. 
We shall always use the adjective elementary. Some examples of the type of diagrams we have 
just defined can be found in Fig. HI 

If we call N the contribution of all the nodal diagrams, and E that of the elementary ones, 
we can write the TBDF as: 

g{ri2) = firi2)exp[Nin2) + E{ri2)] (22) 
= [1 + h{ri2)] [1 + iV(ri2) + E{ri2) + . . .] 

= 1 + N{ri2)+X{ri2) . (23) 

The above equation defines the diagrams contained in X{ri2), which are usually named non- 
nodal diagrams since they have no nodes. 



A nodal diagram can be considered composed by parts which are linked at the nodal point. 
Every nodal diagrams can be obtained by integrating on the node the two functions representing 
the parts of the diagram. Let us consider, for example, the nodal diagram of Fig. [5] having i 
and j as external points and k as node. If we call a{rii.) and h{r}^j) the two functions describing 
the two parts of the diagrams, the contribution of this diagram to the TBDF (]19p is: 



drka{rik)b{rkj)p{rk) = p j drka{rik)b{rkj) = (^a{rik) p{rk)b{rkj)^ , (24) 

where a density function p(rfc) has been associated to the integration point in order to recover 
the proper normalization, and since the density is constant in this case, it has been factorized 
out of the integral. The last term defines the symbol ^ ^ we shall use henceforth to indicate 
the folding product or chain connection. 

By using the above considerations we can write a closed expression which allows the evalua- 
tion of all the nodal diagrams. The global contribution N[rij) of all the nodal diagrams between 
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Figure 5: Example of nodal diagram. We indicate with a{rik) the contribution of the part of diagram 
to the left of the nodal point k, and with b{rkj) the contribution of the right part. 



the two points i and j can be obtained as a folding product at the node of all the irreducible 
diagrams which can be constructed between i and k and k and j: 

N{nj) = [x{nk)\p{rk)[N{rkj) + X{rkj)]) . (25) 

Every nodal diagram has at least one node and any path between its external points i and j must 
pass through all the nodes. The above equation tells us that the part of the diagram between 
i and the first node k, which is a non-nodal diagram, has to be folded to: i) the non-nodal 
diagrams producing in this case nodal diagrams with only one node and, ii) the nodal diagrams 
producing nodal diagrams with more than one node. The folding of two nodal diagrams at the 
Yj, point is forbidden since it would produce many times the same diagram. 

The set of Eqs. (f22]l . (f23l) and (f25]l are known as HyperNetted Chain (HNC) equations. 
Equation (j22p allows one to express the TBDF in terms of the simple diagrams after summing 
in a closed form the composite diagrams and ()25p allows the evaluation of the contribution of 
all the nodal diagrams in a closed form. However, there is no closed expression to evaluate the 
contribution of the elementary diagrams which must be calculated one by one. Calculations of 
the TBDM without the contribution of the elementary diagrams are labelled as HNC/0. When 
the contribution of the first elementary diagram is included the calculation of the TBDF is 
called HNC/4, since this diagram, shown in Fig. SI has four particles. These equations are 
usually solved with an iterative procedure starting from the ansatz A^(ri2) = E{ri2) = 0, then 
X{ri2) = /^(ri2) — 1 and we can get new nodals using (125]) . 

2.2 Fermions 

In the description of a system of fermions we have to deal with the Pauli exclusion principle. 
The IPM wave function $ to be used in the trial wave function (El), is now a Slater determinant 
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of single particle wave functions (/>: 



^>(xi, ....,xa) 



(\>\{xa) 



(26) 



For an infinite system we can write the single particle wave functions as: 



1 



V 



(27) 



where we have indicated with s and t the projections on the z axis of the spin and isospin and 
with Xs and xt the Pauli spinors. In the fermions case, the generalized coordinate x indicates 
position r, spin and isospin third components, in addition to the total spin and isospin values. 

Before attacking the problem of the calculation of the TBDF (jlOp we discuss some property 
of which we write as: 



mi,2,...,A)\' 



po{xi,xi) po{xi,X2) ... po{xi,xa) 

Po{x2,Xl) Po{x2,X2) ... Po{x2,Xa) 



po{xA,xi) po{xA,X2) ... po{xa,xa) 
where we have defined the various elements of the above determinant as: 

Po{xi,Xj) =J24>*a{Xi)4>a{Xj) ■ 



(28) 



(29) 



In the above expression the sum runs over all the occupied single particle states of the system. 
We have defined in Eq. (j29p the uncorrelated One-Body Density Matrix (OBDM) which is the 
basic ingredient of the calculation of the TBDF in the fermion case. A fundamental property of 
the uncorrelated OBDM, due to the orthonormality of the single particle wave functions, is: 



dxjpo{xi,Xj)po{xj,Xk) = po{xi,Xk) , 



(30) 



where in the above integral sign we include both the space integration and the sum on the spin 
and isospin third components, their trace. 
We define the sub-determinant as: 



Ap(l,...,p) 



Po{xi,Xi) Po{xi,X2) ... Poixi,Xp) 
Poix2,Xl) Poix2,X2) ... Poix2,Xp) 



p< A . 



(31) 



Po{Xp,Xi) Po{Xp,X2) ... Po{Xp,Xp) 

Because of the property (j30p of the uncorrelated OBDM the sub-determinants have the property: 

(32) 



j dxp+iAp+i{l,....,p + l) = {A- p)Ap{l,...,p) 
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and, by iterating it, we obtain: 

j dxp+i...dxAAA{l,.:;A) = {A-py.Ap{l,...,p) . (33) 

The above expression implies that: 

Ap = , p> A . (34) 

The property (|34p will be extremely useful in the application of the cluster expansion technique 
to both finite and infinite fermion systems. 

The properties of the uncorrelated OBDM and of the sub-determinats we have just presented, 
depend only on the orthonormality of the single particle wave functions, and not on their explicit 
expressions. For this reason, they will remain valid also in the case of finite fermions systems. 
The expression ()27p of the single particle wave functions has been chosen to describe a infinite 
and homogeneous system. In this case, we obtain for the uncorrelated OBDM the expression: 

po{xi,Xj) = -£{kFnj)^xt i^xt {i)Xs{j)xt{j) ■ (35) 

s,t 

In the above equation we have indicated with u the spin-isospin degeneration of the system, 4 
in the nuclear matter case, and with kp = {Gtt'^p/i')^^^ the Fermi momentum. In the literature 
the function i{x) is called Slater function [39], and has the following explicit expression: 

3 

£{x) = ^(sinx — xcosx) . (36) 

In the description of fermion systems, it is necessary to include in the Mayer diagrams a new 
graphical symbol identifying the presence, and the role, of po{xi,Xj), which, in the calculation 
of the TBDF, forms closed non overlapping loops. This is an oriented line connecting the two 
points Xi and xj. These lines are called statistical correlations to distiguish them from the 
dynamical correlations, f{rij). In the calculation of the TBDF for the infinite system a term 
—i{kFrij) jv should be considered for each statistical line joining the i and j points, and a factor 
—V for every closed statistical loop which is related to the spin and isospin trace |35j . 

There is a basic difference between dynamical and statistical correlations. While any number 
of dynamical lines may arrive at a given point only none or two statistical lines may arrive at 
that point. 

By using the trial wave function ([2]) with the Jastrow ansatz ([3]) and the definition (jl4p of 
the /i- function we write the TBDF (jlOp as: 

9{xx,X2) 

A{A-1) dx3....dxA{l + ^hij+ hijhjk + ■■.■)\'^{ xi, xa)\ 

*<i i<j<k 

I dxi....dx^(l + ^/ly + ^ hijhjk + ■■■■)\^{ Xl, XA)\ 

i<j i<j<k 

with hij = h{rij). By using the definition of sub-determinat (j3ip the numerator and the de- 
nominator of the above equation can be expressed as sums of terms identified by the number of 
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/i- functions: 



A{A - 1) . 



/ (^2) / dxz-.-.dxA 



Kj i<j<k 



V 



/ dxi....dxA \'^ + '^hij+ ^ /lij/ijfc + .... Aa 

\ i<j i<j<k J 



(37) 



(38) 



We rewrite the expressions of M and T> by grouping the terms with the same number of points, 
p, and we indicate them as X^P\l,2,3, ..,p). For example 

X(3)(l,2;i) = hu + h2i + hiih2i ■ 

The expression of the TBDF we obtain is: 



gixi,x2) 



^^^^firu) I dx^.-AxA^A 



A 

E(p_2)!(A-p)!- 



1 + 



I dx^...dxAAA^l + Y:^^^^x(^\l,...,p) 



The factorials factors which multiply the X^p^ functions, take into account the fact that permu- 
tations of the p internal points do not change the value of the diagram. 

By using the property ([33]) of the sub-determinants, we can integrate the above expression of 
the TBDM on all the coordinates not involved by the correlations, i. e. not present in the X^^^ 
functions. So we obtain for the numerator and the denominator of the TBDF, the expressions: 



^,/Vi2)^_2_ I dx3...dxj,X^P\l,2;...,p)Ap{l,...,p). 



V = A\Y,-J dxi...dxpX^P\l,...,p)\{l,...,p) 
p=0 ^' •' 



(39) 
(40) 



We extend up to infinity the upper limits of all the sums of the above expression by using 
the property (j34p of the sub-determinants. Each cluster term (diagram) can be divided in linked 
and unlinked parts. Let us call £„(!, 2, 23, i„) the linked parts of the various cluster terms 
containing the external points 1 and 2. In these diag rctins 63.ch. internal point ^3, is connected 
to the points 1 and 2 by at least one continuous path of dynamical and/or statistical correlations. 
We call Up-n{in+i-, ■■■,ip) the unlinked parts of the cluster terms. In this case none of the p — n 
points is connected to 1 and 2, or to another point of Cn- The contribution of does not 
change for a permutation of some of its internal points. The same property holds for Up-n and 
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its internal points. For this reason every diagram of J\f separated in £„ and ^,_n parts, give 
{p — 2)\ / [n — 2)\{p — n)\ times the same contribution. We can then express the numerator if we 
define q = p — n as: 



firu) 



°° 1 f 
■ V- — - / dx3...dxnCnil,2;...,n) 



^ — / dXl...dXqUq{l, 
9=0 ^' 



(41) 



We extend the above considerations to the denominator (I40p . Since in this case there are no 
external points, the diagrams we have defined as linked ones, are not present. Only the lAn 
diagrams contribute to the denominator: 



V = A\ 



°° 1 f 

/ dxi...dxnUn{l,...,n) 

n=0 



(42) 



This expression is identical to that giving the contribution of the unlinked terms of the numer- 
ator. In the calculation of the TBDF the denominator compensates all the unlinked diagrams 
of the numerator, and we can write: 



g{xi,X2) = giru) 



If 

A2(l,2) + J^^— j dx3...dxpCp{l,2;...,p) 



(43) 




Figure 6: Example of cancelation between two FHNC diagrams. The statistical loop in the upper part 
produces a minus sign, and therefore the total contribution is zero. 

The above result shows that the TBDF can be obtained by calculating linked cluster terms, 
only. As in the case of bosons, we define the reducible diagrams as those linked diagrams 
containing a point, the articulation point, which allows us to write the total contribution in two 
or more separated contributions. An example of reducible diagrams is given in Fig. [6j Like in 
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the bosonic case, the factorization of the reducible diagrams in two or more subdiagrams is due 
to the translational invariance of the system. Also in the fermionic case it is possible to show 
that the reducible diagrams do not contribute to the calculation of the TBDF [18j. However, 
in the fermionic case, the mechanism which allows the elimination of the contribution of the 
reducible diagrams, is very different from that of the boson case. Furthermore, the cancelation 
of the reducible diagrams is exact, not limited to 1/A power terms. The rigorous proof of this 
cancellation is given in |18j . and we present here only the basic idea of how the cancellation 
mechanism works. This discussion will become useful to present the vertex corrections in the 
finite fermion systems case. Let us consider, as example, the case of the diagrams shown in the 
upper part of Fig. [H These diagrams differ only because the second diagram has an additional 
statistical loop. Because the system is translationally invariant, and for the properties of the 
Slater function (j36p . the contribution of the the two diagrams is identical but with different 
sign. Therefore, as is shown in the lower part of the figure, the global contribution of the two 
diagrams is zero. 




Figure 7: The various types of nodal diagrams required by the FHNC equations. The subindexes 
classify the diagrams with respect to the type of correlations reaching the external points 1 and 2. 

The fermionic TBDF (j43p can be calculated by considering the contribution of the irreducible 
diagrams only, in analogy to Eq. (|19p for the bosons. Again in strict analogy with what has 
been done for the bosonic case, see Eqs. p9l - [2T]) . it is possible to show that the contribution 
of all the composite diagrams can be obtained by considering simple diagrams only, which are 
classified in nodal and elementary ones. The elementary, and nodal, diagrams in the fermionic 
case are defined in analogy to those of the bosonic case, but both statistical and dynamical 
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correlations should be considered. The presence of the statistical correlations hinders the pos- 
sibility of writing a single integral equation which allows the evaluation in closed form of the 
contribution of all the nodal diagrams, as Eq. (j25p is doing. However, it is possible to find a set 
of integral equations relating the contribution of the nodal diagrams characterized by the type 
of correlations reaching the external points 1 and 2 |18[ [39] . 

Graphical examples of the type of diagrams required to obtain the various integral equations 
of interest are given in Fig. \7\ In the A and B diagrams only dynamical correlations reach 
the external points. These diagrams are labeled with the dd (dynamical-dynamical) subscripts 
i^dd)- The C and D diagrams have only dynamical correlations reaching the external point 1 
and two statistical correlation lines reaching the external point 2. In this case, we label the nodal 
diagram with a de (dynamical-exchange) subscript {N^e)- The E and F diagrams are labelled 
with a ee (exchange-exchange) subscript since to both external points arrive two statistical 
lines. Up to build those ee diagrams with the external points in the same statistical loop, it 
is convenient to define diagrams where a statistical correlation starts from the external point 1 
and arrives to the external point 2, forming an open loop. We label these diagrams with the cc 
(cyclic-cyclic) subscript and we remark that they do not contribute directly to the TBDF. 

As discussed in the bosonic case for Eq. (I25p . also in this case the total contribution of the 
nodal diagrams can be obtained by doing the folding product of various parts of the diagrams at 
the nodal point. However, the Pauli exclusion principle, prohibits some of the possible folding 
products. It is not possible to fold cc diagrams with diagrams of different type and if one of the 
diagrams to fold has the e type at the nodal point the other one has to be d type at this point. 
This restriction is caused by the afore mentioned fact that only two statistical lines may arrive 
at a point, in this case the nodal point. 

In analogy to the bosonic case, we call N the sum of all the nodal diagrams and X the sum of 
all the irreducible non-nodal diagrams. Of course now N and X are classified by the subindexes 
dd, de, ee and cc, and, for the nodal diagrams, we obtain the following set of equations [T8| [39]: 



Nddiru) 



Xdd{ri3) + Xdeiris) p{rs)[Ndd{r32) + ^^^(^32)]) 



+ 




Nde{ri2) 



Xddiris) + Xdeiris) p(r3)[iVde(r32) + ^de(r32)]) 



+ 




Neeiri2) 



Xedins) + Xeeins) p{r3)[Nde{r32) + ^<ie(r32)]) 



+ 




Ncc{ri2) 



Xcc{ri3) p{r3)[N,c{r32) + X, 



cc 



(r32) - iikFr32)/i^]) 



(44) 



The equations for the non-nodal diagrams are: 



Xdd{ri2) 

Xdeiri2) 
Xeeiri2) 



gdd{ri2) - Ndd{ri2) - 1 , 
gdd{ri2)[Xde{ri2) + Ede{ri2)] - Nde{ri2) , 
gdd{ri2){Nee{ri2) + Eeeiru) + [Ndeivu) + Ede{ri2)f 

-i^Wcciru) + ^cc(ri2) - -iikpru)?} - N^eiru) , 
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Xcciru) = 9dd{ri2)[Ncc{ri2) + Ecdru) - -iikpru)] 

- N^ru) + l^ikpni) . (45) 

Finally, the partial TBDF are defined as: 

gdd{ri2) = f{ri2)ey.v[Ndd{ri2) + Edd{ri2)], 
9de{ri2) = Ndeir 12) + Xde{r 12) , 
ged{r 12) = gde{r 12) , 

gee{r 12) = Nee{ri2) + ^ee(?'l2) , 

gcc{ri2) = Ncc{ri2) + Xcc{ri2) - -Kkpru) . (46) 
The total TBDF can be written in terms of the partial ones as: 

9(^12) = gdd{ri2) + ged{ri2) + gdeiru) + geeiru) (47) 




Figure 8: Graphical representation of an elementary diagram. 



The set of Eqs. dM]), (gS]), ([Ml) and (gTD forms the Fermi HyperNetted Chain (FHNC) 
equations. As we have already mentioned in the case of bosons, also the FHNC equations allow 
the evaluation of the contribution of all the composite and nodal diagrams in a closed form. 
Again, the contributions of the elementary diagrams, such as that shown in Fig. [8l should be 
included by calculating them one by one. In analogy to the HNC case, it is common practice to 
use the nomenclature FHNC/0, FHNC/4 etc. to indicate the elementary diagrams included in 
the calculations. 

2.3 The operator dependent correlations 

The use of purely scalar correlations, as it is done in the Jastrow ansatz ([3]), is not adequate 
to deal with the complicated structure of the nucleon-nucleon interaction. For example, the 
tensor terms of the interaction do not contribute to the energy expectation value if only scalar 
correlations are used. For this reason the Jastrow ansatz has been extended by introducing 
correlations which have the same operator structure of the NN interaction, and are called in the 
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literature [21] state-dependent correlations. The general expression of these type of correlations 
is: ^ ^ 

:f{i,...,a)=s[ fl F,,)=s[ fl Y.fp(^^^H) ■ (48) 

j>i=l j>i=lp=l 

The operators Ofj are defined as: 

Off^'^ = l,Ti • Tj,(Ti ■ O-j, {ai • (Tj){Ti ■ Tj),Sij,Sij{Ti ■ Tj) , (49) 

where 

Sij = 3{ai ■ fy)(crj • irij) - en ■ CTj , (50) 

is the tensor operator. The symmetry operator S is required to guarantee the antisymmetrization 
of the wave function ^'(1, ■■■,A) since, in general, the operators do not commute. In Eq. (|49|) 
we have indicated only the channels up to p = 6 since this is the correlation we have used in our 
numerical calculations. State-dependent correlations constructed by considering a larger number 
of channels have been used in nuclear matter [JOj and in variational Monte Carlo calculations 

m- 

The order that we have introduced in the operators will be useful for the finite system since 
we can write: 

^-i+i^pk^^.^^ ^ (51) 

with / = 0, 1, A; = 1, 2, 3 and P^ - = 1, <Tj • cTj, Sij. This allows us to separate clearly the spin and 
isospin parts of the operators. 

The evaluation of the energy functional ([9]) requires the calculation of the expectation value 
of two-body operators related to the NN interaction which are written in terms of the operators 
(|49l) . In general, these operators can be expressed as: 

B{l,...,A)= (E^'MOfi) ' (52) 

j>i=i \p=i J 

and this suggest to define state-dependent TBDFs as: 

A{A - 1) / dx^... dxA^*{l, A)0^2'^(l, ...,A) 

S'p(ri,r2) = -7^ , (53) 

/ dxidx2...dxA^*il,---,A)^{l,...,A) 

where we understand that all the spin and isospin traces are done, including those of the external 
points 1 and 2. With the above definition, the expectation value of B can be calculated as: 

1 ^ r 

<B>=-p^Y. dridr2BPiru)gp{ri,r2) . (54) 

With the help of the sub-determinants (|3ip we express the state-dependent TBDF as: 
5p(ri,r2) = 
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A{A-l) / dx^...dxAS(^ W Fi,)o?2'5( n Fi,^^A{l,---,A) 

(55) 



j>i=l j>i=l 



A A 

j>i=l j>j=l 



In the calculation of the TBDF, we find it convenient to rewrite the correlation function as: 



Because of the non commutativity of the operator dependent terms, in the cluster expansion 
we have to consider also the ordering of the various terms. Only the scalar term p = 1, which 
commutes with all the other ones, can be treated as we have indicated in the previous section. 
By using the commutativity property of the scalar term we can rewrite the correlation function 
as: 

"a 

:f{i,...,a)=s[ n ^*^) = I n 

j>i=l \j>i=l 



j>i=i 



(57) 



This expression shows that each operator dependent term Hij can be multiplied by any con- 
tribution from the central correlation functions, /i, without changing the operator structure 
of the correlation. In the many-body jargon when we incorporate into the operator terms all 
the contributions from the central correlation, we say that the Jastrow correlations dress the 
operator terms. 

The general treatment of the state dependent correlations for nuclear matter was first pro- 
posed in [2i]. This is the basic reference for the interested reader. In the following we shall recall 
the basic steps of the procedure, and we point out the features of interest for the treatment of 
finite nuclear systems. 

Although the notation in the demonstration would be more involved than in the purely Jas- 
trow case, it is still possible to show that the compensation between the unlinked diagrams of the 
numerator and all the diagrams of the denominator holds for the state-dependent correlations. 

We have shown in the previous section that the second step in obtaining the FHNC equations 
consisted in eliminating the contribution of the reducible diagrams. In the present case, this is 
no longer possible. We explain the problem by using the example of Fig. [9] which is analogous to 
Fig. [6] but with the scalar correlations substituted by state dependent correlations, indicated by 
the wavy lines. We should remark that in the graphical representation, the wavy lines indicate a 
generic operator dependent term of the correlation. Diagrams with various wavy lines, consider 
all the possible combinations and ordering of operators. Back to the case of Fig. [9l we should 
consider that all the statistical lines in a closed loop, but one, carry a spin-isospin exchange 
operator: 

1 1 ^ 

W^'^iiJ) ^ -(1 + • a,){l + T, • Tj) = -J2 r^Of, , (58) 

p=i 

on the corresponding pair of particles [21j and with F^ given in Tab. [9]of AppendixjAj The spin 
and isospin dependent part of this operator is linear in cr and/or r, therefore its trace is zero 
when only scalar correlations are used, as in the case of Fig. [6l In that case, the contributions 
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1 2 
A 



1 2 
B 



Figure 9: Diagrams analogous to those of Fig. [5] with state dependent correlations. 

of the two diagrams were identical, with a sign difference, therefore, the global result was zero. 
In the present case, the trace of the operator dependent part of the exchange term is not always 
zero, but its value depends on the type of operators linked to the points at the edges of the 
exchange loop. For this reason, the traces of the A and B diagrams of Fig. [9] are in general 
different, therefore the global result is, in general, different from zero. 

This loss of irreducibility in the expansion of the TBDF joined to the non commutativity 
among the operators makes it not possible to calculate all the contributions of state-dependent 
correlations. This obliges us to use approximations. The difficulty in eliminating the reducible 
diagrams is overcome by using an approximation, consisting in considering operator dependent 
diagrams with specific topological properties. These diagrams are classified as Single Operator 
Chain diagrams. Since the calculations of these diagrams requires the evaluation of traces of 
non commuting operators, we present first the technique used to calculate these traces, and then 
the FHNC equations in the SOC approximation. 

2.3.1 Traces 

In the calculation of the TBDF, or of the energy functional, we have to calculate expectation val- 
ues of operators between the trial wave functions ([2|)- These calculations require the integration 
on the space coordinates and the sum on the spin and isospin coordinates. We have called trace 
this last operation. There are no specific strategies for evaluating the space integrals, which are 
done numerically. On the contrary, there are strategies to obtain the spin and isospin traces. 
Here we shall present them for the specific case of the infinite symmetric nuclear matter, and, 
later, we shall generalize them for finite nuclei. We have emphasized the characteristics of the 
system under discussion, because we want to point out that it is a spin and isospin saturated 
system, therefore the traces of terms linear in cr and/or r operators are zero. 

We start our discussion with the simplest possible case, the matrix element between two 
points only, 1 and 2, having a single operator acting between them. Because of the Pauli 
principle we have direct and exchange terms. The direct term is: 

Of2/'o(l,l)po(2,2) = 
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P'^ E X+(l)xi(2K(l)Xtt(2)0?2Xt.(lk*.(2)x.,(l)Xs.(2) 



Sl,S2,tl,t2 



12/ ' 

where we have used the expression (I35p for po and the fact that the hmit of the Slater function 
()36p when its argument goes to zero, is one. In the above expression we have indicated with 
(7(0^2) the value of spin and isospin trace relative to the operator divided by the number 
of states, v'^ in this case. We call it the C-trace. The exchange term for the case under study is: 

0?2Po(l,2)po(2,l)=p^^^^^ 

E (1)X+ {l)xt (2)0?2Xt. (l)Xt. (2)x.i (1)X.. (2) 

= pH^kpru) C Q(l + cri . o-2)(l + Tl • T2)0?2 

where the term multiplying the C-trace is the spin and isospin exchange operator which 

has to be on the left of the rest of the operators. 

The operators ([19|) are built to be scalar in the Fock space formed by the product of config- 
uration and spin, and isospin, spaces, therefore they are constructed as scalar product of spin, 
coordinates and isospin operators. For this reason, their contributions can be evaluated by using 
the Pauli identity: 

(qi • A)(ai • B) = A • B -FiQi • (A X B) , (59) 

where a = a, t and A and B are generic vector operators. By using this identity, we can 
isolate the terms linear in t or cr which do not contribute in infinite and symmetric nuclear 
matter, as we have already stated. 

The two examples we have discussed are the easieast ones to calculate. In evaluation of the 
TBDF, or of the energy functional, we have to deal with more complicated situations. Following 
Ref. [21] we consider three type of situations. 

a) Products of operators acting on the same pair, such as: 

b) Products of operators acting on different connected points forming a ring, such as: 

'^12 ^23 ^n-\n^n\ ■ 

We call this situation Single-Operator Ring (SOR) 

c) The situation when more than two operators act on an internal point. These are multipole 
operators terms. 

\2. 3.1[ a Products of operators 
We analyze the trace algebra of the products of operators acting on the same pair of 
coordinates 1 and 2. 

The C-trace of a single operator is: 

cm = 6p,i (60) 
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since all operator linear in cr and/or r have zero C-trace in a spin and isospin saturated system, 
as we have already mentioned. 

The C-traces of the product of two operators 012*^0^^^ are obtained by using the relations 

(cti • 0:2)^ = 3 — 2(y.i ■ 0L2 , 

S12CT 1 • <T2 = 512 , 

= 6 + 2<Ti • cr2 - 2^12 , (61) 

calculated by using the Pauli identity ()59p and with a = a^T. The values of the C-traces in 
this case can be summarized as: 

c{0l^0l^=BP5,,, (62) 



where the values of are given in Tab. [10] of Appendix |Al In [21], these are called but we 
shall use this name for their spin parts. 

The knowledge of the values of the C-traces of one operator and of the product of two 
operators is enough to calculate the values of the C-traces for the product of any number of 
operator. The relations ([6T|) indicate that the product of two operators O^^^Oy^^ can be written 
as sum of operators 0^2 multiplied by a coefficient. We can write: 

0?20?2 = E^"''0^2, (63) 

r=l 

where the values of the matrix i^P'?'' are given in Tab. [11] of Appendix [A] For example for three 
operators we have: 

6 

m=l 

We would like to point out that since the operators 0^=^'^ acting on the same points com- 
mute, their ordering does not matter in the calculation of the C-trace. This means that: 



\2. 3. li b Single operator rings 
In Fig. [To] we show an example of Single Operator Ring (SOR). In the SOR diagram each point 
is reached by two operators only. The ordering of these operators does not change the value 
of the C-trace of the SOR. This is because, due to the Pauli identity (I59p . the non commuting 
terms are linear in cr or r, therefore their trace is zero. 

The basic step to evaluate of the C-trace of a SOR it is to calculate the C-trace of two 
operators acting on a common point. We call ^^'^ ^23 operators, and 2 is the 

common point. All the variables relative to the common point 2 should be summed or integrated. 
We sum on all the spin and isospin third components and integrate on the azimuthal angle </>: 

Y: [ 0?2 02^3 = E / ^-^2 01, . (64) 

(T2T2 r=l 
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Figure 10: Example of Single Operator Ring (SOR). 

The ^^^23 functions depend on the angles of the triangle formed by the 1,2 and 3 points, and 
have the following properties: 

^2ki-l+h 2k2-l+l2 2k3-l+l3 _ AfeifefesA A 
?123 ~ ''123 '^h^2"h^3 ' 

+ P2{ri3 ■ r2z){^k^2 - 4i3)'^fc23433 

+ ^^2(^12 • rn)5ki3{5k22 - <^fc23)433 

+ ^2(^^12 • r23)4i3423(25fc32 - 5fc33) 

- 2 (^'^^^^ ■ ^23) (n2 • n3)(?'l2 • ''23) + 1) (5/013423433 , (65) 

with P2{x) = (3x^ — l)/2 the Legendre polynomial of second degree and we have used the sepa- 
ration in spin and isospin parts of the operators presented in Eq. (j5ip . The global contribution 
of the SOR is calculated by using Eqs. and (j65|) for all the points of the ring. 




Figure 11: Example of a multipole operator diagram where more than one operator acts on the same 
two points. 



\2. 3. li e Multiple- operators diagrams 
The last type of diagrams we discuss is represented by the diagram of Fig. [TTJ In this diagram 
two operators act on the same points i and j. The diagram of the figure represents the product 
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of the operators: O^-, Of^ and Of^-. In this case, the value of the C-trace depends on 
the ordering of the operators. It is possible to prove that this value is unchanged by a cyclic 
permutation of the operators of the same type [21j. That is: 



c(ol,{olpi^ )) =c((o^o^ )0- 



As a consequence, for the evaluation of the C-trace we have to consider only two different 
orderings of operators. A first one where Of^ and Ol^ are close to each other, and a second one, 
where these two operators are separated by another operators of the type 01^ or 0\y In the 
first case, by using Eqs. (i63]l and (1651) we obtain: 



t=i ■' 

For the second case we have: 

I d<\>,c(0%0lfllp{^ = ^ I d^^ill] , (67) 

where we have defined: 
The + sign is assigned if 

and the — sign if 



t=i 



c{of^[o!^,ouot^ 
c7(of^.{of^,oijoy = 0, 

where we have indicated with the symbols [, ] and {,} the commutator and anticommutator 
respectively. The values of matrix U"^"^ are given in Tab. [12] of Appendix lAl 
Another important trace is that of 



c oi.oimo. 



q 



which represents two SOR's linked at the point i. Also in this case the result depends on the 



In the first case, we obtain: 



ordering of the operators and we distiguish the case when O^^ and O^^ are close together or not. 



To evaluate the C-trace of the second case we consider the fact that: 

^ 0\piPl = 5,,,B^{1 + E,s)Oi, , (69) 

where, in the case of tensor operators, the above equation assumes an integration over the angle 
between rij and rij [2T]. Then we obtain: 

C(0%0l,0tpi) = 5%,,(1 + Eps)B%,r , (70) 

with the values of E, given in Tab. [13] of Appendix [XI In [21] these are called D but we shall 
use this name for their spin parts. 
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2.4 The Single Operator Chain (SOC) equations 

The strategy to attack the problems arising when state-dependent correlations are used, consists 
in separating the purely scalar, Jastrow, terms from those depending on the operators for 
p > 1. The Jastrow part is treated by using the set of FHNC equations ([11]) , (|15|) and (06]). 
With respect to the operator dependent part, we have learnt that we need at least two operators 
arriving at a given point to get a C-part different from zero. As the operators are dynamical 
correlations, there is no limitation in the number of them that can arrive at every point. An 
increasing number of operators makes more complicated the evaluation of the traces so the Single 
Operator Chain (SOC) approximation is adopted. This supposes that only a pair of operators 
arrive at every internal point of the diagrams, this makes the operators form closed single chains 
and it allows the formulation of closed expressions to calculate all the nodal operators and those 
composite ones within the approximation. The reliability of the SOC approximation is tested 
afterwards, by controlling the validity of sum rules exhaustion. Examples of SOC diagrams, of 
nodal type, are given in Fig. [T2l 

The SOC diagrams do not have limitations in the number of particles. The discussion made 
in Sect. 12.3.11 has clarified that the contribution of the SOC diagrams to the TBDF's and/or 
to the energy, is independent of the ordering of the operators. The single operator between two 
points of the SOC diagram, can come from the correlation, from the hamiltonian, or from an 
exchange line, whose contribution is considered by inserting the spin-isospin exchange operator 




Figure 12: Some nodal diagrams considered in the FHNC/SOC compuational scheme. 



The choice of considering SOC diagrams only, eliminates the problem of the reducible dia- 
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grams, since all the SOC diagrams are irreducible. The equations we should consider contain 
only irreducible diagrams. Also in this case it is possible to show [2l] that the composite di- 
agrams can be obtained in terms of simple diagrams classified again in nodal and elementary 
ones. 

As in the FHNC case, the contribution of a nodal diagram Nmn,pifij), can be obtained by 
doing the folding product of diagrams at the nodal point, as can be deduced by observing the ex- 
amples given in Fig(12j Here we used mn to indicate the type of diagram mn = dd, de, ed, ee, cc. 
As already discussed, the folding product in the point should be done between the irreducible 
non nodal diagrams Xmn,p{rik) and Xmn,q{i'ki) + ^mn,q{fkj)- From the discussion made in the 
previous section, it appears clear that not all the possible combinations of operators acting on 
the integration point provide contributions diff^erent from zero. The allowed combinations are 
ruled by the values of ^^'^'^ given in Eqs. (|65p . 

The sequence of Eqs. ()44p giving the nodal diagrams, is generalized for the state-dependent 
terms as: 



mn,r 



712) 



m' ,n' p,q=l 



X. 



mm' ,pV 13 j ?132 



P{'C?.)[Xn'n,q{r2?,) + N^' n,q{r2?,)] 



(71) 



for m, n = d, e and with m'n' = dd, ed, de. If we neglect the contribution of the elementary 
diagrams, i.e. in approximation FHNC/0, we obtain the following expresssions for the TBDFs 
for p > 1: 



9p{r 12) 
9dd,pir 12) 

9 de,p{r 12) 
9ee,p{r 12) 



9dd,p{ri2) + 2gde,p{ri2) + 9ee,p{ri2) , 
hpir 12) hc{r 12) = Xdd,p{ri2) + Ndd,p{r 12) , 

(hp{n2)Nde{ri2) + fl{ri2)Nde,p{ri2))hc{ri2) 

Xde,p{ri2) + Nde,p{ri2) , 

hpiru) (^Ndeir 12) Nedir 12) + iVee(ri2)) 

+/i (n2) ( - i^L\n2)rP + Nee,p{n2) , 

+2Nde,p{ri2)Ndeiri2))] hc{ri2) 



(72) 
(73) 

(74) 



= Xee,p{ri2) + Nee,p{ri2) , 

where the N and X factor without index of the operator channel indicate the scalar p 
TP is given by Eq. (jSHj) . and we have defined: 

hp{ri2) = 2/p(ri2)/i(ri2) + f!iri2) Ndd,p{n2) , 

hciri2) = eyip[Ndd{ri2)] , 
L{ri2) = N^c{ri2) - l{kFri2)/i^ , 



(75) 
1 term, 

(76) 

(77) 
(78) 



For the calculation of the cyclic nodal diagrams within the SOC approximation, it is con- 
venient to distinguish two cases [21]. Since all statistical lines but one carry one spin-isospin 
exchange operator, we have to add a dynamical operator to complete the operator chain. This 
no exchange operator may be added to the left or to right of the chain. Example of these two 
cases is given by the D and E diagrams of Fig. [T2l In the D diagram the point 1 is connected 
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with an operator dependent correlation while the point 2 is connected by the spin-isospin ex- 
change operator. We label L this case. The situation is reversed in the case of the E diagram, 
and we call R this case. By using this classification we define for the cyclic case the quantities: 



(79) 



Xcc{ri2) = [fKri2)hc{n2) - l\L{n2) , 

NccAn2) = iVc^.(r-12) + A^c^r(n2) , 



(80) 



(81) 




E (^i:p(n3)C?|2n|p(r3)[^cc(r23)+i(r23)]) , 



(82) 



p,g=l 



6 



<r(n2) 



E (Mee(ri3)e?f2|p(r3)[Xi^,(r23) +<(r23)]) 



(83) 



p,q=l 



The set of equations we have presented is called FHNC/SOC. In this case the contribution of 
the elementary diagrams is neglected. Some example of diagrams considered by these equations 
is given in Fig. [T^ The diagram A is a -/V^(i.p(ri2) nodal diagram; the diagrams B and F are 
examples of N(i(,^p{ri2) diagrams and the diagram C of the A'^ee.p(''i2) diagram. Finally, the 
diagrams D and E are of Nccpi^^^) type. 
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3 Finite nuclear systems 



Two of the basic hypotheses done in the previous section, infinite number of particles and 
translational invariance, are no longer valid in the description of finite nuclei. In the literature 
the extension of the FHNC theory to finite nuclear systems was first done by considering systems 
with equal number of protons and neutrons with single particle wave functions produced by 
a unique Mean-Field (MF) potential, within the Is coupling scheme [27l [28l [30l El]. This 
situation allowed a straightforward use of the spin and isospin trace techniques developed to 
describe symmetric nuclear matter [21]. The treatment of nuclei not saturated in isospin and 
described in the jj coupling scheme, was done in following works [29\ I33j. Here we do not follow 
the historical development of the theory but we present directly the formulation of the FHNC 
theory for double closed shell nuclei not saturated in isospin and in the jj coupling scheme. 

The changes of the OBDM expression, due to the loss of the translational invariance, are 
presented in Sect. 13.11 We should point out that we consider doubly magic nuclei only, which are 
spherically symmetric. In Sect. 13.21 we present the calculation of the TBDF and we introduce 
vertex corrections and the Renormalized FHNC (RFHNC) equations, for purely scalar correla- 
tions. The extension of the theory when state dependent correlations are used, is presented in 
Sect. EH 



3.1 The single particle basis 

The nuclear system under study has Z protons, neutrons and, therefore, A = Z + N nucleons. 
The set of single particles wave functions used to describe this system is produced by solving 
the one-body Schrodinger equation: 



(84) 



where the one-body hamiltonian is composed by the kinetic energy term and a spherical mean- 
field potentials different for protons (t = 1/2) and neutrons {t = —1/2): 



h 



t=±l/2 



2mt 



(85) 



where we have indicated with mt the nucleon mass. Our calculations have been done with single 
particle wave functions generated by a Wood-Saxon potential of the form: 



UUr) 



1 + exp 



i?* /a* 



exp 



1 + exp 



(86) 



where is the pion mass and the Coulomb term V^(r), active only for protons, is that produced 
by a homogeneous charge distribution. 



V(j (r) 



(Z-l)eVr 



IRc 



r>Rc 
r<Rc 



(87) 
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The values of the parameters Vq, V^*, Oq, aj^, Rq, and Rc are fixed by the variational 
principle ([T]). In an infinite system, the variational parameter related to the single particle basis 
is the density of the system. 

The eigenfunctions of the hamiltonian (jSSp are also eigenfunction of and jz operators, 
where we have indicated with j the total angular momentum of the single nucleon. The single 
particle wave functions are conveniently expressed as: 

4>iijm{xi) = Rnijin) < ln^s\jm > Yi^{ni)xs{i)xt{i) 

IJ,,S 

= <t>il,„^{r^)xt{i) = Riij{n)Y^{ni)xt{i) , (88) 

where we have indicated with the symbol < | > the Clebsh-Gordan coefficients, and with Yi^j_{^i) 
the spherical harmonics. Here we used the symbol to indicate both polar angles 6i and ipi, 
characterizing the position of the nucleon with respect to a fixed center of coordinates chosen 
to be the center of the spherical nucleus. We have also defined the spin spherical harmonics 

Y-(a) m- 

The uncorrelated OBDM (1291) can be written as: 



Po{xi,Xj) = Y pQ'\ri,Tj)xt{i)Xs'{j)xt{i)Xt{j) 

s,s',t 

where the spatial part is defined as: 



Y,Riij{n)Riij{r,) 



nlj 



< Ip^sljm X Ifi'^Um > Yil{ni)Yi^>{nj) . (90) 

fifi'm 

We find it useful to consider separately the uncorrelated OBDMs of pairs of particles with 
parallel or antiparallel third components of their spins. For these OBDMs we obtain respectively 
the expressions: 

^ ^(2j + l)Rii^{ri)Rii^ir,)Pi{cos9,,) , (91) 



nlj 



1 ':t 



Po ' (ri,rj) = Poj{ri,rj) 

nlj 

where we have called 6ij the angle between and rj, and we have indicated with P/(x) the 
Legendre polynomial of Ith degree, and with P'i{x) its first derivative with respect to x. Some 
useful symmetry properties of these OBDMs are: 

iit — i— it 

Pc)(ri,ri) = Mri,rj) =Po' Mri,rj), (93) 

i -i t -iit 

Poji^i^^j) = Po ' (ri,rj) = -/3o ' ' (ri,rj) , (94) 
Po(ri,ri) = Po(rj,ri) , (95) 
Poj(ri,rj) = -Poj(ri,ri). (96) 
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The uncorrelated OBDM's describing finite nuclei do not depend only on rij = |rj — rj|, 
as in the infinite systems case. However, the properties (I30p . ()32p and (I33p . relevant for the 
construction of the FHNC equations, remain valid. 

The Is coupling can be recovered by switching off the spin-orbit term VJ* = in Eq. (j86p . In 
this case, the single particle energies, ej, and the radial functions do not depend on j, therefore 
the OBDM can be expressed as: 

Po{xi,Xj) = J2Poi^i'^j)xtii)XsU)xti^)XtU) , (97) 

s,t 

with 

p*(ri,r2) = + l)RUri)RUr2)Pi{cose,2) , (98) 

nl 

showing that only the parallel spin OBDM survives. The calculations of Refs. [271 128| [30l 131) 
have been done by using an Is coupling scheme and by assuming equal number of protons and 
neutrons moving in a unique MF potential. With these assumptions the expression of the OBDM 
can be further simplified as: 

p|(ri,r2) = /9o^(ri,r2) = ^^.{21 + 1) Rniiri)Rniir 2) Pi (cos 9 , (99) 

nl 

i.e. the spatial part of the uncorrelated OBDM is also independent of the value of the isospin. 



3.2 The vertex corrections 

The construction of the FHNC equations for the finite systems follows the steps used for the 
infinite systems. The minimization of the energy functional, Eq. ([T|), with the ansatz ([2]) on 
the wave function and ([3]) on the correlation, leads to the requirement of evaluating the TBDF 
(jlOp . The loss of translational invariance obliges us to calculate also the One-Body Distribution 
Function (OBDF), which, in the present case, depends on the isospin third component: 

\ r f 

P*Hri) = ^^^ J dx2...dxA ^*{xi,X2,...,xa)P^ 

X '^{xi,X2,...,xa) , (100) 

where A/j indicates the number of protons (t = 1/2) or of neutrons {t = —1/2), and the projector 
operator P* selects the particle with isospin third component t. By using the above definitions, 
we can express the operator dependent TBDF as: 

^(r„r2) . pr^^Hr„r2) = ^^^f;- ^ ^^^^^^ 

X j dX3...dxA'^*ixi,X2,...,XA)Pl'Ol^P^''^{xuX2,...,XA), (101) 

where the O'^ operators have been defined in Eq. ()49p . In the remaining part of this section we 
shall be concerned only with the calculation of the scalar TBDF, q = 1. The evaluation of the 
other operator dependent TBDFs is treated in Sect. 13. 3i 

The first steps to be done to calculate the one- and two-body density functions defined above, 
are analogous to those used in the infinite systems case. We start by defining an /i-function as in 
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Eq. ()14p and we use it to make the cluster expansion of the numerator and the denominator of 
both OBDF and TBDF. The various terms of the cluster expansions can be analyzed by using 
Mayer diagrams. The topological analysis of these diagrams is done in analogy to what we have 
discussed in the case of infinite systems. 

The arguments used in Sects. 12.11 and 12.21 to show that the contributions of the unlinked 
diagrams of the numerator are simplified by the denominator, can be repeated also in the finite 
systems case [26]. The demonstration is done by formally extending up to infinity all the sums 
of the various cluster terms, since the property (133p of the sub-determinant Ap ensures that 
diagrams containing a number of particles greater than the number of particles forming the 
system, do not contribute. 

In the infinite systems case, the next step was the elimination of the reducible diagrams. We 
have already said that this elimination is only approximated for boson systems, up the 1/A order, 
but it is exact for infinite fermion systems. The basic point of the demonstration for this latter 
case, was the possibility to associate to each reducible diagram, another diagram containing 
only one additional exchange loop. The contributions of these two diagrams to the TBDF differ 
only by a sign, therefore they cancel each other. This cancellation mechanism is produced by 
two specific characteristics of the infinite system. The fact that for a given reducible diagram 
it is always possible to find another diagram having one additional particle, and one additional 
exchange loop, is ensured by the presence of an infinite number of particles. The translational 
invariance is instead responsible for the fact that the additional exchange loop contributes only 
an overall minus sign. In the finite nuclei the number of particles is limited, and the translational 
invariance is lost, therefore there is no cancellation of the reducible diagrams. 

However, even in finite systems it is possible to recover the irreducibility of the expansion 
by introducing the so-called vertex corrections [43^ I26j. A graphical representation of this idea 
is given in Fig. [T3l Every reducible diagram can be thought as composed by two parts, as 
indicated by the diagrams A and B of the figure. A first part contains the external points and 
is the irreducible part of the diagram. A second, reducible, part contains only internal points, 
and it is linked to the irreducible part through the articulation point a. The total contribution 
of these connected, and reducible, diagrams to the TBDF can be written as the folding integral 
of the irreducible part with a function taking into account the contribution of all the diagrams 
connected to the articulation point and it is directly related to the OBDF (jlOOp . 

It is necessary to distinguish the case where the irreducible part is linked to the articulation 
point only by dynamical correlations, as in the A diagram of the figure, from the case when 
there are statistical correlations joining the articulation point, as in the diagram B. To simplify 
the drawing, we show in the A diagram only a single dynamical correlation line connecting the 
irreducible part to the articulation point. In reality, there are no limitations on the number of 
dynamical correlations. In the case of the B diagram we show only the statistical lines connecting 
the articulation point, but also dynamical correlations may be present. 

The fact to be considered is that the Pauli principle allows each point to be reached by 
no more than two exchange lines. When the articulation point is of type d, i.e. linked to 
the irreducible part of the diagram only by dynamical correlations, the Pauli principle is not 
active. In this case, the reducible part of the diagram can reach the articulation point with both 
dynamical and statistical correlations. We call C*(a) = /5*(a) the sum of all the possible linked 
diagrams containing the articulation point a which has isospin third component t. This is really 
the OBDF (jlOOp of the nucleons with isospin third component t. 

The situation changes when the articulation point is of type e, i.e. linked to the irreducible 
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Figure 13: Graphical illustration of the vertex corrections. Since the Pauli exclusion principle allows 
each point to be reached by no more than two exchange lines, we have to distinguish the reducible parts 
of the diagrams. In d the articulation point a can be reached only by dynamical correlations, while in 
C* also by statistical correlations. 
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part of the diagram also by statistical correlations. In this case, because of the Pauli principle, 
the reducible part of the diagram can reach the articulation point exclusively with dynamical 
correlations. We call Cd{a) the sum of the diagrams reaching the articulation point a with 
dynamical correlations only. 

The evaluation of Cd{a) can be performed by extending the diagrams classification done 
in Sect. 12.11 to the case of a single external point. All the linked diagrams, both simple and 
composite ones, contribute to Crf(a). As an example, the contribution of the {II) diagram of 
Fig. [13] is obtained by squaring the contribution of the (/) diagrams and dividing by two, in 
order to avoid double counting. The procedure used in Sect. 12.11 to calculate the contribution 
of the composite diagrams can be applied also in this case. If we call Ud{a) the sum of all the 
simple irreducible diagrams connected to the point a by dynamical correlations only, we can 
write: 

Q(a) = 1 + Udia) + ^Ujia) + ic/|(a) + • • • = exp[Ud{a)] (102) 

It is understood that all the diagrams contained in Ud{a) are renormalized by the vertex correc- 
tions, therefore they must be irreducible in each internal point. 

For the calculation of C*(a) we have to consider also the diagrams linked to the articulation 
point a with statistical correlations. We call U^{a) the sum of all these simple irreducible 
diagrams. Because of the Pauli principle, one can construct composite diagrams with C/*(a) 
combining it only with any number of Ud{a) that produces Cd{a). By definition, C*{a) is given 
by all the diagrams contributing to Cd{a), i.e. all those reaching a with dynamical correlations, 
plus the diagrams constructed by associating those with Ul{a): 

C\a) = Cd{a) [pl{a) + [/*(a)] = p\a), (103) 

where /9o(a) indicates the uncorrelated one-body density for nucleons with isospin i, and p^{a) 
is the corresponding OBDF. In the absence of correlations, Ud and C/* are zero, therefore C* is 
equal to the uncorrelated density, as is expected. The construction of the functions Ud,e{^) is 
done by integrating the composite daigrams over the coordinate 2. This procedure requires a 
careful attention to the avoid the possible overcounting problems [l3]. The explicit expressions 
of Ud,e{X) given in Appendix iBl 

By introducing the vertex corrections idea, we recover the irreducibility of the cluster ex- 
pansion. Obviously the evaluation of the nodal diagrams is more involved than in the case of 
infinite systems. However, the basic ideas used to calculate the nodal diagrams in the infinite 
system, are still valid in the present case, and the expressions of the Ndd{Yi,Y2), A^^g*^(ri, r2), 
A^gi*2 (rj^^ r2) diagrams, given in Appendix |Bl are a rather straightforward extension of those 
presented for the infinite systems. We have shown the dependence on isospin associated to the 
exchange. The only relevant differences are related to the cyclic-cyclic A^*^(ri,r2) diagrams and 
they are worthy of a short discussion. 

In analogy to the infinite system case, the nodal A^*^(ri,r2) diagrams are generated by 
the folding products of X*i(ri,r2) or of pQ^(ri,r2). In the finite system case, the presence of 
the vertex corrections generates the possibility of having nodal diagrams where there are two 
consecutive statistical correlations pQ^(rj,rj). We show in Fig. [14] an example of this situation. 
In the diagram A the point 3 is reached by a statistical correlation on the left hand side, and 
by a dynamical correlation on the right hand side. In the diagram B, the point 4 is reached 
on both sides only by statistical correlations. In infinite systems, i. e. in the absence of vertex 
corrections, because of the property pop of the uncorrelated OBDM, the two diagrams give the 



33 



same contribution, except for a minus sign. For this reason, in order to avoid overcounting, we 
did not consider diagrams of the B type in the evaluation of the N^^ contribution to the TBDF. 
In finite systems each point is vertex corrected, therefore Eq. (|30p cannot be apphed to describe 
the integration over the point 4 in the B diagram, and, consequently, the contribution of the 
diagram B is different from that of the diagram A. 

n 

1 2 
B 

Figure 14: Two A^cc(ri,r2) diagrams. In infinite systems they give opposite contributions since the 
integration over the point 4 in the diagram B recovers the diagram A, except for a minus sign. In finite 
systems, because of the presence of vertex corrections in the 3 and 4 points, the two diagrams give 
different contributions. 

We take care of this fact by separating the contribution of N^^ in four different terms: 

iV*^(ri,r2) = iV,7*Hri,r2)+iV,rKri,r2) 

+ NPf-{r,,r2)+NPP'-ir,,r2) . (104) 

The superscripts x and p refer to the type of correlation reaching the external points. With the 
label p we indicate that the point is reached by a statistical correlation only, like the point 1 in 
Fig. [T4l while with the label x we indicate that also dynamical correlations are present, as in 
the point 2 of the figure. 

The full set of FHNC equations containing the vertex corrections, called Renormalized FHNC 
(RFHNC) equations, is given in Appendix [Bl 

3.3 Treatment of the operator dependent correlations 

In the discussion done so far for finite nuclei, we have considered scalar correlations only. As in 
the infinite system case, the use of the operator dependent correlations is treated in the SOC 
approximation. However, the treatment of these correlations presented in Sect. 12.31 cannot be 
straightforwardly extended to finite nuclear systems. First, the systems we want to describe are 
not saturated in isospin, and this changes the treatment of the isospin dependent terms. Second, 
the jj coupling of the single-particle wave functions modifies the calculation of the spin traces. 
In any case, we tackle the problem by using the strategy outlined in Sect. 12.31 We first reduce 
the product of operators to a sum of single operators, and then we calculate the appropriated 
traces. We find it convenient to treat separately the spin and isospin operators. For this reason, 
we use Eq. (|5T]l . to express the operator dependent correlation (j48ll as: 

6 13 

P^J = E /P(^^^)0f, = E(^^ • E f2k-l+l{U,)Ptj . (105) 

p=l 1=0 k=l 
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As we have already pointed out in the case of infinite systems, we remember that, in addition 
to the product of operators coming from the dynamical correlations and from the interaction, 
we have also to deal with the spin and isospin exchange operator Eq. (j58p . coming from the 
statistical correlation. Each exchange loop formed by n statistical correlations pq{xi,X2) carries 
n — 1 spin-isospin exchange operators. By using the symbols defined in the above equations, we 
express the operators generated by a single exchange loop as: 

p=l 1=0 k=l 

where the values of are those given in Tab. [9] of Appendix [Al from which we see that 
A'^ = 1 - 6k,3- 

3.3.1 Spin traces 

We consider a closed loop of statistical correlations (|89p . involving only two-particles: 

Po{Xi,Xj)po{Xj,Xi) = 

/?o''''*'(ri,rj)x^i(«)Xs;(i)Xtt(Oxti(j) x 

si,s[,ti 

/?o'^'*' (rj , ri)xt2 (^)Xit (•?')^*2 (^) 

^2,52*2 

= Yl 5Z ^o'''*'(ri,rj)Po'''*'(ri>ri) X 

Si,s'^,ti 52,32*2 

xX ij)xt2 i^xt (i)xtt (On'^'^(«> j) Xs[ (i)xs^ (Oxti (j)xt2 (0 , 

where we have used the operator (jSSp . to exchange the spins and isospin of the bra. Fixing the 
isospins, which will be treated appart, the kernel of the above equation is: 

Ypo' '\ri,rj)xt{j)Xs'ij) 

s,s' 

= Po(ri>rj)^X^(i)Xs(j) 

s 

+pl.{r,,r,)Yi-^r-'^\tU)X-s{j) , (107) 

s 

where we used the uncorrelated OBDMs for parallel and antiparallel spins defined in Eqs. 
()91l92p . The treatment of the parallel spin term is similar to that of the infinite system case. The 
antiparallel term, which appears only in the jj coupling scheme, should be treated differently. 

Following the scheme presented in Sect. 12.31 we evaluate the spin traces by considering three 
cases: the product of operators acting on the same nucleonic pair, the product of operator 
forming a ring (SOR) and the product of more than one operator acting on an internal point. 
Consistently with the definitions (15ip of the operators, in the following expressions the upper 
indexes i,j,k can assume the values 1,2 and 3, only. 
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\3. 3.1[ a Parallel spin traces 
In analogy to the infinite system we find that the trace of the product of two operators P*, 
acting on the same pair of nucleons is: 

= B'^^^-'Sk^k, = A'^^Sk.k, , (108) 

where the values of are given in Tab. [TOl Prom this table we see that A'' = k{k + l)/2. 

We find that the product of two operators can be reduced to the following sum of single 
operators: 

Pi2Pi2 = T.l''''Pl2. (109) 

k=l 

where the matrix P^^ is constructed by selecting only the values of the odd indexes of the matrix 
KPir given in Tab. dHof Appendix H 

00\ /010\ /oo 

= ( 3 , P^"^ =1-2 , P^^ = 1 I . (110) 
06/ \002/ \11 

In analogy to Sec. 12.3.11 we can use recursively Eqs. (jlOSp and (jl09p . For example, we find 
for the trace of the product of three operators the expression: 

3 

c(p{2Pl2Pi2) = Yl P^^^C{P^:^P'^^ = P^^A'' . (Ill) 
fci=i 

The evaluation of the product of operators forming a closed loop (SOR) follows the steps 
outlined for the infinite system case. In analogy to Eq. (j64p . we find that also in the present 
case this product of operators can be written as a sum of single operators as: 

3 

j d4>2P{2Pk = E / ^-^2 Cfs PlZ . (112) 

0-2 k=l 

where the values Ci2z Actors are given in Eqs. (j65p . By using the above equation we find that 
the trace of multipole operator diagrams, such as those of Fig. [TTl can be calculated as: 

k-i=i 

The matrix J*-'' is built by using the odd index values of the matrix L given in Tab. [12] of 
Appendix [Al 

\ / 3 

J'Ji = ( 3 , P^^ = 3 6 I , J'J3 = I -6 I . (114) 
6/ \ 

\3. 3. li b Antiparallel statistical function. 
Since the Eqs. ()109p and (jll2p involve only operators, they do not depend upon the spin 
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structure of the wave function, therefore they are vahd also in the antiparallel spin case. The 
change with respect to the parallel case, is in the basic trace value (llOSp which is no longer valid. 
For the antiparallel statistical function of Eq. (jl07p we write the trace of a single operator as: 

CjiP^2) = J E(-l)''^''x^i(l)x^2(2)A'2X-s,(l)x-s.(2) , (115) 
and we find the values 

Cj{Pi\) = dk,2 . (116) 

We see that the contribution of the scalar operator is zero and that of the spin operator is one, 
just the opposite results of those of the parallel case. The result ()116p for the tensor operator 
has been obtained under the hypothesis of a spherical symmetry of the system. 



3.3.2 Isospin expectation values 

In a system not saturated in isospin, we should not sum on the isospin third components, since 
the various diagrams, nodal, elementary, vertex corrections, etc., depend on these quantum 
numbers. This means that the values of these diagrams are different when they are calculated 
for protons, neutrons, or mixed, clusters. We do not calculate isospin traces, but the expectation 
values of products of isospin operators. In order to obtain these expectation values we use the 
properties of the Pauli matrices which allows us to express the product of n isospin operator 
pairs acting on the same pair of nucleons as: 

(Ti • Tj)" = On + (1 - a„)Ti • Tj , (117) 

with: 

ttn+i = 3(1 - an) , and ao = 1 . (118) 

The recursive relation (|117p expresses the product of isospin operator pairs as a sum of a 
scalar term plus a term depending from a single isospin operator pair. The expectation value of 
the operator sequence (I117p is: 

X'^'- = xt (l)Xit (2) (^1 • Xi. (1)X*. (2) ■ (119) 
By using Eqs. (jll7p and (jllSp we obtain: 

X[,^*^ = l, and xf^ = 2<5M2-l, (120) 
and by applying the recursive relations, we have the more general result: 

X'^'^ = 2a„ - 1 + 2(1 - an)bt,t2 . (121) 

This result is the contribution to the cluster expansion of terms like that represented by diagram 
A of Fig. [TSl In this figure we show the various types of diagrams which appear in the calculation 
of the energy expectation value, and more precisely in the calculation of the expectation value 
of the two-body interaction: 

< > < > 



37 



In the figure we indicate with the black area joining the external points 1 and 2, the product of 
the three operators O12O12O12, where and come from the dynamical correlations J-^ while 
O'^, always in the middle, comes from the two-body interaction ¥2- All the isospin expectation 
values necessary for the calculation of the energy expectation value are taken into account by 
Eq. (fml) . 




Figure 15: Example of diagrams contributing to the energy expectation value. The black areas indi- 
cates the product of three operators 0^2*^12 '^12- The two operators , C" come form the dynamical 
correlations, while comes from the interaction. The dashed wiggly line of the diagram B indicate an 
operator dependent vertex correction M . This implies the presence of a single operator coming from the 
correlation, because of the SOC approximation, eventually associated with another operator coming from 
the spin-isospin exchange term. The wiggly dashed and dotted lines of diagram C indicate the correlated 
TBDF containing operators of the correlations and those related to the spin-isospin exchange terms. 



The case of the A diagram of Fig. [15] is an example of our procedure. We find general 
recursive relationships connecting the expectation values of isospin operators, and use them 
to obtain expectation values. The case of the A diagram is the easiest one, but we apply 
an analogous procedure to calculate the cases of the vertex correction, represented by the B 
diagram, and the case of the nodal diagram, represented by the C diagram of the figure. The 
calculation of the isospin expectation values for these two, more involved cases, is presented in 
detail in Appendix ICl 

The rules that we have presented in this section to calculate the spin traces and the isospin 
matrix elements have been used to evaluate the OBDF and the TBDF in the SOC approximation. 
The expressions of the RFHNC equations which allows us to calculate the cluster expansion of 
the TBDF in the SOC approximation are given in Appendix [B] by equations (j2 1411245]) . The 
cluster expansion of the TBDF allows us to calculate the energy expectation value. The details 
of this calculation are given in next section. 
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4 The evaluation of the ground state energy 



The discussion of the previous sections was devoted to the calculation of the scalar TBDF defined 
by Eq. (jlOip when q = 1- The TBDF contains all the many-body effects independent from the 
two-body interaction, whose expectation value can be obtained as indicated by Eq. (fTTI) for the 
infinite system. For the calculation of the TBDF we developed the cluster expansion techniques 
and we built a set of hypernetted chain equations which allows us to evaluate in a closed form 
the contribution of all those diagrams we have called nodal and composite. In this section, we 
apply these techniques to the calculation of the energy expectation value ([9]) with a hamiltonian 
of the form: ^ a A 

^ = -e£v'+ E ^^^■+ E ^^^>^- (123) 

i=l j>i=l k>j>i=l 

The two-body interaction is expressed as: 

y^,=^2'^'ir^3)0';,, (124) 

p=l 

where vP{rij) are scalar functions of the distance rjj of the two interacting nucleons, and the 
p = 1, . . . , 6 operators are those defined in (j49p . We consider in addition the spin-orbit operators 

Of=''^ = (L • S),„ (L • S),,(t, • Tj) , (125) 

where L is the relative angular momentum of the two interacting nucleons, and S is the sum of 
their spins. We give in Sect. 14.31 the explicit expression of the three-body interaction Vijk- 
The calculation of the energy expectation values is done by using the trial wave function: 

^il,...,A)=T{l,...,AMl,...,A)=s(^llFij')'^{l,...,A) , (126) 

i<j 

where the two-body correlation function consider only the first six operators as indicated by 
Eqs. (dSD and (jlUSp . The uncorrelated state, is a Slater determinant composed by all the 
single particle wave functions (j88p lying below Fermi surface. 

In the following, we treat together the kinetic energy and the two-body interaction up to 
the tensor channels. All the other parts of the hamiltonian, two-body spin orbit interactions, 
three-body and Coulomb interactions, are treated separately. In the following expressions the t 
indexes run on protons and neutrons, the p, q, r, s labels may assume values from 1 up to 6, and 
are used to generically identify the different operator channels as in Eq. ()49p . When we separate 
the spin and isospin dependence as in Eq. (fSTI) we shall use the indexes ^=0,1, and k = 1,2, 3. 

4.1 Kinetic energy and part 

We evaluate the expectation value of the kinetic energy, by using the Jackson-Feenberg separa- 
tion scheme [SI US] as suggested in Ref. f39]. We obtain: 

(T) =T^ + Tf- T,,m. , (127) 
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where we have defined 

< ^ > - < ^(Vi$*) • (^^V,^) > I , (128) 

^2 



Am 



i=l i=l 

A A 



Tf = -— < 
4m 



i=l 1=1 

and for the contribution of the center of mass term we have: 



1=1 

In the above equations we have used the symbol <>, which has been defined as |27j : 

< > ■ 

Before attacking the problem of calculating (T) we define some useful quantities: 

P%M) = E K/}n.(ri)V?<^l^,,„(ri) - yi4>%mi^i) ■ ^^i^^Jri)] , (132) 

nljm 

p^*2(ri,r2) =/9*Hri,r2)V?p*^(ri,r2) - Vi/5*Hri,r2) • Vip*^(ri,r2) , (133) 
p%{ri,r2)=2Vlpl'{ri,r2) , (134) 

PT4{ri,r2) = ViPo^(ri,r2) • V2P0 (ri,r2) 

- Po (ri , r2) Vi • V2PQ (ri , r2) 
+ Vi^o^j-(ri,r2) • V2Poi(ri,r2) 

- p*M^i'r2)Vi-V2/9^^,(ri,r2) . (135) 



The above expressions involve neither the interaction nor the correlations, therefore they 
depend only on the uncorrelated many-body state, The expressions of these quantities in 
terms of single particle wave functions (j88p are given in Appendix[Dj In terms of these quantities 
the center of mass contribution can be expressed as: 



Tern = ^ (/ *i/'ti('^i) ~ / dridr2p^r^4^{ri, r2)J . (136) 

ti 



The operator structure of Tp can be easily associated to that required by the calculation 
of the interaction expectation value. For this reason, we calculate together Tp and V2 = 
< v^O^ >. The contribution of Tp + V2 = is also called interaction energy |30] . 

p=l,6 

The structure of the SOC approximation generate various cases that we separate in four 
parts, i.e. : 

W = Wo + Ws + Wc + Wcs . (137) 
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Figure 16: Graphical representation of the W terms of Eq. (|137p . 



We sketch in Fig. [16] the characteristics of the four terms. The black bands between the 
interacting points 1 and 2, indicate O^O'^O^, operators and coming from the correlation 
and O'^ from the interaction. 

In the term Wq we calculate the sum of all the diagrams connected to the interaction points 
by scalar correlation functions only. The interaction point are vertex corrected. If the product 
O^O'^O^ produces a scalar term, the vertex correction contains all the operators types. If, 
instead, the product O^O'^O^ generates operator terms, the vertex correction contains scalar 
terms only. The explicit expression of Wq is given in Appendix [El 

With Ws we consider the sum of operator rings touching a single interaction point. These 
diagrams may include also the presence of scalar operator chains such as those considered in 
Wq. The Wc term is the sum of all the diagrams forming a SOC between the two interacting 
points. We present in Appendix [E[ the derivation of the explicit expressions of the Wg and Wc 
diagrams. 

The contribution of the Wcs term is given by the sum of all the diagrams with operator rings 
reaching a single interaction point, and, in addition, the SOC between the two interacting points. 
The Wcs term is obtained by the combination of the topologies of the Wg and Wc terms. Because 
of the large number of operators present in Wcs, we do not calculate explicitly its contribution, 
but we rather estimate it by using the prescription proposed in [30] : 



Wc^ 



WsWc 

Wo 



(138) 



Nuclear matter calculations [21] where, more refined computational schemes are used, indicate 
that the largest contribution of Wcs is two orders of magnitude smaller than those of the other 
W terms. We have compared nuclear matter estimations of the values of Wcs obtained with our 
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prescription with those calculated more accurately in [2T], and we have found agreement up to 
the second significant figure. 

In the evaluation of the interaction energy W, the Tp part of the kinetic energy is included. 
We describe now how we calculate the contribution of the first term of the kinetic energy (I127P , 
the Tfj, term, where the Vi operator acts on the mean-field wave functions. In Ref. we found 
it convenient to separate the contribution of in three parts: 

r^ = rj^) + rf + rf . (139) 

where each part is characterized by the type of statistical correlations reaching the interacting 
point 1. In Fig. [17] we show some diagrams which identify each term. We have denoted the 
interacting point 1, by an open circle. The nomenclature interacting point is due to the fact that 
this is the point on which the differential operator is acting. In the T^^^ type of diagrams, the 
interacting point is connected to the other points by means of dynamical correlations only. In 

(2) 

, the statistical correlations reaching the point 1, form a closed loop involving only another 

single point. In the diagrams of Fig. [17] we call 2 this other point. Finally, in the type of 
diagrams the statistical correlations reaching the point 1, form close loops which involve at least 
two other internal points. 




Figure 17: Graphical representation of the terms of Eq. (|128p . 

The expression of the T^^^ is the expectation value of a one-body operator dressed by the 
vertex correction, and is given in Appendix lEl 
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The calculation of the other two terms is more involved, since we have to calculate the 
expectation value of a two-body operator. Nevertheless, the structure of these operators is the 

(2) (3) 

same as the ones in the exchange terms of VF. So we further classify the and terms in 
the same way that we did for the W terms. 

(2 3) 

We consider q terms where the points 1 and 2 are connected to the other nucleons by a 
scalar operator chain. As indicated in Fig. [T7]this means that these points are reached by scalar 
correlations only, whose contribution can be calculated in terms of RFHNC diagrams. This is 

(2 3) 

indicated by the fi label in the gray box. A second class of diagrams is that we labeled as ^ . 

(2 3) 

In this case, we consider that, in addition to the scalar chains of T^q , the interaction point is 

(2 3) 

reached by a ring of p > 1 operator dependent dynamical correlations. Finally, with T^^ we 
indicate the sum of the diagrams where the points 1 and 2 form a SOC. In the figure, this is 
indicated by the wiggly lines and by the /p>i label in the gray box. As in the case of W terms, 

(2 3) (2 3) 

we may have also diagrams which are the topological combination of ^ and T^^ . We have 
estimated that their contribution is negligible. 

The separation of the contribution in three parts, was proposed, and used, in [27], where 
the finite nucleon systems treated were saturated in both spin and isospin, and only scalar 
correlations were considered. The presence of operator dependent correlations, requires a further 
classification of the various terms. Clearly, depends on the isospin third components of the 
particles 1 and 2. Furthermore, since we work in a jj coupling scheme, we have to distinguish in 

(2 3) 

the calculation of the ' terms, the cases when the statistical correlations have parallel and 
antiparallel spin components. The complete list of expressions of the various terms composing 
Tfj, is given in Appendix [El 

4.2 Spin-orbit and Coulomb terms 

The contribution of the Coulomb interaction is: 

^ 2 

< vcoul >=< ^* Yl pI'^p]'^—"^ > , (140) 

]>i=l ■' 

1/2 

where the projection operator selects the protons. The Coulomb interaction is added to the 
scalar part of the Vij interaction when two protons interact. This means that its contribution is 
consistently calculated following the methodology described in Sect. l4.1l for all the proton-proton 
W terms. 

We calculate the contribution of the spin-orbit terms of the potential, i.e. the p = 7, 8 chan- 
nels in Eq. (11241) . by considering only diagrams containing scalar chains between the interacting 
points. In other words, for the spin-orbit interaction, we calculate only the Wq term of Eq. 
(jl37[) . The explicit expression of the spin-orbit contribution is given by Eq. (ISlOp in Appendix 

E 



4.3 The three-body potential 

In our calculation we use three-nucleon potentials of Urbana type |46j . The explicit expression 
of this potential is: 



V^Jk = vf^, + V^^, . (141) 
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Figure 18: Fujita-Miyazawa term of the three-body force. 



The vfjj^ term, describes a three-nucleon force produced when one of the interacting nucleons 
is transformed into a A by a first interaction with another nucleon, and it turns back to a nucle- 
onic state after interacting with a third nucleon (see fig. [T8]l . This term called Fujita-Miyazawa 
\A7\ describes the long range part of the three-body interaction and produces an attractive con- 
tribution. The second term of Eq. (jl4ip . vf-^^, has a shorter range and a phenomenological 
origin. The explicit expressions of the two terms of Eq. (jl4ip are: 



^123 



^123 



M-K ^ ({-'^ai, ^32}{t3 • Ti, Ts • T2} 

cycl 



+ ^[-^^31 , -^^^32] [t3 • Tl , T3 • r2 

[/oJ]r2(r3i)r2(r32) . 

cycl 



(142) 
(143) 



where the sums run on all the possible cyclic combinations of the 1,2 and 3 indexes. In the 
above equations we have used the terms T and X defined by: 



3 3 
1 + — + 



[l-e 



k=2 



where: 



fir 



[I 



(144) 
(145) 

(146) 



where, as usual, Sij indicates the tensor operator ([50]) . and the symbols {,} and [, ] indicate 
the anticommutator and commutator operators, respectively. The values of the constants of the 
■^123 term are: /x = 0.7/m~^, and c = 2/m~^. The parameters ^2711 and Uq of the ^(23 term are 
fixed to reproduce the binding energy [46j. 

In analogy to the calculations in nuclear and neutron matter [211 [l6l [2H BOl [25] , we evaluate 
the contribution of the three-body interaction by considering only the sum of the five diagrams 
presented in Fig. [191 In Ref. [l6] it was shown that the diagrams (2.1), (2.2) and (2.3), provide 
the relevant contribution to < ^^23 >• The other two diagrams are those important for < 



123 



>. 



In the diagram (2.1) the pairs of nucleon connected by operators X^ (pairs 31 and 32) are 
dressed by scalar correlations, whereas the remaining pair (12) is also dressed by all the other 
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3 3 3 




(2.1) (2.2) (2.3) 



3 3 




(3.1) (3.2) 



Figure 19: Cluster diagrams considered for the three-body force expectation value. The 2.1, 2.2, 2.3 
diagrams are related to the < vf^j, > part of the force, and the 3.1 and 3.2 diagrams are related to 
< w^j, > . The points denote the particle coordinates. The dashed, wavy, and double-wavy lines represent 
generalized scalar, operator and single-operator ring correlation bonds, respectively. 

operator dependent correlations p > 1, in the SOC approximation. In the (2.2) diagram, the 
pairs 31 and 32 are linked by operator dependent correlations in the SOC approximation, while 
the pair 12 is dressed by the scalar correlation. In the (2.3) diagram there is a cyclic permutation 
of the operator dependent correlations. 

The < Vi2^ > term is calculated by evaluating the (3.1) and (3.2) diagrams. We calculate 
the case where all the pairs are dressed by scalar correlations at all orders, diagram (3.1) and 
the case when there is in addition a SOC correction for a single pair of nucleons, diagram (3.2). 

The detailed derivation of the expressions of the three-body potential is given in Appendix 

m 
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5 Specific applications 



We have applied the formaUsm presented in the previous sections to the description of the C 
, ^^O , ^'^Ca , ^^Ca and 208p|-, nuclei. The only inputs required by our calculations are the two- 
and three-body nuclear interactions. In Sect. 15. 1^ we present those chosen for our studies. The 
single particle wave functions, and the correlation functions, fixed by the minimization proce- 
dure ([I]), are presented in Sects. 15.21 and 15.31 respectively. The theoretical and computational 
reliability of our calculations has been tested by verifying the exhaustion of some sum rules. 
This discussion is developed in Sect. 15.41 After that, we discuss in Sect. 15.51 the results regard- 
ing the most important observable of our calculations: the binding energy. We compare results 
obtained by using two different interactions. We continue our discussion by presenting a set of 
quantities which we have chosen to investigate the effects, and the relevance, of the short-range 
correlations. These quantities are: matter and charge density distributions, Sect. 15.61 momen- 
tum distributions. Sect. 15.71 natural orbits, Sect. 15.81 two-body density distributions, Sect. 15.91 
and, finally, quasi-hole wave functions and spectroscopic factors, Sect. 15.101 

5.1 The nuclear interaction 

The definition of the hamiltonian (I123p requires the definition of both two- and three-body 
forces. We use two-body forces constructed to reproduce the data of the phase-shifts analysis of 
the large body of nucleon-nucleon scattering [3| l48]. 

We have used nucleon-nucleon interactions of the Argonne-Urbana family. These interac- 
tions are local and non-relativistic, and are expressed as a sum of operator dependent terms as 
indicated in Eq. (jl24p . The most recent interaction of this type, fitting the phase-shifts of Refs. 
[HlllH], is the Argonne vis interaction (AV18) [5], containing 18 operator terms, some of them 
breaking the charge symmetry. 

In our calculations we have considered interactions containing up to eight operators channels, 
see Eqs. (I49p and ()125p . For this reason, we used a truncated version of the AV18 potential, 
called Argonne v'g (AV8'), and introduced in Ref. because its simpler parameterization 
allowed a simplification of the numerically involved quantum Monte Carlo calculations. This 
interaction is not a simple truncation of the full AV18 interaction, but its parameters have been 
slightly modified to simulate the effects of the missing channels. The AV8' interaction reproduces 
the results of the full interaction for the S and P scattering waves and also the ^Di wave. The 
details of the construction of the AV8' interaction are given in Ref. [7]. 

The major part of our calculations have been done with the AV8' interaction. However, 
in order to test the sensitivity of our results to the nuclear interaction, we have also used the 
Urbana vi4^ interaction, fixed in Ref. [l9], to reproduce the set of phase shifts data available at 
the beginning of the '80s. In reality we have used only the first eight channels of the interaction, 
without any readjustment of the parameters values. For this reason, we shall refer to this 
interaction as UV14, understanding that we used the operator channels only up to the spin- 
orbit ones. 

The two-nucleon interactions are implemented with the three-nucleon interactions fixed to 
reproduce the binding energy. This means that associated to each two-nucleon potential 
there is a three-body force. Even though more elaborated versions of the three-nucleon forces 
have been recently proposed [50], in our calculations, we use the original expression ^46j, as 
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Figure 20: The AV8' and UV14 nucleon-nucleon interactions as a function of the relative distance, in 
the channels used in our calculations. The continuous thin line, show the behavior of the Argonne vis 
potential. 
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presented in Sect. 14. 3i 

We use the Urbana UIX three-nucleon interaction [7] together with the AV8' potential, and 
the Urbana UVII interaction [5lj together with the UV14 potential. The values of the parameters 
of Eqs. (1142^ . . ■ )145p . fixing these two forces are given in Tab. [H 





UVII 


UIX 




-0.03330 


-0.02930 


Uo 


0.003700 


0.004800 



Table 1: Parameters used in Eqs. (|1421 . . ., I145P to fix the three-body UVII and UIX interactions. In 
both cases fi = 0.7/m~^ and c — 



5.2 The single particle wave functions 





12C 


16q 


4UCa 


4«Ca 






-62.00 


-52.50 


-57.5 


-59.50 


-60.40 




-3.20 


-7.00 


-11.11 


-8.55 


-6.75 


p 

% 


0.57 


0.53 


0.53 


0.53 


0.79 


p 

% 


0.57 


0.53 


0.53 


0.53 


0.79 




2.86 


3.20 


4.10 


4.36 


7.46 


K 


2.86 


3.20 


4.10 


4.36 


7.46 


Rcoul 


2.86 


3.20 


4.10 


4.36 


7.46 




-62.00 


-52.50 


-55.00 


-50.00 


-44.32 




-3.15 


-6.54 


-8.50 


-7.74 


-6.08 


Og 


0.57 


0.53 


0.53 


0.53 


0.66 


n 
"-Is 


0.57 


0.53 


0.53 


0.53 


0.66 


nn 
Kq 


2.86 


3.20 


4.10 


4.36 


7.46 


T>n 
Rls 


2.86 


3.20 


4.10 


4.36 


7.46 



Table 2: The values of the parameters of the Woods-Saxon potential well, Eqs. ((5^ and ([57)1 . The 
superscripts p and n indicate protons and neutrons respectively. The values of Vb and Vi^ are expressed 
in MeV, those of all the other parameters in fm. 



In our calculations, the search of the minimum is done by making variations on the corre- 
lation function and on the mean-field potential generating the set of single particle states. We 
have already said in Sect. 13.11 that we used a mean-field basis generated by two Woods-Saxon 
wells, one for protons and another one for neutrons, both containing spin-orbit terms. The an- 
alytical expressions of the Woods-Saxon wells, Eqs. (j86p and (j87p . involve thirteen parameters. 
Variational calculations done by changing all these parameters would be extremely heavy, from 
the computational point of view. In reality, we found [3T], that the energy minimum is more 
sensitive to the correlation function, than to the single particle basis. To be more precise, we 
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found that when the correlation functions provides an energy minimum, changes of the poten- 
tial do not produce large differences of this value. In Ref. [31] it has been shown that, in ^^O 
, variations of 47% of the central well changed by only 1.2 % the energy value. In calculations 
of the ^^Ca nucleus we found a change of 9% in the energy value by doubling the depths of the 
neutrons and protons wells. 

These findings induced us to make calculations by using, for each nucleus, a fixed set of 
Woods-Saxon parameters. The values of these parameters, given in Tab. [21 are taken from 
the literature [29]. They have been fixed to reproduce the charge root mean square radii and 
the single particle energies around the Fermi surface. One could consider the requirement of 
reproducing these data, a further variational constraint. In any case, we have further verified 
that our energy minima are only slightly modified by large changes of the potential. 

5.3 The correlation functions 

The correlation function is fixed by the minimization procedure ([1]), and the result is independent 
of the starting expression of the correlation function. In practice, however, in order to minimize 
the computational effort, it is convenient to choose expressions of the correlation functions 
containing a limited number of parameters, and behaving at large interparticle distances as 
intuitively expected. This means that asymptotically /i should reach the value 1, while the 
other correlation functions, should be zero. A commonly used expression for the scalar term of 
the correlation is the gaussian form: 

fpir) = (5p,i + ap e"^''^" 

where Up and bp are the free parameters to be changed in the variation procedure. For example 
in Refs. [271 1221 129] correlations of this type have been used. 

In our calculations, we found more convenient, from the physical point of view, and also 
in terms of number of variational parameters, to fix the correlation functions by using what 
we called the Euler procedure. The basic idea is to use as variational parameters the distances 
where the correlation functions fp{r) reach their asymptotic values. For fixed values dp of these 
distances, that we call healing distances, the functions fp{r) are obtained by doing a minimization 
of the energy calculated up to the second order cluster expansion. This means that we solve the 
FHNC/SOC equations when all the nodal diagrams are zero. We give a detailed description of 
the Euler procedure in Appendix iFl 

The application of the Euler procedure involves a single variational parameter, the healing 
distance dp, for each operator channel of the correlation, therefore in our calculations we have 
to deal with six variational parameters. On the other hand, it is know from nuclear and neutron 
matter calculations [2H [Ml [25] that the healing distances of the four central channels p < 4 are 
rather similar, as are those of the two tensor channels p = 5,6. For this reason, we performed 
our calculations by using only two variational parameters, a healing distance dc for the four 
central channels p < ^, and another one dt, for the two tensor channels. 

The values of the healing distances providing the energy minima for the five nuclei considered, 
are given in Tab. [31 In this table, we compare the values obtained by using both interactions. 
In Fig. [21] we show the two-body correlation functions obtained for the AV8'-|-UIX interaction 
as a function of the two-nucleon distance. 
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Figure 21: Two-body correlation functions fp, obtained with the Euler procedure, as a function of the 
two-nucleon distance. In this calculations the AVS'+UIX interaction has been used. 
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it>o 4uca 4«Ca ^^^«Pb 


AV8' + 
UIX 


dc 
dt 


1.20 2.10 2.15 2.10 2.20 
3.30 3.70 3.66 3.70 3.60 


UV14 + 
UVII 


dc 
dt 


1.40 2.10 2.15 2.10 2.20 
3.30 3.80 3.86 3.90 3.80 



Table 3: Values, in fm, of the healing distances, which minimize the energy functional with the two 
different interactions we have adopted. 

Various remarks are in order after observing these results. The most evident one is the scarce 
dependence on the type of nucleus, with the only exception of ^^C . We shall see that, also in 
the calculations of other quantities, this nucleus always produces anomalies in the general trend 
of our results. We think that this is due to the fact that ^^C is not really a good doubly closed 
shell nucleus, therefore, open shell effects, not included in our calculations, are relevant. The 
following observations are done by excluding the ^^C results. 

The results of Tab. [3] show the scarce dependence of the healing distances on the interaction 
model. The values of dc are identical for both the interaction models we have used. There 
are small differences in the dt values, those of the UV14+AVII, are slightly larger than those 
obtained with AV8'+UIX. 

Also the dependence on the nucleus is rather weak. The variations of the healing distances 
are very small, and, the curves of Fig. [21] relative to each nucleus are rather similar for each 
considered channel. 

The values of the scalar correlation functions /i , are one order of magnitude larger than those 
of the other correlations. The behavior of the /i functions reflects the presence of the repulsive 
core in the scalar channel of the nucleon-nucleon force. The correlation function hinders two 
nucleons from approaching each other too much. 

The healing distances of the tensor correlations are larger than those of the central correla- 
tions. Also this effect reflects a characteristic of the interaction where the tensor channels have 
slightly larger interaction range than the central ones (see Fig. [20]) . 

5.4 The sum rules 

The numerical solution of the FHNC/SOC equations is not trivial at all. We have to deal with 
a set of interrelated, hypernetted, integral equations. The numerical technique used to solve 
this set of equation is based on an iterative procedure. In analogy to what we have discussed 
in Sect. [2J] for the HNC equations (|22|23|25p . we started the calculation of the FHNC/SOC 
equations by setting the nodal diagrams to zero. The various integrals equations are calculated 
and they provide new values of the nodal diagrams that are used again to solve the FHNC/SOC 
equations. The convergence test is done on the energy, and we stop the iterative procedure when 
the energy calculated in two different iteration loops differ by less than 1 keV. Every calculation 
is done by using a fixed correlation function. The numerical convergence of the solution does not 
ensure that this solution is acceptable from the physical point of view. For example we found 
numerically convergent solutions which provided a wrong number of nucleons. 

In addition to these computational problems, we should remember that our calculations do 
not solve exactly the many-body Schrodinger equation. The solution of the FHNC/SOC set 
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of equations does not include the elementary diagrams. Furthermore, we do not consider the 
contribution of those operator dependent terms which are beyond the SOC approximation. 

An important tool used to verify the numerical, and theoretical, accuracy of the calculations 
is, the test of the sum rules exhaustion. For this purpose, we have evaluated the following sum 
rules: 

^ drip*Hri) = l, (147) 



'•*^*^(ri,r2) 



1 [ ^ ^ /2fci-l+ii(^^12)/2fc2-l+«2(^12) 

av\aT2- 



{4:l(ri,r2)^^Mfc,..xJ^^,, 

+ [p*2:L(ri,r2)^''^/'=^'=^'^ +P*2:eL,(ri,r2)/'^''^''^/'^^'^' 



/3=0 

S2,a = ^ j dridr2P2'*'*'(ri,r2) 



*1*2=P,"- ' 



1 /" ^^^^^^/2fci-l+«i(^12)/2fc2-l+«2(''l2) 



+ [p2'eL(ri,r2)^''3/^l^^^/^*^^''^ + p'i%,j{YuV2)l''''^^^I^^''''^H^'^^^'^' 

^''ilx\\%,^,,} = -l. (149) 

«3=0 

In the above equations, a sum on the repeated indexes is understood. The operator dependent 
TBDFs P2 are defined in Eq. (llOip . The sum rule 82,(7 is related to the two-body density 
spin function. This sum rule is valid only for spin saturated systems and for correlations not 
containing tensor operator terms. For this reason, we expect 52,0- to be exhausted in ^^O and 
^°Ca nuclei only, and in the absence of the O^^^'^ operator terms in the correlation. 

To obtain the expressions (11471 . • - 1149^ of the sum rules, we have considered all the possible 
types of correlation operator between the 1 and 2 coordinates, which have been vertex corrected 
by the scalar correlations. This is the same approach used to calculate the Wq term of the 
interaction energy Eq. ()137p . 

In Tab. HI we show the sum rule values calculated with the AVS'+UIX interaction, by 
using different types of correlations: purely scalar correlations, /i, central correlation, /4, and 
correlations containing also tensor terms, /g. To simplify the presentation of the results, we 
indicate with 5*2 the quantity: 



^2 = Ail- 1) - 



t2 



(150) 
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1.000 


1.000 


1.000 


1.000 


0.999 




1.000 


1.000 


1.000 


0.999 


0.999 




0.998 


0.996 


0.993 


1.005 


1.008 




0.998 


0.996 


0.993 


1.004 


1.004 




0.997 


1.006 


1.008 


0.994 


1.002 




0.997 


1.006 


1.008 


0.996 


1.000 


S2ifl) 


1.004 


1.003 


1.001 


1.000 


0.998 




0.995 


0.999 


0.989 


1.012 


1.014 




0.996 


0.998 


0.978 


0.994 


1.003 


S2MI) 




-0.95 


-0.929 






S2.{h) 




-1.080 


-1.101 







Table 4: Sum rules exhaustion calculated for the AV8'-|-UIX interaction. The indexes fp indicate the 
number of operator terms of the correlation. 




Figure 22: Elementary diagram included in the RFHNC-1 calculations of Ref. [27] . 

which must be equal to one. The /i results give an indication of the error made by neglecting 
the elementary diagrams. The differences with the other sum rules, is a measure of the validity 
of the SOC approximation. 

Apart from the S2o- values, which we shall discuss separately, the various sum rules are 
satisfied at the level of few parts per thousand. The /i sum rules, give the best results which 
are only slightly spoiled by the other correlations. 

Different remarks should be made for the spin sum rule 82(7- As already pointed out, in our 
calculations, these sum rules have to be satisfied only for the ^^O and ^'^Ca nuclei, and for the 
/i and /4 cases. The values of Tab. [H indicate that the 82(7 sum rule is satisfied only at the 
5-10% level. In Ref. [23, we have verified that the inclusion of the elementary diagram shown 
in Fig. [22] improves the exhaustion of the 520- sum rule up to few parts per thousand. This is 
the same level of accuracy obtained for the other sum rules. 
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5.5 The ground state energies 

The most important results of our work are summarized in Tab. [5]where we give the values of the 
binding energies per nucleon, for all the five nuclei considered. We show the results obtained with 
the two interactions AV8'+UIX and U14+UVII, and we compare them with the experimental 
energies [52]. We present the various terms contributing to the total energy: the kinetic energy 
T, the two-body interaction, where the contribution of the first six channels V^^-body ^^^^ ^^^^ 
of the spin-orbit interaction Vls are separately given, the Coulomb interaction Vcoui and the 
three-body force Vs-body In the kinetic energy term the spurious contribution of the center of 
mass motion, Eq. (jl30p . calculated as discussed in ^7J, has already been subtracted. 

The various terms show some saturation properties. For example, the values of the kinetic 
energies per nucleon, T, increase up to ^''Ca and then they remain almost stable. An analogous 
behavior is shown by the y2-body terms whose contribution per nucleon increases with increasing 
number of nucleons up to ^'^Ca , and afterword it remains almost constant. 

We have mentioned the fact that the spin-orbit terms are not treated consistently in the 
FHNC/SOC computational scheme, but they are evaluated by using some approximation. In 
any case, in all the nuclei considered, their contributions are of the order of a few percent with 
respect to the V2-body contributions. We have done calculations in ^^O and ^''Ca after switching 
off the spin-orbit terms in the mean field potential. In this case the spin-orbit partner single 
particle wave functions are identical. The differences in the total spin-orbit contributions, with 
respect to the values given in Tab. [5] are within the numerical uncertainty. 

As expected, the results of Tab. [5] show that the binding is obtained by a subtle subtrac- 
tion between the repulsive kinetic energy term and the attractive contribution of the two-body 
potential. The sum of only these contributions for the AV8' interaction, provide -2.25, -6.20, 
-8.30, -7.31 and -9.32 MeV for the , ^^O , "^"Ca , ^^Ca and ^ospb nuclei respectively. The 
sum in the U14 model provide -2.40, -6.56, -9.25, -8.19 and -10.17 MeV. It is evident that the 
U14 interaction is more attractive than the AV8'. This depends on the intrinsic structure of the 
interaction and its parametrization. 

The contribution of the Coulomb term VqouI is evaluated within the complete FHNC/SOC 
computational scheme. As expected, the behavior with increasing size of the nucleus does not 
show saturation because of the long range nature of the interaction. The Coulomb terms behave 
as expected, their contributions increase with increasing number of protons. The apparent 
inversion of this trend from ^'^Ca to ^^Ca is due to the representation in terms of energy per 
nucleon, which in this case is misleading, since the proton number is the same for the two nuclei. 
In this case it is better to compare the total values of the Coulomb energies, 78.80 MeV for 
^'^Ca and 75.36 for ^®Ca . The 4.4% difference between these two values is due to the different 
structure of the two nuclei. The inclusion of the Coulomb repulsion reduces the nuclear binding 
energies. 

In addition, there is the contribution of the three-body force. As discussed in sect. 14.31 
the two terms composing this interactions provide contributions of different sign; the Fujita- 
Miyazawa term w^23 is attractive, while the other term ^(23 is repulsive. In our calculations, 
the total contribution of the UVII and UIX three-body interactions is always globally repulsive. 
This feature is common to the FHNC/SOC nuclear matter results \24:\ |25J. 

The comparison with the experimental energies indicates a general underbinding of about 
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AV8'+UIX 








4«Ca 


2U8p]^ 


T 


27.13 


32.33 


41.06 


39.64 


39.56 


^2-body 


-29.13 


-38.15 


-48.97 


-46.60 


-48.43 


Vcoul 


0.67 


0.86 


1.97 


1.57 


3.97 


Vls 


-0.25 


-0.38 


-0.39 


-0.35 


-0.45 


T + V2 


-1.58 


-5.34 


-6.34 


-5.74 


-5.35 


y^i-body 


0.67 


0.86 


1.76 


1.61 


1.91 


E 


-0.91 


-4.48 


-4.58 


-4.14 


-3.43 


U14+UVII 


i^C 
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2U8p]^ 


T 


24.63 


29.25 


37.32 


36.12 


36.07 




-27.08 


-35.84 


-46.65 


-44.40 


-46.28 


Vcoul 


0.68 


0.88 


2.01 


1.59 


4.00 


Vls 


0.05 


0.03 


0.08 


0.09 


0.04 


T + V2 


-1.72 


-5.68 


-7.24 


-6.71 


-6.17 


y 3 -body 


0.54 


0.69 


1.46 


1.32 


1.61 


E 


-1.18 


-4.99 


-5.77 


-5.27 


-4.55 


E 

^exp 


-7.68 


-7.97 


-8.55 


-8.66 


-7.86 



Table 5: Energies per nucleon in MeV, obtained by using the AVW + U IX and UlA + UVII interactions. 
We have indicated with T the kinetic energy, with V2_},ody contribution of the first six channels of the 
two-body interaction, with Vls the spin-orbit contribution, with Vcoui the contribution of the Coulomb 
interaction and with Vs-tody the total contribution of the three-body force. The rows labeled T + V2 
show the energies obtained by considering the two-body interactions only. The experimental energies are 
fromRef. IMl. 



4.0 MeV per nucleon. This is roughly the same underbinding obtained, at the saturation density, 
by the most recent FHNC/SOC nuclear matter calculations [25j. 

The behavior of the ^^C nucleus is anomalous in this general trend. This nucleus is barely 
bound in our calculations. Some crucial physics ingredient, relevant in ^^C , but negligible for 
the other nuclei, is missing in our approach. Probably, this has to do with soft deformations of 
the ^^C nucleus, effects which we are unable to treat. 

The comparison between the two interactions indicates that the UV14-I-UVII interaction is 
more attractive than the AV8'-|-UIX force. This fact is already present when only the two-body 
interactions are considered, and it is enhanced by the inclusion of the three-body force. The 
contributions of the spin-orbit term in the two cases have different sign, they are attractive for 
AV8' and slightly repulsive for UV14. Globally, the differences in the total energies, calculated 
with the two interactions, vary from a minimum of 5% (^^O ) to a maximum of 18% (^'^^Pb ). 

We have done a detailed study of the relevance of the various terms contributing to the 
energy, as they have been presented in Sect. [H As an example, we show in Tab. [6] the various 
contributions obtained in the calculation of the ^osp^^ g^ergy with the AV8'-|-UIX interaction. 
We have obtained analogous results for all the other nuclei investigated and also for the other 
interaction. 

The larger contributions to the energy come from the terms calculated with what we have 
called the (0) approximation in Eq. (jl37p . This is the contribution of those diagrams containing 
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-1.49 
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3.282 
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-1.368 



Table 6: Contributions, in MeV, of the various terms forming the ^'^'^Pb energy, calculated with the 
AV8'+UIX interaction. The various terms are defined in Sect. [H The terms T of the kinetic energy are 
defined by Eqs. (|127[ 11281 11291 fT39)) . The terms indicate the six-channels of the two-body interaction. 
The three-body terms are defined in Eqs. (|142|143|) . The subscript j indicate the contribution produced 
by antiparallel spin densities. The labels (0), (s) and (c) indicate the various approximation of the energy 
related to the pieces Wq, Wg, and Wc as defined in Eq. (|137p . 



all the scalar dressings of the (1, 2) pair. The (s) and (c) terms, more difficult to calculate, give 
much smaller contributions. 

It is interesting to observe that, the contributions given by the terms depending on one-body 
densities with antiparallel spins, i. e. the terms labeled with the j subscript in the table, are 
very small. In Tab. [6] we show their largest values, since we found their contributions to be even 
smaller for nuclei with saturated / shells such as ^^O and ^''Ca . 

The study of the contributions of the various channels of the two-body interaction, indicates 
that the spin-isospin [p = 4) and tensor-isospin {p = 6) terms are the main source for the 
binding. This is a common feature for all the nuclei we have considered [53]. 

The results for the UV14-I-UVII interaction are analogous to those shown in Tab. [6] for the 
AV8'-|-UIX interaction. In the remaining part of Sect. [5l we shall present some quantities with 
the aim of studying the effects of the Short-Range Correlations (SRC). We have found that the 
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Figure 23: Neutron density distributions. The full lines are the IPM distributions, the dotted ones have 
been obtained by using scalar correlations only, and the dashed lines show the results of the complete 
calculation. 

differences between the results obtained with the two interaction models are smaller than the 
effects we are looking for. For this reason, henceforth, we shall present only the results obtained 
by using the AV8'+UIX interaction. 

5.6 The one-body distribution functions 

The OBDF, /9*(r), has been defined by Eq. (jlOOp . The physical meaning of this quantity is the 
probability density of finding a nucleon of type t, in the position r with respect to the nuclear 
center. Since we have assumed that our systems are spherical, this probability depends only on 
the distance from the center of the nucleus. The expression of the density distribution in terms 
of the FHNC/SOC quantities is given by Eqs. (fT03D and (|239]1 . 

We show in Fig. [23] the neutron density distributions for the five nuclei we have considered. 
The full lines show the IPM results, the dotted lines those obtained by using scalar, /i, cor- 
relations only, and the dashed lines show the results of the full FHNC/SOC calculations. In 
an analogous way we show in Fig. [25] the charge distributions, obtained by folding the proton 
distributions with the proton electromagnetic form factor. We have used a dipole form for this 
form factor. In this figure the continous thick lines indicate the empirical charge distributions 
extracted from elastic electron scattering experiments |54j . 

The results obtained with scalar interactions produce distributions which are smaller at the 
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Figure 24: Charge density distributions, obtained by folding the proton distributions with the electro- 
magnetic form factor. The full thick lines show the empirical distributions [54 . The other curves have 
the same meaning as in Fig. 1231 
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Figure 25: Elastic electron scattering cross sections on "'^Ca and ^°*Pb target nuclei, as a function 
of the scattering angle. The theoretical cross sections have been calculated in distorted wave Born 
approximation [SHI [57]. The full lines show the IPM results, the dotted lines the results obtained with 
the /i correlations, the dashed lines with the /e correlations. We have indicated with e the electron 
energy. 

center of the nucleus with respect to the mean-field distributions. This effect is reduced when 
all the correlations are included in the calculation. These findings are in agreement with the 
results of Ref. [55] where a first-order cluster expansion was used. 

We used the charge density distributions to calculate the elastic electron scattering cross 
sections within a distorted wave Born approximation [561 [57]. We compare in Fig. [25] our 
results with the experimental data of the ^^Ca and 208p^ nuclei |58[ 159] . In Ref. [31j a similar 
figure for the ^^O and ^''Ca nuclei is shown. It is evident that the main discrepancies with 
the data appear at large scattering angles. The use of the correlations slightly improves tha 
agreement with the experiment. The effects of the /e correlations go in the opposite direction 
of those of the fi correlation (as has already been remarked) in the neutron density, and in the 
charge distributions cases. 

We would like to point out here that the density distributions are the only quantities amongst 
those we have investigated where the use of operator dependent correlations reduces the effects 
of purely scalar correlations. In all the other cases, as we are going to show, the effects of the 
/e correlations are larger than those of the /i correlation. 
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5.7 Momentum distributions 



The momentum distribution is related to the probabihty of finding a nucleon with a certain 
value of the momentum. This quantity is related to the Fourier transform of the OBDM, which 
is defined in analogy to the OBDF (jlOOp but for different values of the non integrated variable: 



p{xi,x[) 
A 



s,s't 



dXA *^(X1,X2,... ,XA)'^ix'i,X2, . . . ,xa) 



(151) 



The evaluation of this quantity merits some comments. For simplicity, we consider in this 
discussion only scalar correlations. In order to perform the cluster expansion of the OBDM, it 
is necessary to define a new dynamical correlation function: 



hwinj) = fi{rij) - 1 



(152) 



where i can be either 1 or 1', and j = 2, . . . ,A. This takes into account the fact that the 
coordinate xi is present only in the bra, while the coordinate x'^ is present only in the ket. The 
coordinates describing the other particles can appear in both bra and ket states, and generate 
the usual h = — 1 dynamical correlation function. For the calculation of Eq. (jl5ip . we found 
it convenient to define a new type of sub-determinant, in analogy to what has been done in Eq. 
(El: 

Po{xi,Xi') Po{xi,X2) ... po{xi,Xp) 
Poix2,Xl') Pq{x2,X2) ... Poix2,Xp) 



a;(i,i', 



(153) 



Po{Xp,Xi>) Po{Xp,X2) ... po{Xp,Xp) 

Also for this new sub-determinant the properties (1321 . . [Ml) remain valid, therefore we can apply 
the usual cluster expansion techniques, developed for the calculation of the TBDF. However, 
we have to consider that the separation of the xi and xy coordinates which refer to the same 
particle, implies the absence of dynamical correlations between these two coordinates, as is shown 
in the diagrams of Fig. [26j In addition, we should take care of the fact that the statistical loops 
containing the coordinate xi must contain also xy and must be open in these two points, see 
again the diagrams of Fig. [26l 



The calculation of the OBDM proceeds in analogy to that of the TBDF. The denominator of 
Eq. ()15ip simplifies the contribution of the unlinked diagrams of the numerator. At this point, 
we encounter a difficulty, since it is not possible to cancel those reducible diagrams where the 
articulation point is one of the points in the open statistical loop containing to 1 and 1'. This 
is not a specific problem of the finite systems. Its solution is based on the use of the vertex 
corrections which have been introduced for the first time to calculate the momentum distribution 
of an infinite system of particles [l3] . 

Coming back to our specific case, we found that the calculation of the OBDM requires the 
use of a new vertex correction, which we call Cu;(rj), understanding that i can be 1 or 1'. New 
types of diagrams appear. We call them wd, we, ww, WcC, and WcWc- The label w is associated 
to the new dynamical correlation /t^ and, as d, denotes that no statistical line arrive at the 
corresponding point. With the label Wc we indicate diagrams analogous to those we have called 
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Figure 26: Example of diagrams used to calculate the OBDM (|15ip . The dotted lines indicate the new 
dynamical correlation function = /i — 1, while the dashed line indicate the usual dynamical correlation 
function h = — 1. 

c in the calculation of the TBDF. These diagrams have an open statistical loop with the 
dynamical correlations, as the diagrams of Fig. [26l The diagrams contributing to the OBDM 
are of WcWc type and the rest of classes diagrams are auxiliary quantities needed to calculate 
them. The RFHNC/SOC expressions of the different quantities involved in the calculation of 
the OBDM, are presented in Appendix |B] where the parallelisms with the TBDF expressions are 
also pointed out. 

Because of the spherical symmetry of the systems we are describing, the quantity of interest 
in our calculations is: 

p\r,,r[)= [p'''-'\rur[) + p^'-^-'\r,,r[)] , (154) 

s=±l/2 

whose diagonal part, = ri, is the OBDF. We obtain the momentum distributions of protons 
or neutrons as: 

""^^^^ = ji^Jft I ^rirfr;e^'^-(r-r'i)/(ri,r;) , (155) 

which is normalized: 

j dkn\k) = 1 . (156) 

The uncorrelated OBDMs are obtained by inserting in Eq. (j89p the Slater determinant, 
formed by the single particle wave functions ([55]) . We obtain the expressions (pi]) for s' = s and 
(pj) for s' / s. 

For the correlated OBDM we obtain the expression: 

p\r,,r[) = -2C4,ii(ri)C^,ii(r'i)5Uc(ri,r'i) 

-2C4,22(ri)C^,22(ri')E^'^'5U.,p(ri,r'i) • (157) 
P>1 

with p = 2k — 1 + I. All the FHNC/SOC quantities have been defined in Appendix iBl 

In Fig. [27] we compare the ^^C , ^^O , ^''Ca , ^^Ca and ^osp]^ momentum distributions 
calculated in the IPM model, with those obtained by using fi and /g correlations. 
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Figure 27: The proton momentum distributions of the ^^C , ^^O , '^''Ca , ^^Ca and ^°*Pb nuclei calculated 
in the IPM model (full lines), by using the scalar correlation only (dotted lines) and the full operator 
dependent correlations (dashed lines). 
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Figure 28: Protons (full lines) and neutrons (dashed lines) momentum distributions of the ''^Ca and 
208p|-, 'Yhc thick lines show the results of our calculations, the thin lines the IPM results. 



The general behavior of the momentum distributions, is very similar for all the nuclei we have 
considered. Correlated and IPM distributions almost coincide in the low momentum region up to 
a precise value, when they start to deviate. The correlated distributions show high momentum 
tails, which are orders of magnitudes larger than the IPM results. The value of k of which 
uncorrelated and correlated momentum distributions start to deviate is smaller the heavier is 
the nucleus. It is about 1.9 fm~^ for ^^C , and 1.5 fm~^ for ^''^Pb . We recall that the value 
of the Fermi momentum of symmetric nuclear matter at the saturation point is 1.36 fm^^. The 
results presented in Fig. [27] clearly show that the effects of the scalar correlations are smaller 
than those obtained by including the operator dependent terms. 

In our calculations, we have found that the proton and neutron momentum distributions 
for nuclei with N = Z are very similar. For this reason, we show in Fig. [28] the proton and 
neutron momentum distributions of the two nuclei we have investigated with N ^ Z: the ^^Ca 
and ■^''^Pb nuclei. The thicker lines show the results of our RFHNC/SOC calculations, while the 
thinner ones the IPM distributions. The main differences between the two distributions are in 
the zone where the n{k) values drop by orders of magnitudes. This zone, corresponds, in the 
infinite systems, to the discontinuity region of the momentum distribution, related to the Fermi 
momentum. In a finite system, the larger number of neutrons implies that the neutron Fermi 
energy and, consequently, the effective Fermi momentum, is larger than that of the protons. 
For this reason, the discontinuity regions of the neutron momentum distributions are located at 
larger values of k with respect to the protons momentum distributions. 

After the discontinuity region, the behaviors of the distributions are dominated by the SRC 
effects, and the protons and neutrons results are very close. In terms of relative difference, the 
SRC effects are essentially the same for protons and neutrons |60j . A discussion about the role of 
SRC effects on the proton and neutron momentum distributions in asymmetric nuclear matter 
is open, and our results are in agreement with the findings of Ref. ^61j) but disagree with those 
ofRef. [62]. 

The increase of the momentum distribution at large k values, induced by the SRC, is a well 
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Figure 29: Occupation numbers of the proton natural orbits of the ''^Ca nucleus. The dashed line 
indicates the IPM values. The black bars show the values obtained with the scalar corrrelation and the 
gray bars those values obtained with the full correlation. 



known result in the literature, see for example the review of Ref. [63] . The momentum distribu- 
tions of medium-heavy nuclei, have been usually obtained by using approximated descriptions 
of the cluster expansion, which is instead considered at all orders in our treatment. We found 
in [60j that the results of the approximate treatment provide only a qualitative description of 
the correlation effects. They produce a high- momentum enhancement, but they underestimate 
the correct results by orders of magnitude. 

5.8 Natural orbits 

We have studied the effects of the SRC on the natural orbits which are defined as those single 
particle wave functions forming the basis where the OBDM is diagonal: 

p\r,,r[) = • (158) 

nlj 

In the above equation the c^;^- coefficients, called occupation numbers, are real numbers. In 
the IPM the natural orbits correspond to the mean- field wave functions of Eq. (|88|) . and the 
c^^ij numbers are 1, for the states below the Fermi surface and for those above it. 

In order to obtain the natural orbits we found it convenient to express the OBDM of Eq. 
(fT58]) as: 

p\ri,rr) = A\ri,rr)pi{ri,ry) + B\ri,rr) , (159) 

where /9Q(ri, ri/) is the uncorrelated OBDM of Eq. (|9ip . and the other two quantities are defined 
as: 

^*(ri,r;) = 2C4,n(ri)C^,ii(ri)5"^„(ri,r;) + 
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Figure 30: The same as Fig. [55] for tlie neutron natural orbits of the ''^Ca nucleus. 



P>1 

i?*(ri,r;) = -2C^,n(ri)C4,iiW)5"«,(ri,r;X^^^(ri,r;)- 

2C4,22(ri)C^,22«)E^'^'[^^--,p(ri'r'i)<c-.(^^ 
P>1 

5L(ri,r'iX^^^_p(ri,r'i)j . 



(160) 



The expressions of the various FHNC/SOC quantities used in the above equations are given in 
Appendix [BJ 

We expand the OBDM on a basis of spin spherical harmonics YJj defined in Eq. (jSSp . 

p\r^,r[) = [4i(n,rl) +4(ri,rl)] Y^pimi^m (161) 

Ijm 

where Q and 0,' indicate the polar angles identifying ri and r'^. The explicit expressions of the 
A and B coefficients are: 



A%{r,,r[) = (2/ + 1) J] (2;2 + l)(2i2 + l)f Q ^ 



with 



and 



21 + 1 



dnA\Yi,r[)Pi{cos9 



11' 



An 
21 + 1 



d(cos eir)B\ri,r[)Pi{cose 



11' 



(162) 
(163) 
(164) 
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Figure 31: Natural orbits for some neutron states in ^*Ca . The full lines indicate the IPM orbits, the 
dotted lines those obtained with scalar correlations only and the dashed lines those obtained with the 
complete operator dependent correlation. 
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Table 7: Protons (p) and neutrons (n) natural orbits occupation numbers for ^^O . The PMD values 
are those of Ref. [SS] . The /i values have been obtained with scalar correlations only, and the /g values 
with the complete state dependent corrrelation function. 

In the above equations we have used the 3j and 6j Wigner symbols [42j and we have indicated 
with Oil' the angle between ri and r'^,. The term A depends on both orbital and total angular 
momenta of the single particle, / and j respectively, and B depends only on the orbital angular 
momentum /. 

As has been done in Refs. [64j and [65] we identify the various natural orbits with a number, 
a, ordering them with respect to the decreasing value of the occupation probability. The general 
behavior of our results, is analogous to that described in Ref. [64] where a system of ^He drops, 
composed by 70 atoms, have been studied. The orbits corresponding to states below the Fermi 
level in the IPM picture, have occupation numbers very close to unity for a = 1, and very small 
in all the other cases. 

As example of our results, we show in Figs. [29] and [30] the protons and neutrons occupa- 
tion numbers for the natural orbits with a = 1 of the ^^Ca nucleus. In the figures, the IPM 
results are indicated by the dashed lines. The black bars show the values obtained by using 
scalar correlations only, the gray bars those obtained with the complete operator dependent 
correlations. 

The correlated occupation numbers are smaller than one for orbits below the Fermi surface, 
and larger than zero for those orbits above the Fermi surface. This effect is enhanced by the 
operator dependent correlations. We observe that for the states above the Fermi surface the gray 
bars are larger than the black ones, indicating that also for these states the operator dependent 
correlations, produce larger effects than the scalar ones. 

We show in Fig. [31] some a = 1 natural orbits for three neutron states in ^'^Ca . In this 
figure, we compare the IPM results (full lines) with those obtained with scalar correlation only 
(dotted lines), and with the full operator dependent correlation (dashed lines). The effect of 
the correlations is a lowering of the peak and a small widening of the function. Despite the 
small effect, it is interesting to notice the inclusion of operator dependent terms diminishes the 
correlation effect. This fact is consistent with the results on the density distributions we have 
presented in Sect. 15.61 

In Tab. [7] we show the occupation numbers of the ^^O protons and neutrons natural orbits 
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also for a > 1, and we make a direct comparison with the results of Ref. [65]. As already said 
in the discussion of the ^^Ca results, the inclusion of the state dependent correlations increases 
the differences with respect to the IPM. The occupation numbers of the orbits below the Fermi 
surface are smaller than those obtained with scalar correlations only. The situation is reversed 
for the orbits with a > 1 or above the Fermi level. For the states below the Fermi surface, 
our full calculations produce correlation effects slightly larger than those found in [65], whose 
results are closer to those we obtain with scalar correlations only. For orbits above the IPM 
Fermi surface, our occupation numbers are always smaller than those of Ref. [65] . 

5.9 The two-body distribution functions 



We have already defined the state dependent TBDF in Eq. (jlOlh . In our FHNC/SOC compu- 
tational scheme, we calculate the TBDF as: 



P2 ''''(ri,r2) 



f2ki-i+h {ri2) f2k2-i+l2 (^12) 
fKri2) 

il''''''^''^l%2+lsPttiri,r2) 



+ 



Jfc4fclfc5jfc3fc2fc5^fc5^ili2^(ri,r2) 



14=0 



(165) 



where a sum is understood on every repeated index. The FHNC/SOC quantities are defined in 
Appendix [BJ 

We first discuss the case q = 1 which, apart from some constant factors, gives the probability 
of finding a nucleon in ri, and at the same time, another nucleon in r2- The IPM two-body 
densities are obtained by inserting a Slater determinant in Eq. (jl67p . For q = 1 TBDF we 
obtain the expression: 



+ 



P0j{T^l,T^2] 



(166) 



where the /9*^^_^.^(ri, r2) are the uncorrelated OBDMs defined in Eqs. (|9ip and (j92p . 
The relevant information about the TBDF is contained in the function: 



pf'''in2) = I dRi2pf^*^(ri,r2) , 



(167) 



where ri2 = |ri — r2| is the relative distance and R12 = (ri -|- r2)/2 the center-of-mass of the 
nucleonic pair. 

We show in Fig. [32] the proton-proton scalar, q=l, TBDFs, Eq. (|167p . for all the nuclei 
considered as a function of the relative distance of the pair. The dotted lines represent the 
uncorrelated joint probability density of finding the two nucleons at a certain distance, as often 
given in the literature (see e.g. Ref. [66]). These lines have been obtained as a product of 
the uncorrelated one-body densities. This definition of the uncorrelated two-body densities 
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Figure 32: Proton-proton scalar (q=l) TBDFs (|101[) as a function of the distance between two nucleons 
The dotted hncs show the uncorrelated TBDF obtained as a product of two uncorrelated OBDMs. The 
full lines show the uncorrelated TBDF obtained by using in Eq. (llOip Slater determinants. The other 
lines have been obtained with the scalar (/i) and full (/e) correlations. 
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Figure 33: Proton-neutron operator dependent (q>l) TBDFs (|10ip as a function of the distance between 
two nucleons for the ^ospj-, ^.^cleus. 



can be meaningful from the probabilistic point of view, but it is misleading in our framework, 
since it corresponds to using only the first, direct, term of Eq. (jl66p . In fermionic systems 
the uncorrelated two-body density is given by the full expression ()166p which contains also the 
exchange term. These complete uncorrelated TBDF are shown by the full lines. The effects of 
the SRC can be deduced by comparing these lines with the dashed-dotted lines obtained with 
the /i scalar correlations only, and with the dashed lines showing the full FHNC/SOC results. 

In all our results the correlations reduce the values of the TBDF at short internucleon dis- 
tances. The exchange term of the uncorrelated density already contributes to this reduction, but 
the major effect is produced by the SRC, and mainly by the scalar correlations. We found simi- 
lar results for the neutron-neutron TBDF. When the TBDF are composed of different particles, 
the results change only slightly. Beside a strong reduction at small distances, the correlations 
produce enhancements, with respect to the IPM results, around 2 fm, in all the nuclei considered 

To discuss the effects of the correlations on the other operator dependent TBDFs (llOip we 
show in Fig. [33]the TBDFs for the 208p|-, nucleus, for the various operators. We first notice that 
the tensor TBDFs, p = 5,6, are different form zero only when the /e correlation is used. This 
occurs because, in order to get a spin trace different from zero in Eq. (jlOip . at least two-tensor 
operatosr are needed. Also the spin term p = 3 differs from zero in the /g case only. Here there 
are different reasons why both terms of Eq. (jl66p are zero. The first term is zero because of the 
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spin trace, while the second term vanishes because the two nucleons have different isospin. 

The most remarkable feature is the range of the various TBDFs. The scalar TBDFs of 
Fig. [32] and the isospin TBDF of Fig. [33] extend up to relative distances comparable to the 
dimensions of the nucleus, 15 fm. In contrast, all the other density functions have much smaller 
ranges, of the order of 3-4 fm. We found similar results for all the nuclei considered [33\ |53] . 

5.10 Quasi-hole wave functions and spectroscopic factors 

The quasi-hole wave function is defined as: 

I'nljniy-^) >l/2< >l/2 ' 

nljml nljm ^ I ^ 

where ^nij^^i^^ ...,A — 1) and ^{1, are the states of the nuclei formed by A — 1 and A 

nucleons respectively, and is the isospin projector. In analogy to the ansatz ([2]), we assume 
that the state of the nucleus with A — 1 nucleons can be described as: 

Kijmih A-l)= A - A-1) , (169) 

where ^nij^i^, A — 1) is the Slater determinant obtained by removing from ^{1, ...,A) a single 
state characterized by the quantum numbers nljmt. For the system of ^ — 1 nucleons we use 
the same correlation function fixed for the system of A nucleons. In an uncorrelated system the 
quasi-hole wave functions coincide with the hole mean- field wave functions (j88p . 

We are interested in the radial part of the quasi-hole wave function, which we obtain by 
multiplying equation (jl68p by the vector spherical harmonics Yj*™'(J7) and, then, by integrating 
over the angular coordinates O, and summing over m. It is useful to rewrite the radial part of 
the quasi-hole wave function as |32j : 



= ^^j;/df^Y,7(i7)V'*,,„(x) 

= — E 

2j + 1 ^ 



where we have defined: 

w , , _ ^ < KumiA - mtj^imr - rA)P'A\^{A) > 

^nljmdr) - — T^i , 

nljmi nljm 



and 



_ <Kia^^-^A-mi,^^{i,...^A-i)> 

"'J™* < *(1,...,^)I^(1,...,A) > ■ ^ ' 



Following the procedure outlined in Ref. [32], we consider separately the cluster expansions 
of the two terms AA* and . For X^ we have: 



5^:,(r,ri)C*;^,(ri) 



3 terms J\ and X . _bor /C we have: 
K{r) = Ct'^ui^)f^Rii^{r)+Jdr^Rl^i^{n)Pi{cose) 

+p*''^(r,ri)-iVi^)*'"(ri,r) | 

+C^r22(r) / driRii^{n)Pi{cose)X'socir,ri) , (173) 
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and for we obtain: 



6i(r)|2 + j dri0i*(r)0i(ri)2 



9cc (r,ri)C '^^(n 



ri,rj 



+ I dvct^':{r)C%{v) j dri</>i(ri)AA*oc(r,ri) , 



(174) 



where we have indicated with a the set of the nljm quantum numbers. All the FHNC quantities 
have a superscript a since these equations must be built by using: 



A-^i (r,ri) 



■^sbc(r>ri) 



(175) 



Po"(r>ri) = Po(r,ri) 
instead of po(r,ri)- The expressions of M^Q(j{r,ri) and A'_gQ(^(r, ri), are: 

,2fc-i(r.ri) 

A;=l t2=p,n 

+Xf ^ ri) + ^2t2.-i(r. ^i)^ 

+Xf^'^'2t2fc(r,ri)} , (176) 

E {(l-'^M)AA2f42.-i(r,ri) 

A;=l t2=p,n 

+xf^[AA*^i\,2.(r,ri) 

+'^2t2.-i(r. ri)] + xrW*^,5,(r, n)} . (177) 
where the indexes t refer to the isospin and the functions and we have defined: 
X^l'Hr,r,) = i{<*,^'-(r,ri)5i^*-"(r,ri)C*-:;(ri) 
"-/,*-"(r,ri)+<^j,"(r,ri)" 
+5t^*-"(r,ri)C*-^(ri)<',^,(r,rO- 

iVS'"(ri'r)}A^^ , (178) 

A/;VHr,ri) = {/^*^*-"(r,ri)5*:,*-"(r,rOC*-^(rO 

(r,ri)+iV*-"(r,ri)" 

+5^;,*^(r,ri)C*-^(ri)iV*,^;;(r,ri) - iV(^'J^'"(ri, r)} A^^ , (179) 

with q = 2/c2 — 1 + L The other terms are defined in Appendix [Bl 

The knowledge of the quasi- hole functions allows us to calculate the spectroscopic factors: 



'-'nlj 



(180) 
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Figure 34: Protons spectroscopic factors of tlie ^osp^ nucleus as a function of the separation energies. 



The inclusion of the correlations produce spectroscopic factors smaller than one, the mean- 
field value. In general, this effect increases together with the complexity of the correlation. The 
/e results are smaller than those of f^, which are smaller than those obtained with /i. 

We found that in our calculations the effects of the correlations are larger on the internal 
shells [60]. This fact emerges by observing that for a fixed set of Ij quantum numbers, the 
spectroscopic factors increases with n and at the same time, that values of the spectroscopic 
factors become larger when n and Ij values increase. This effect is well represented in Fig. 
[Ml where we show with the black points the proton spectroscopic factors of the ^'^^Pb nucleus 
as a function of the separation energies, defined as the difference between the energy of the 
A-nucleons system and that of the correspondent A-l-nucleons system. We have associated the 
spectroscopic factor of each level to its empirical separation energy. The behavior of the black 
points of the figure indicates that, in our calculations, correlation effects are larger on the more 
bound levels. 

In Fig. \35\ as example of the correlations effects on the quasi- hole wave functions, we show 
the squares of the proton 3si/2 and neutron 3pi/2 quasi-hole wave functions. The global effect is 
a lowering of the wave function in the nuclear interior, and this effect increases with increasing 
complexity of the correlation. 

In Fig. [36]we show with a gray band the difference between the empirical charge distributions 
of 2°^Pb and ^°^Tl The dashed dotted line, labeled as IPM, has been obtained by using 
the fact that the difference between the two charge distributions can be described as a single 
3si/2 proton hole in the core of the lead nucleus. This curve has been obtained by folding the 
IPM line of Fig. [35] with the electric proton form factor in its dipole form. In a slightly more 
elaborated picture, the ground state of the ^os^pj jg composed of the 3si/2 proton hole in the 
206p|j ground state, plus the coupling of the 2^5/2 and 2^3/2 protons levels with the first 2+ 
excited state of ^'^'^Pb [68l[69]. This description of the ^^^Tl, charge distribution, labeled IPM*, 
and shown by the dotted line, is still in a IPM framework. The dashed line has been obtained 
by adding to the dotted line the core polarization effects produced by long-range correlations. 
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Figure 35: Square of the ^"^^Pb proton 3si/2 and neutron 3pi/2 quasi-hole functions. Tlie various lines 
show the results obtained by using different type of correlations. 
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Figure 36: Differences between charge density distributions of 206p]-, ^'^^Tl. See the text for the 
explanation of the various hnes. 



These effects have been calculated by following the Random Phase Approximation approach of 
Refs. \70\ I71j. The full line has been obtained when our SRC effects are also included. 

The various effects presented in Fig. [36] have been obtained in different theoretical frame- 
works, and the final result does not have any pretense of being a well grounded and coherent 
description of the empirical charge differences. The point we want to make by showing this 
figure is that the effects of the SRC are of the same order of magnitude as those commonly 
considered in nuclear structure calculations. 
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6 Perspectives 

In this section we give a short review of the possible developments of the theory. We present 
only those topics which have been already formulated. Some of these subjects have been already 
well studied in nuclear matter, and the work to be done consists in adapting the formulation to 
the case of finite nuclei. Other topics are still at the level of a very abstract general formulation, 
valid for any kind of many-body systems. We first present some possible applications of the 
theory, which do not require changes in its basic hypotheses. Then, we discuss some extensions 
of the theory. 

The main goal of the RFHNC/SOC equations is the evaluation of the hamiltonian expecta- 
tion value, in order to apply the variational principle. Once the parameters of the correlation 
function and of the mean- field potential which minimize the energy functional, Eq. ([T]), have 
been found, it is relatively straightforward to apply the RFHNC/SOC equations to evaluate 
expectation values of other operators. This is the strategy used in Sect. [U to calculate all the 
quantities other than the energy. 

So far our approach is aimed at the description of the nuclear ground state, therefore we can 
only obtain expectation values of ground state observables. On the other hand, a clever use of 
completeness relations allows us to get information on excited states by calculating expectation 
values of operators between ground states. This is, for example, the case of the sum rules. The 
enhancement factor of the electric-dipole sum rule has been calculated in nuclear matter |72j . 
and the same approach can be applied to finite nuclei. 

Dynamical response functions and hole spectral functions have been calculated in nuclear 
matter by using the FHNC/SOC formalism [231 CSl EH ES]. The responses of the system for a 
momentum transfer q and an excitation energy u) have been evaluated by using the expression 



where and Eq are the ground state wave function and energy respectively, H the nuclear 
hamiltonian and p the external operator exciting the system. Analogous expressions have been 
used to calculate the spectral functions. Also in this case the completeness of the excited 
states has been used, and the response is expressed as the ground state expectation value of an 
operator. Formally, the cluster expansion and the RFHNC/SOC resummation techniques can 
be applied without any major problems to evaluate these expectation values. However, from 
the pragmatical point of view, we have to consider that the expression of the global operator 
is extremely involved, as is shown by Eq. (118ip . This global operator combines the operator 
describing the effect of the external probe, usually a relatively simple one-body operator, with 
the hamiltonian (11231) composed by one- two- and three-body terms. The effort will be rewarded 
by the results, since the evaluation of the response functions gives direct access to the calculation 
of the cross sections, and of other observables. 

The RFHNC/SOC theory can be applied to describe hypernuclei by considering the hyperon 
as an impurity in the nucleonic fluid. The FHNC equations for an impurity in homogeneous 
matter have been derived in Refs. [771 178j to describe the presence of atomic ^He in liquid ^He. 
A first application of this theory to the case of single A hypernuclei has been done in Ref. [79]. In 
these calculations we used simple interactions and correlations. The nucleon-nucleon interaction 
contained only the first four central channels, and the A-nucleon interaction the scalar and spin 




(181) 
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channels only. All the correlations were purely scalar functions. 
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Figure 37: Binding energy of the A hyperon for the Is, Ip and Id states as a function of A^^^^. The 
white circles are the energies calculated as indicated in Ref. [79j by using a Woods-Saxon mean-field well 
for the hyperon. The experimental energies are from Ref. [5D] (triangles) and ^81j (squares). The full 
lines connecting the theoretical values have been drawn to guide the eye. 

Despite the simple ingredients used in the calculations, the agreement with the experimental 
A binding energies of Refs. [80\ [HT] is rather good, as is seen in Fig. [371 The extension of the 
formalism of Ref. [7S] to the case of operator dependent correlations, which allows us to deal 
with realistic hamiltonians, is technically rather involved. On the other hand, the results of Fig. 
[571 are very encouraging and the potentialities of our technology to make predictions are wide, 
and rather unique. 

In the introduction, we have mentioned that our approach is the lowest order approximation 
of the CBF theory formulated at the beginning of the '60s [82 1183118^ . The starting point of the 
CBF theory is the construction of a basis of normalized, but in general non-orthogonal, state 
vectors of the type: 



where F is the many-body correlation function acting on a complete orthonormal set of model 
states >• The CBF theory constructs a coherent perturbation theory on the correlated 
basis formed by the states given by Eq. (|182p . [13j. In the limit F — > 1 the application of 
the variational principle would provide the Hartree-Fock equations, and the excited states (|182p 
would be constructed as particle- hole excitations of the single particle basis. The so-called 
residual interaction would mix these single particle excitations in a perturbative expansion. It 
appears evident how our calculations, based on the application of the variational principle ([T]), 
can be considered, in the CBF framework, as the first step of a perturbation expansion. This 
first step is necessary to fix the correlation function F and the wave functions basis. 

An improvement of our work consists in the inclusion of higher-order perturbative terms 
in the CBF expansion. This has already been done for nuclear matter, where perturbative 
corrections to the binding energy [85], momentum distributions [22], responses |23 1 174 1 175]. and 



^^(1,...,^) > 




(182) 



< «>„|FtF|«>„ >V2 
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spectral functions [73l[86], have been evaluated. With the help of the perturbative expansion, 
the nucleon-self energy in nuclear matter has also been calculated and the optical potential has 
been evaluated [871 ES] . 

The use of the CBF theory allows the description of nuclear excited states. In this respect, 
they can be treated within the correlated Random Phase Approximation theory, whose basic 
equations have been obtained by using a time-dependent Hartee-Fock approach \88 \ \89 \ [90l I13j . 

Another line of development of our theory consists in modifying the Jastrow anstaz ([3|) . For 
example, correlation functions composed of scalar Jastrow function and linear state-dependent 
correlations have been proposed in Ref. [JT]. These correlations have been applied to the 
description of light nuclei using Variational Monte Carlo techniques with promising results [911 
I92j . With linear state-dependent correlations, the structure of our RFHNC equations becomes 
simpler. However, this approach has serious drawbacks, since it requires the use of up to six-body 
distribution functions. 

The developments described so far are thought to be applied to doubly closed shell nuclei. 
The description of other nuclei is related to a change of the single particle basis, which should 
consider open shell and deformations. Also in this case, there are no problems in principal in 
applying our theory, but the technicalities are rather involved. 
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7 Closing remarks 



The theory we have presented aims to describe the properties of medium-heavy nuclei by solving 
the many-body Schrodinger equation with microscopic nucleon-nucleon interactions. This is an 
ab initio calculation, since there are no free parameters to be adjusted. However, the technique 
of solving the many-body Schrodinger equation requires some approximations. 

The most basic approximation is related to the use of the variational principle. The search 
for the minimum of the energy functional, Eq. ([T|), is done by spanning the space of many- 
body wave functions which can be expressed as a product of a correlation function times an 
independent particle model wave function, see Eq. ([2]). The many-body correlation function 
is expressed as a product of two-body correlation functions which hinder two nucleons from 
approaching each other for distances smaller than the range of the strongly repulsive core of 
the nucleon-nucleon interaction. The use of the ansatz ([2]) allows an explicit treatment of a 
great number of correlations. However, these are not all the possible correlations, and they are 
specifically related to the strongly repulsive core of the interaction. Other types of correlations, 
such as those produced by collective motions of the nucleons, are not so well described. Since this 
is the starting hypothesis of the theory, we cannot test its validity within the theory itself. The 
only possibility of making this test, is a comparison of our results with those obtained by other 
microscopic theories which are approximation free, such as Quantum Monte Carlo calculations. 
The validity of the variational ansatz has been tested in the literature by comparing Variational 
Monte Carlo and Quantum Monte Carlo results for light nuclei. This comparison, done for nuclei 
with A<8, indicates a minimum difference of 0.1 MeV per nucleon in ^He and a maximum one 
of about 1.0 MeV per nucleon in ^Li |i93j. Hence, this difference seems to increase with the 
number of nucleons. 

The evaluation of the variational energy functional ([T]) in Ref. [93] is done with a Monte Carlo 
integration technique, i. e. without approximation. In our approach we calculate the energy 
functional by doing a cluster expansion. After a topological analysis of the various diagrams, we 
sum two categories of diagrams, the nodal and the composite ones, in closed form, by solving 
the RFHNC integral equations. The procedure we have used does not consider a certain type 
of diagrams, the elementary ones. This is another approximation of our computational scheme. 
The role of the elementary diagrams has been studied in quantum liquids and in nuclear matter, 
and they have been found to be more important in the former systems than in the latter ones. 
This is because, in many-body jargon, liquid helium is denser than nuclear matter. This means 
that the number of particles in a volume characterized by the range of the repulsive core of the 
interaction, is larger in liquid helium than in nuclear matter. 

We have calculated in Ref. [27] the contribution of the elementary diagram of Fig. [HI the 
simplest one. The calculation has been done for a simple interaction, the Bl force of Brink 
and Boeker with scalar correlations and in model nuclei with isospin degeneracy and single 
particle wave functions treated in the Is coupling scheme. We obtained for ^^O and ^'^Ca a 
repulsive contribution of about 1.1 MeV per nucleon. This is about 2% of the contribution of 
the potential energy produced by the interaction. However, the total energy is obtained as a 
difference between kinetic and potential energies, and, in this case, 1 MeV is not negligible. The 
contribution of the elementary diagrams should be further investigated. 

The third approximation of our computational scheme is the SOC. It is clear that the ad 
hoc SOC approximation breaks the formal completeness of the FHNC theory which holds when 
only scalar correlations are used. We have difficulties in making estimates of the contribution to 



79 



the energy of the diagrams excluded by the SOC approximation. We can only say that the sum 
rules exhaustions are slightly worsened when the SOC approximation is used. We are talking in 
any case of few parts on a thousand. On the other hand, we cannot exclude that the diagrams 
not considered are irrelevant for the sum rules, but not totally negligible for the calculation of 
the energy. 

The three approximations just discussed are somehow intrinsic to the computational scheme. 
We now talk about the simplifications we have done in the specific applications of the RFHNC/SOC 
formalism. In our calculation we have used a nucleon-nucleon interaction limited to the first 
eight channels. Modern interactions, with isospin symmetry breaking terms, have up to eighteen 
channels. Our correlation functions are also limited to the first six channels. We have estimated 
the contribution of these missing terms of the interaction and of the correlation by doing a 
bold extrapolation of nuclear matter results. The results of this estimate are shown in Tab. 
[HI The contribution of the neglected channels is small when compared with the contribution of 
the interaction term, V2_i,g^y in Tab. [5j However, since the energy is obtained by a subtracting 
potential energy to the kinetic one, their contribution is not negligible on the final result. 





i«o 4uc^ 48ca ^u«pb 


< kp > 
E 

AE 

E+AE 


1.09 1.09 1.19 1.19 1.21 
-0.91 -4.48 -4.58 -4.14 -3.43 
-0.93 -0.93 -1.37 -1.37 -1.48 
-1.84 -5.41 -5.95 -5.51 -4.91 



Table 8: Average Fermi momenta < kp >, in fm~^, used to estimate the corrections to the binding 
energies produced by the interaction channels beyond the spin-orbit ones, from nuclear matter calcu- 
lations. All the energy values are expressed in MeV. The E rows gives the results of Tab. [5] for the 
AV8'+UIX interaction. The AE rows values have been obtained by doing a local density approximation 
interpolation of the nuclear matter results given in Tab. Ill of Ref. [31) . 

Our minimization procedure has been done on two parameters only, the healing distance for 
all the central correlation channels, and that of the two tensor correlations. In principle, our 
model could handle six variational parameters for the correlation, plus other thirteen describing 
the mean field Woods-Saxon potential generating protons and neutrons single particle states. 
The choice of using only two variational parameters, dictated mainly by computational reasons, 
could seem rather limiting. In reality, we are confident that our minima are rather close to the 
minima that a full minimization with all the variational parameters could obtain. Concerning the 
correlations, there is a nuclear matter experience [21] indicating that the four central correlations 
heal at the same distance, noticeably different from the tensor healing distance. We have done 
some calculations in '^'^Ca by changing the healing distances of all the correlations, and we did 
not find significant improvements of the energy value. The case of the single particle basis has 
been discussed with some detail in Sect. 15.21 Also in this case, we have tested that, when the 
minimization on the correlation function has been done, even large changes of the mean-field 
potential do not produce sensitive changes in the energy minimum. 

At present, the most evident problem of our calculations is related to the behavior of the 
three-body force. In few body and light nuclei the contribution of this part of the hamiltonian 
is attractive, while in our calculations it is repulsive. This happens also in variational nuclear 



80 



matter calculations |24[|25j . There are two possible solution of the puzzle. The one which should 
be explored first, is that our present treatment of the three-body force is not accurate enough. 
For example, we could find that the set of diagrams we are now considering, see Sect. 14. 3^ is 
not sufficient, or that the inclusion of a three- nucleon correlation function is necessary. If after 
improving the description of the three-body force, its contribution remains repulsive, we have 
to deduce that the three-body interactions we use, are tailored to provide good descriptions of 
few body systems and light nuclei, but they are not adequate for medium, heavy and infinite 
nuclear systems. 

The three-body force puzzle, is an example of the potentiality of our approach, and, in 
general, of all the microscopic calculations, in nuclear structure. These calculations allows us to 
investigate phenomena that mean-field based effective theories cannot study. We have pointed 
out a few examples, by showing the effects of the SRC on momentum distributions, natural 
orbits, quasi-hole wave functions and spectroscopic factors. But the relevance of microscopic 
calculations goes beyond that, since these calculations have reached such an accuracy to put 
constraints on the nuclear hamiltonian itself. Microscopic nuclear structure calculations will 
become more and more important in the near future, and, among them, the RFHNC/SOC 
computational scheme will play a relevant role. 



81 



Final note 



We had just started to write this article, when Adelchi Fabrocini passed away. The reader who 
had the privilege of knowing him is certainly aware of the absence of his touch in the writing 

of this article. Adclchi's contribution to the realization of the work presented here, has been 
enormous and fundamental. We miss his talent, his leadership, his experience and, not least, 
his subtle sense of humor. 
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A Matrices of the spin and isospin traces 

In this appendix we give the numerical values of the various matrices used to calculate the spin 
and isospin traces, when state dependent correlations are used. 
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5 
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1 


1 









Table 9: The values of defined in Eq. (j58|) . The values of p indicate the operator channel. 
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18 



Table 10: The values of BP defined in Eq. 
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Table 11: The values of the matrix K^"^^ defined in Eq. ([63]) . The values oi p,q and r indicate the 
operator channel. 
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Table 12: The values of the matrix L^'''^ defined in Eq. (f67|) . The values oi p,q and r indicate the 
operator channel. 




Table 13: The values of the matrix Epq defined in Eq. (|69p . 
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B The RFHNC/SOC equations for nuclear finite systems 



We present here the set of RFHNC equations for finite nuclear systems with different number 
of protons and neutrons, and in a jj coupling scheme. The upper index t distinguishes the 
contributions of protons and neutrons. In this appendix, in addition to the terms necessary 
to calculate the one- and two-body density functions, which are used for the evaluation of the 
energy of the system, we present also those terms which are required in the evaluation of the 
OBDM, necessary ingredients for the calculation of the momentum distribution. 

We show the expressions of the various parts necessary to build the one- and two-body density 
functions used in the calculation of the energy of the nucleus, and also those parts required to 
obtain the OBDM necessary for the momentum distribution evaluation. 

The various terms can be written in general as (rj,rj). An important simplification in 
the writing of the equations is the property: 



Yr 



titj 



y-tjti 
nm 



(183) 



We present first the set of RFHNC equations for the scalar part of the correlation /i, and, in 
a second step, the equations involving operator dependent correlations, i.e. the RFHNC/SOC 
equations. 

We start our presentation by considering the dynamical diagrams, i.e. those where the 
external points are reached only by dynamical correlations. We used two types of dynamical 
correlations, hd = fi — 1, and = /i — 1. This last correlation appears in the calculation of 
the OBDM, where it connects only the external points. The dynamical TBDFs can be written 



as: 



5td'(ri>r2) = /i(ri2)exp[<^,*^(ri,r2) + Std'(ri,r2) 
gll,'j (ri , ra) = exp (n , rs) + E'^'J (n , rs) 



(184) 
(185) 
(186) 



These distribution functions are related to the non-nodal diagrams X*!*^ by the equation: 

9iki'{rl,r2) = l + N^^i^{r^,r2)+X'^i^{r^,r2) , (187) 

with m,n = d,w. In the presence of purely scalar correlation functions, the dynamical functions 
do not depend on the isospin of the external particles. However, in the above equations, we wrote 
the explicit isospin dependence in view of the treatment with the state dependent correlations. 

The next step is to consider the case when the dynamical correlations link only one external 
point, let's say point 1, while the other point is reached by the statistical correlations, forming 
an exchange loop. In this case, the functions to be considered are: 



trt2 
yme 



(ri,r2) 



(ri,r2) 



N, 



tlt2 I 



ri,r2) + i?^^*=^(ri,r2 



= iV^*2(ri,r2) + X^*Hri,r2) 



(188) 



with m = d,w. 
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In the case when both the two external points are reached by exchange loops, the equations 



are: 



5*^e!L(ri,r2) = glf{r,,r2)\Nll'Hr,,r2)+EXUri,r2) 



(189) 



+ 



tlt2 

Hee.exc 



(ri,r2) 



ri,r2)+^:;,^^(ri,r2 
(iV*f (ri,r2)+4f (ri,r2)) 
ff^f (ri,r2)[^*^*^,e(ri,r2) 
(nII (n , r2) + E% (n , r2) - p*^ (ri , r2 
(n'^{v^,V2)+E%{v^,V2) - P^^(ri,r2)) 



(190) 



(191) 



,exc] 



e: 



1*1*2 

ee,excj 



(ri,r2) 



A^cci (r 1 , r2 ) + (r 1 , r2 ) - p^) (ri , ; 
(iV,%. (ri , r2) + E% (n , ra) - pg- (ri , ra)) 



(192) 



In the above expression the subindex dir denotes the case when the two external particles are 
linked to different statistical loops, and exc{j) when they are connected to the same loop. In 

this last case, a spin-isospin exchange operator is present, whose trace in the isospin space is 
25t^t2- III the exc(j) case, we further distiguish the parallel spin case, exc, from the antiparallel 
one, excj. 

Finally, in the construction of the exchange parts, we have to introduce the contribution of 
diagrams with open statistical loops. These open loops appear in both the calculation of the 
energy and of the OBDM. We have to distinguish, also in this situation, the case where the 
external particles have parallel or antiparallel spins. In either case, the exchange loops may be 
combined with the two kinds of dynamical correlations, therefore we define: 



<?*i(^)(ri,r2) 



9 



WcWcU) 



(ri,r2) 



(ri,r2) 



5^^inri,r2) 

iV*^(,.) (ri , r2) + E'^^.^ {r^ , r2) - pl\.^ (n , r2] 
5^^i^(ri,r2) 

9ww {^1,^2) 

(n, r2) + E'l^^^.^ir,,r2) - p\^){vuV2 



(193) 
(194) 
(195) 



where with (j) we have indicated the possible presence of the label j. For all these six partial 
distribution functions we can write: 



9%^^j){ri,r2) = iV^„(^.)(ri, r2) + X^„(,)(ri, r2) - P*^(^.)(ri, r2) 



(196) 
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The RFHNC expressions of all the nodal diagrams which do not have open statistical loops, 
can be written in a compact form as: 



E E 

ts=p,n m'n' 



(197) 



where m,n = d, w, e. In the previous equation, due to the diagrammatic rules, the sums are 
limited to the values (m'n') = dd, de, ed. Furthermore we have defined: 



(i) 

' mn\ J 



C*^ (i) for (mn) = dd 
(i) otherwise 



(198) 



The expressions of the nodal diagrams with open statistical loops in the external particles, 
such as the diagrams of Fig. [T^ follow the classification presented in Eq. (jl04p : 



iVitHri,r2) 



^^2?(ri,r2) 



A^S?(ri,r2) 



+ 
+ 

+ 



X^,(ri,r3)C*i(r3) 



9ck{r 3,^2] 



^X(ri,r3)C*Hr3) 5*^,(r3,r2)) , 

- (p*^ (ri, r3)C*^ (rs) |iVi;?^*^ (r3, r2) + X*Ur3, r2) 
p*^(ri,r3)(C*^(r3) - l)|iVi^)*nr3, ra) - pl'{r,,r2: 
4(ri, r3)C*^ (r3)|A^i:f (rg, rs) + X*;^.(r3, rs)) 
p*;,.(ri,r3)(C*Hr3) - l)|ivi:J*^ (rg, r2) - 4(r3,r2 
X^,^.(ri,r3)C*^(r3)|5*Ur3,r2) 
X^,(ri,r3)C*^(r3)|g*^^.(r3,r2)) , 

- (4(ri, r3)C*^ (r3)|iVii5-^*^ (r3, rs) + X*J,(r3, r2 
-*Uri,r3)(C*^(r3) - 1) iVi^^^^rs, r2) - pI'{ 



(199) 



(200) 
(201) 



Po 



r3,r2. 



^(r3,r2)+X*J,^.(r3,r2) 

(p)t 



^PoKri,r3)C*^(r3) 

;p*^(ri,r3)(C*^(r3) - 1) iV.^^f (r3, rs) - /.*^(r3,r2; 
where m,n = c,Wc- In the above equation, we have used the definitions: 



(202) 



^2-.(ri,r2) = NZ%ir,,r2) + C^^^, (ri , r2 



mn{j) 



mn{j) 



mn(j) ^ 



and 



iV^''^*A(ri,r2) = NP'''}.Jri,r2) + N""'] .Jri,r2) 



It is worth to point out that Eq. (j200p indicates that N^^^\~^{ri,r2) does not depend on m. 



The results of Sect. 13.21 indicate that the vertex corrections are given by: 
C7^(ri) = exp[[/^(ri)] , 

C*Hri) = C*^(ri)[C/iHri)+p*^(ri)]=p(ri) , 



(203) 
(204) 
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with m = d,w. In order to simplify the writing of the RFHNC expressions for Ud,w,e, we have 
defined the quantity: 



(205) 



and we obtain: 



t/^(ri) 



tlt2, 



t2=p,n 

-Sl'i{r,,v,){g'^'i{r,,r,)-l] 

+C*^(r2)[x^*^(ri,r2)-i?^^*Hri,r2) 
-5^*Hri,r2)(<7^^^(ri,r2)-l) 
-5^^^(ri,r2)5^^*^(ri,r2)]}+i^^(ri) , 

t2=p,n 

-Slf{v,,v,){g\f{v,,r,) - 1) - S'i^-{v,,v^)glf{v,,v,] 

+C*^(r2)[x*^*^(ri,r2)-i?*^*^(ri,r2) 
-Sll'^{v,,r,){g\f{r,,V2) - 1) - Si^^\v,,V2)g\t {v.^v^] 
(ri,r2)5*;,*^(ri,r2) - S',f{v^,V2)gi'^-{r^,r2] 



(206) 



-p*Hri,r2)(iVi^')*Hri,r2)-p*^(ri,r2)) 
-p*^(ri,r2)(ivig*^(ri,r2) - p%{r,,v,))\ } + i^^^ri) , 



(207) 



with m = w,d. The expressions (jl 8411207]) form the close set of RFHNC equations, valid for the 
scalar part of the correlation. In terms of these quantities and defining other useful ones, we 
can express the TBDF and the OBDM as: 



PVir(ri>r2) 

P2'eL(ri,r2) 
P2'eLj(ri,r2) 
/>*nri,r2) 



= P^Kri, r2) + 26t,t2 {ptlU^i^r^) + p'iX^{r,,r2)) 

= C*Hri)(C*^(r2)5rf:,*^(ri,r2) +C*^(r2)5^^,*^(ri,r2)) 

+ C*^(ri)(C*nr2)5*^;^(ri,r2) + C*^(r2)5*i5,,(ri, r2)) 

= C'iir,)Cfir,)gil%,,{r,,r,) , 

= C';{r,)Cl^{r,)gll%^^ir,,r,) , 

= -2C4Hri)C^Kr2)5^^Kri,r2) 



(208) 

(209) 
(210) 
(211) 
(212) 



We discuss now the case of state-dependent correlations, in the SOC approximation. We have 
mentioned in Sect. 13.31 the need of separating the spin and isospin dependence of the TBDF in 
order to describe nuclei not saturated in isospin. In these systems, the contribution of the linear 
isospin operators is not zero, therefore we distiguish the proton and neutron dependence of the 
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various RFHNC/SOC terms. This affects the chain equations in the calculation of the nodal 
A^*J^*|,(ri, functions. In order to generate these chains, we have to consider the following 
folding products: those of a function Xxy^p{ri,r3) with another function Xa;j/,g(r3, r2), those of 
^a;j/,p(ri,r3) with 7V^.y,q(ri,r3), and those of Xa;j^,p(ri, rs), or A^a;y,q(ri, ra), with /9o(r3, r2). These 
combinations are present also in isospin saturated systems. In addition, we should also consider 
that the action of a single isospin operator on a single external point of the nodal function, 
produces a non zero contribution. 

We give in the following the expressions of the vertex corrections of the nodal terms in SOC 
approximation, and also those of the various terms of the TBDF. In these expressions we use 
the index k which can assume the values 1,2 and 3 and the index / which can be or 1. The p, 
q and r indexes indicate the operator channels and can assume values from 1 up to 6. 

The expressions of the nodal functions without open statistical loops can be written as: 



:2fci-ilri,r2j 



3 

E 

k2k3=l ti=p,n 



+ (25*,i3-l)A^^Sk-i,2.„2fc3-i(ri'r2) 

+ (25,,i3-lX*f,^_,_2..-l,2fc3(ri'r2) 



^^^;!2fc,(ri,r2) 
where we have defined: 



E E 

^2,^3 = 1 t3=p,n 



(213) 
(214) 



iV^*^*|,(ri,r2) = (x\X{^r.^zX\t'''C%{v^)^^^^ 

(x^*3^(ri,r3)Cil2'^'=^C*^,,(r3)|x*f,(r3,r2) + iVK^,(r3,r2) 
(x^*^^(ri,r3)Ct|2'^'=^Cj^,(r3)|x*0(r3,r2) + N^-iv^^v^) 



+ 
+ 

(215) 

In the above equations we have defined p = 2ki — 1 + ^i, q = 2k2 — 1 + ^2; and r = 2A;3 — 1 + /3. 
The sub-indexes m and n indicate the type of link with the two external points, specifically 
m,n = d, w, e. 

Since in the calculations of ww diagrams, we must include the isospin trace, for this case, 
we substitute Eq. (j214p with: 



<i2.,(ri,r2)= E E ^2''Kt2ku2k2,2kM^^2) ■ 
k2,k3=l t3=P,n 

The expressions of the nodal diagrams with open statistical loops are: 

N'r^n,pi'^ur2) = N^l%{r,,r2)+N^l%{rur2) p=l,...,6 



(216) 



(217) 
(218) 



with m,n = c, Wc- We have used a symbology analogous to that of Eqs. (j200l - [202]) . to indicate 
the parallel spin case, Eq. (j217p . and the antiparallel one, Eq. (j218p . In this last case, we have 
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considered only the contribution of the first two channels of the interaction, p = 1,2, in order to 
simplify our calculations. This is a good approximation since the relevance of the antiparallel 
loops is small. The explicit expressions of the nodal functions present in the Eqs. (j217p and 
(f2lk are: 



iv(^)*\,, i(ri,r2) 



[^mn(j),2fci-l, 2^2-1,2^3-1 



^2,^3=1 t'i=P,n 

+{2Kt, - i)(<T(;;k-i,2.„2.3-i(ri.r2) + 
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iV^^^*\,, (ri,r2) 



where y = x, p. The other terms are: 
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/ A ^2 

+ 5.,iUnri,r3)— Cf3t'=HC*^o.(r3) - 1) r'Uv^.v^) 



+ 

+ 



'5fcil5fc2l'^fc3l 



Poj(ri, r3)Cjgr(r3) X*^^- ^(r3, 



pfe(ri,r3)Cj^,(r3)|A^*:^,Vr(r3,r2)) 



p*Uri,r3)(C*y(r3)-l) 



^cni,r(r3,r2; 



5,,i(p*^(ri,r3)(Cj^,(r3) - 1) pg.(r3,r2) 
X^,Jri,r3)C^,,(r3)i5*?.,(r3,r2) 



(1 - 5o)(x^,(ri,r3)C7j,,(r3)i|5X-,(r3,r2) 
^l.s(ri,r3)C*^,,(r3)^ 



5cn (r3,r2; 



+ (1 - 5o) XX.(ri,r3)Cj,,(r3)- 5c?.,*(r3, r2; 



(219) 
(220) 



(221) 



(222) 



(223) 
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mnjjpst V 



(ri, r3)C*3^,(r3)^|x*3^.,(r3, r^) + ^^^^(ra, r^)) 

- (poUri,r3)(C*^,,(r3) - l)^|iV^g-5(r3, r^)) 

- (p*^(ri,r3)Cj,,(r3)^|x*,«,,,(r3,r2) +ivig^(r3,r2)) 
+ 5*,i(poMri,r3)(Cj,,(r3) - l)^|p*^(r3, r^)) . 



(224) 



In the above equations the symbols s and t can assume the values 1 and 2. 

In the calculations of WcWc diagrams, we must include the isospin trace, therefore we use: 



WcWc,2ki 



fc2,fc3=l t3=P,n 



(225) 



The expressions of the TBDFs are rather similar to those of the symmetric nuclear matter 
case. We define the quantities: 



2/p(ri2) jyM2/ ) 



^S(ri,r2) 

/^Ltp(ri'r2) 

with p = 2A; - 1 + / > 1. We obtain: 
9^n,p(ri,r2) 



(226) 

(227) 
(228) 



5j^e!p(ri,r2) 

5ee!p(ri,r2) 



5ee3ir,p(*'l'''2) 



5r*.m(ri,r2)/?4m,p(ri,r2) 

5^^*^ (ri, r2)/i^*4(ri, r2) + g'^'^ (n, r2)iV^*^p(ri, r2) 

(ri,r2) +iV^*fp(ri,r2) , 
5ee!L,p(ri, r2) + 5ee!Lc,p(ri, r2) + 5ee!Lci,p(ri> r2) 
X*i;^(ri,r2) + iVi^*^(ri,r2) , 
5ee5ir(ri,r2)/?-dip(ri,r2) 



(229) 
(230) 



+ 5,r(ri,r2)iV*^,:|(ri,r2) 
+ (ri,r2)iV*^5(ri,r2) 



+ 5:r(ri,r2)A^^;;^(ri,r2), 



,exc,p 



with m,n = d, w. 



ed 1^1' 
A''<7ee!lxc(ri,r2) , 

A'5SLc,(ri,r2), 



(231) 
(232) 
(233) 
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For the other diagrams we can write: 



5*^(,),p(ri,r2) 



5*^a)(ri,r2)/i^;(ri,r2) 



+ 5^^*^(ri,r2)iV*^(^.),,(ri,r2), 
= 5t'^c(j)('^i'r2)/iL'Sp(ri,r2) 

= 5!,L.,,(j)(ri>r2)/iL'4\p(ri,r2) 



(234) 
(235) 
(236) 



All the TBDFs can be expressed in terms of nodal and non nodal diagrams, and for all of 
them, we can write: 

0*1 (ri,r2)=X*i (ri,r2)+iV*i (ri,r2) . (237) 
The expressions of the vertex corrections for the operator dependent part of the correlations 



are: 



CXv,{^i) = C^(ri) 



1 + Vll^m,50c(ri) 

+ ViiC^dHri) [ul]soci^i) + ^j,50c(ri 
with m = w,d and where we have used: 

A;i = l t2=p,n K. 

+ Xf ^ K2\,-i,2fe,(ri) + C^SV,2.,-i(ri) 
+ Xr^f^r2\,2fc,(ri)| , 

feiA;2=l t2=p,n y 



(238) 
(239) 



(240) 



+ Xr^t^5,2\„2fc.(ri) , 



(241) 



and we have considered m = d,w, e, the relation p = 2ki — 1 + li and q = 2^2 — 1 + l2- The 
expressions of the U coefficients are: 
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+ (ri,r2)Cj^,(r2)(l-(^,,i)iV*^;,^(ri,r2)| , 



(242) 



Utlt2 /j. \ 
^ w,pq\ >-) 



+ g'^l- (ri, r2)Cj^^(r2)] <$,(ri, r^) 

+ 5Lt(ri,r2)C^,(r2)A^4^,(ri,r2)| , 

KtM) = |dr2/i*i*^(ri,r2)|[5*:,*^(ri,r2)C*|(r2) 

+ 5S.r(ri,r2)C^,(r2)]/.*^*^(ri,r2) 

+ (ri,r2)C7*^,(r2)(l - <5,,i)^et.'(ri, r2) + 5^i,*^(ri, r2)C^^(r2) 

X (1 - S,,,) te(ri,r2) + iV*^5(ri,r2)iV*^*^(ri,r2) 



(243) 



ryt 1*2 



+^ei:^(ri,r2)iVr'(ri,r2 

+ {l-dp,,)A''- 1 dr2|/i*i*^(ri,r2)25*i,*^,,(ri,r2)C^,(r2) 

- iV*^,^(ri,r2)5*?,(ri,r2)C*^,(r2) 

+ iV(4*Hr2,ri)[A^(^')*^(r2,ri)-p*^(ri,r2)]| , 

= (1 - <5,,i)A'=^ I dr2|/i*i*^(ri,r2)25:2.c,(ri,r2)C7^,(r2) 

- iV*^,^(ri,r2)5*^,,.(ri,r2)Cjp^(r2) 

+ iV^*^ (r2, ri) [ivg*^ (r2, n) - p%{r,,r2 



(244) 



(245) 



In the above equations we have used the functions 



/.r^(ri,r2) = |M + (i_<5,_,)iv*i|^(r,,r2), 



/i(ri2) ' ' 

In addition we remark that, since there is no contribution in the case p = g = 1, we have 

C'r^,pqir2) = C^,22(r2) . 
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C The isospin matrix elements 



In Sect. 13.3.21 we have shown how to calculate the expectation value of the diagram A of Fig. 
[T5] by using the properties of the Pauli matrices. We found a recursive relation which allowed 
us to express a set of n pairs of isospin operators as a scalar term plus a single isospin operator 
pair, Eq. (jll7p . This recursive relation was used to evaluate the expectation value of a set of n 
isospin operator pairs between two external points, Eq. (112ip . 

In this Appendix we calculate the isospin expectation value for the situations represented by 
the B and C diagrams of Fig. [151 We call these diagrams vertex correction and nodal diagram, 
respectively. 

The starting point is a general expression for the expectation value of a product of isospin 
operators between the external points 1 and 2, and a generic internal point, we label it as 3. 
Since this expression is symmetric in the external points, we write it in a general manner and 
we understand that the i and j points can be either 1 or 2. We define the expectation value: 

(n • T2)'-Hr, • T3)'^(ti • T2t^xtA^)XtA'^)XtA3) 
b'l + b25t^t2 + b'^kit'i + b'^St^ts + b'^^titi^tit-A^tz , (246) 

where the h'^^ 5 coefficients are real numbers depending on Zi,...,5 which indicate the number of 
isospin pairs appearing in the expression. For the properties of the Kronecker's 5 symbol we 
have that: ^ 

Kt2Kt35t2ts = 2 (^t^t2 + Kt3 + ^t2t3 " 1) , (247) 

therefore we obtain: 

KtiXh ii^j)=bi + b2Kt2 + bz6t,t3 + b,5t,t3 , (248) 

where we have defined hi = h'^ — 65/2 and = h'^, + h'^/2 for k = 2,3,4. By exchanging the 
coordinates 1 and 2 in Eq. (1246P we find that: 

-^IlSat/s (3 - ^> 3 - ,) = Tl-^^^ i^, j) . (249) 

Eq. (j246p represents the most general expression for all the cases we want to treat. All 
the possible combinations of isospin operator pairs of our calculations, can be reconducted to 
this expression, by appropriately redefining the values of the I powers, and those of the i and 
j coordinates. The structure of the sequence of operator pairs can be interpreted as follows. 
The isospin pairs with exponent li and I5 represent the operators coming from the dynamical 
correlations, and, eventually those coming from a statistical correlation line applied to the points 
1 and 2. The pair with I3 represents the operator of the interaction. The isospin pairs with 
power I2 and /4 are associated to the operators of the vertex correction in the case j = i, and 
to the nodal diagram in the case j = 3 — i. Since we work in SOC approximation, there is only 
a single pair of isospin operators acting on these particles, therefore = 0, 1. On the other 
hand, there are no limitations on the values that the other I indexes can assume. We discuss 
below the three possible cases, compatible with the SOC approximation. 

When ^2 = ^4 = we have a structure analogous to that treated in Sect. 13.3.21 and therefore: 

'?;:J)SoI(^'J-)=x[J^.3+z, , (250) 
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with x'l%^+i, given by Eq. 

For the case I2 7^ h, obviously one of the / indexes is 1 and the other one is 0. The 
expression of Eq. (j246p can be written, with an appropriated redefinition of the power indexes, 
to an expression of the kind: 

{Ti-T2yHri-T3){Ti-T2y'XtA^)XtMXtM ■ (251) 

These coefficients are related to the T coefficients by the relations: 

Koi%S''3) = Tllltlij)^ (252) 

By using the recursive relation (jllSp . we find for r^*'^^*^*^ (i) the following expression: 

Tl^\f-'{i) = ai,ai,Tl}^^'^{i) + [ai,{l - m,) + (1 - a^Ja,,] T^^^^^ (i) + 

{l-ai,){l-ai,)Ti\'^''^{i) . (254) 



h h 


rptlt2ts /-I ^ 

^hh 

bi 62 ^3 h 




1 

1 1 


-10 2 
-10 2 
-10-2 4 



Table 14: The values of the b coefficients of Eq. (|248p used to calculate the basic T/^/^*^ (1) terms of Eq. 



The values of the 6 coefficients for each term of Eq. (j254p are given in Tab. [U] and have 
been calculated by using the expression: 

■ TjXu ii)Xt, (j) = 2xt, {i)xu (j) - Xu ii)Xt, (j) ■ (255) 
The values presented in Tab. [T^ refer to the i = 1 case. By using Eq. (|249p we have: 

The last case we have to analyze is: I2 =Ia = 1- In this case, we can write: 

'^aSl(^' J') = «'i«'3 k'?'oi^o1o^(^, j) + (1 - ai,)T^l^f^{i,j)\ + (256) 
(1 - ai^)ai, [aij^{^fi{i,j) + (1 - auJT^l^f^i^j)] + 
ai,{l - aQ [ai,T^l^i{i,j) + (1 - ai,)T^l^^\i, j)] + 
(1 - a,J(l - ai,) [ai,Tllifi{i,j) + (1 - aQT^llf^ [i, j)] . 

The values of the matrix elements T of the previous equations are given in Tab. [15] for i = 1 
and j = 1,2. The other cases are calculated using Eq. ()249p . 
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1 


k 


u 


^ ± 


-7-41*2*3 
bo 


(1,1) 
6s 


64 


'7-41*2*3 

-'Zil/slis 


(1,2) 
6s 


64 











5 





-4 





-1 


2 








1 








-1 


6 





-4 


5 


-4 


4 


-4 





1 





-1 


-2 





4 


5 


4 


-4 


-4 








1 


-1 


6 





-4 


5 


-4 


-4 


4 


1 


1 





-7 


4 


4 





11 


-2 





-8 


1 





1 


17 


-12 


4 


-8 


-13 


14 











1 


1 


-7 


4 


4 





11 


-2 


-8 





1 


1 


1 


11 


-14 


8 


-4 


-7 


16 


-4 


-4 



Table 15: The values of the b coefScients of Eq. ()248p used to calculate the basic ^'j';*/ i (1, j) terms 
of Eq. ((2561) with j = 1,2. 



In order to simplify the notation in the calculation of the energy, we rename the matrix 
elements in the isospin space: 

^Sui^i^ = 'T^Xisi',^} (257) 

itfuJ^) = '^XtL.J^-S-.) (258) 

In the above definitions, the limitations related to the SOC approximation in the possibilities 
of linking isospin exchange operators, have been considered by defining L2 = (5ii^2 + (512^4 and 
L4 = 5iiU + Si2l2- This implies that the isospin term of Eq. (|246p with exponent I2 acts on the 
pair 13, and the term with exponent I4 on the pair 23. 
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D The uncorrelated one-body densities 



In this appendix we present the general expressions of the uncorrelated one-body densities 
involved in the calculation of the kinetic energy for systems not saturated in isospin. We use a 
set of single particle wave functions expressed as indicated by Eq. 
The one-body density p^l^ is given by: 



nljm 



where we have carried out the trace in the spin space and we have defined: 



(259) 



(260) 



/^T^ (ri,r2) 



For the one-body density matrices pT2 and PT2,j we find: 

: p*Hri,r2)VVo^(ri,r2) - Vip*Hri,r2) • Vip*^(ri,r2) 

: ^ Yl (2j + l)(2/ + l)<^,,.(r2)i?*V(r2)| 

^ ' nljn'l'j' I 

Rliijiri)D%^,in) - Rl^;^in)R%.,in)]Pi{cos9)P,{cos9) 



sm 



R%{n)Rl^n'f(^i)Pli^ose)P;,icose) 



r2 nlj 



(261) 



PT2j(ri.r2) 



1 ^ (-l)^-+^'-'-''-^i?i^,.(r2)<^,,,,(r2)( 

^ ^ nljn'l'j' (. 

Rni,in)D%.,{n) - R'^l^{n)R%,.,{n)\Qi{cose)Qi,{co^e) 
-^^^no(^i)<^'.,'(ri)QKcos0)Qj,(cos0)| , 



(262) 



where Pi are Legendre polynomials, 9 is the angle between the vectors ri and r2 and we have 
defined: 



Qi{cosi 
Q;(cos ( 



sin 0P/ (cos 9) , 



sin( 



cos9P[{cos9) - l{l + l)P/(cos6l) 



For the pT2, densities we have: 

Prsl^^i'^s) = 2V?/Jo^(ri,r2) 
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nlj 

PT3j(ri,r2) = 2V?4.(ri,r2) 

nlj 

The last density that appears in the calculation of the center of mass energy can 



as: 



PT6{ri,r2) - pI' (ri, r2)p3i5(ri, rs) 
Poj(ri,r2)/)ji5j(ri,r2) , 



where we have defined: 



pT(i{ri,r2) = 2(Vi/9*nri,r2)-V2P*o'(ri,r2) + 
Vip*'^.(ri,r2)-V24(ri,r2)) , 



2Vi-V2/9*^(.)(ri,r2) 



The explicit expressions of the above defined quantities are: 



nlj 



Rn!jin)Rl^!j(.r2)cosePi{cose) + 



Kljir2) 



Rnljir,) 



n 



+ Klj{n) 



sin^ e Pi {cos 9) + 



r2 



nlj^ 



rir2 



sin^ eP/(cos e) + l{l + 1) cos eP,(cos 6)^ 



nlj 



R%{ri)R'^l^{r2)coseQi{cose) + 



ri 



r2 



sin^ eQ'i{cose) + 



R'^l,iri)R%{r2) 



r\r2 



sin^ 6'Qj(cos6') + 



+ 1) cos eQi{cos 9) 

V sin 6 ' 

^ ' nlj 
n'l'j' 

cose R'^l.{r2)R%,^,{n)Pi{cose)Pv{cose) - 
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sm 8 i/ (cos 6')/;/ (cos y I 

rir2 ' ^ ^ " 



+ 



sin^e 



R%'l, An) Pi {cos 6) Pi, {cos 0) + 



ri 



Pl{cos9)Pl,{cos0) 



+ 



L_ y: {-iy^^'-'-'-'R'^i,{r2)R%,{n) 



2(27r) 



COS0 



nlj 
n'l'j' 



R%{r2)R'^n'j'{ri)Qi{cos e)Q,{cos 6) - 



sm 6 Qi{cos B)Q If {cos 0) 



rir2 



+ 



sin^^ 



r2 



R'^!i.Ari)Q'i{cose)Qr{cosO) + 



ri 



Qi{cos e)Q'i,{cose) 
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E The expressions of the energy expectation value 



The interaction energy defined in Sect. [Has Tp + V2 = W, can be expressed in terms of: 

H^jriru) = -^2(^|-|^<^,a[/p(n2)VVr(ri2 
6 

2~/p(''l2)/r-(ri2)('^r,5 + 5r,6)(l + Sp,5 + ^pfi) 

+/p(ri2)t''^(ri2)/r(ri2)} • (271) 

In the following, we shall use the separation of the spin and isospin operators as was done in 
Eqs. ()5ip and pOSp . We shall identify the various isospin independent operators by using the k 
and u indexes, which can assume the values 1,2, and 3. The I indexes can assume the values 
and 1. 

As in Sect. [H the expectation value of the interaction energy is calculated in four parts, see 
Eq. (|137p and Fig. [HI The first contribution, called Wq, is given by: 



1 



Wo = dr,dr2Hf^''^'^''''''^'''^'''''+''{r,2) 

Tklk2k3 Aks Jlt2 ^ \^tit2 

I Tk4kik5 Tk2k3k5 Aks Ak4 tlt2 ( i-„W*i*2 



(272) 



_|_ jk4kik5 jk2kskfi j-ksk(i2 ^k4 ptit2 ^2)Xl^ 



tlt2 
+l2+h+k 



(273) 



and corresponds to the case when p > 1 operators act between the external particles only. These 
operators are associated to fp, v'^ and fr functions in Eq. (|27ip for the direct terms. In addition, 
we should recall that the exchange terms have additional spin and isospin dependent operators 
because of the presence of the exchange operator n'^'^(l,2), see Eq. (j58|) . By definition, this 
operator may act on the bra or on the ket. In our conventions we make it act on the ket so 
it is always to the left of the rest of the operators. The sequence of the three operators we 
have mentioned above, plus the exchange operators, are present in all the terms we are going 
to analyze, therefore we shall always use the same set of indexes. Furthermore, we should point 
out that in Eq. ()273p a sum on all the k, I and t indexes is understood. This convention will be 
used in all the equations of this appendix. 

The expressions of the densities used in Eq. (|273p are given in Eqs. (|209ti21ip defined in 
AppendixlHl and include only state-independent vertex corrections. In this way, the p functions 
consider all the direct and exchange central dressing linked to the external points. As discussed 
in Sect. El ^'==1-2,3 = 1,3,6 and A'^ = 1 - 4,3. The values of I^^^^^'^ and j'^i'^zfca are given by 
Eqs. (jllOp and (|114p . The X;^*^ functions give the isospin traces and their values are given by 
Eq. dm]). 

We now discuss the effects produced by one SOR linked to one of the interacting particles. 
This is the Ws diagram of Fig. [TBI The operator structure that we must analyze, for the direct 
case, is: 

1 1 

-{/ci,ui}A:2-{A;3,it2} = 
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kiUlk2k3U2 /*=i^2'=3^fc3(l + 
klUik2U2k3 Jfeife2fc3^fc3(l + 



Table 16: Spin traces of the operators in Eq. (|274p . 



- {uikik2k3U2 + uikik2U2k3 + kiUik2k3U2 + kiUik2U2k3) , (274) 

In the above equation, the k and u indicate the operators and P", Eq. (15ip . The fci, ^2 and 
/ca operators act on the pair of particles 1 and 2. The ui and U2 operators act, instead, on the 
particles pair 1 and 3, or 2 and 3. The symbol {, } indicates the anticommutator. The values of 
the traces of the various terms of Eq. ()274p . are given in Tab. [T6l 



Ulklk2k3U2 /'«4felfe5jfe2/C3fe5^fe5(l + 

Ulklk2U2k3 /fc4fclfc5/fc2fc3fc5^fc5(l + Dk^uM'^'K 

klUlk2k3U2 /'=4fclfc5jfc2fc3fc5^fc5(l + Dk^^^^A^'^Su, 

klUlk2U2k3 Jfc4fclfc5jfc2fc3fc5^fc5(l ^ 



Table 17: Spin traces for the parallel spins case of the operators of Eq. (|274p . 



kiUikik2k3U2 Ii^ii^^>'5ii^2k3k,,jk5ke2^^ ^ Dk^u,)A'''5u^ 

kiUlklk2U2k3 /fc4fclfc5jfc2fc3fc6jfc5fc62j;L _^ Dfegm ) (5^^ 

kikiUik2k3U2 Ik4k,ksjk2k3k6jksk62^i ^ P'fcsnJA^Mni 

kikiUik2U2k3 lkAk,jk2k3k6jk5ke2^-^^ ^ Dfc^^JA"! 



Table 18: Spin traces for the case of antiparallel spins between particles 1 and 2 in Eq. (|274|) . 



In the exchange case, an additional A;4 operator must be included on the left hand side, 
following our conventions. Since we work with single particle wave functions expressed in a 
jj coupling scheme, antiparallel spin terms are usually contributing. For this reason in the 
exchange case, we have to distinguish the parallel and antiparallel spin situations. For the case 
of parallel spins between the interacting points, by following Ref. [21], we obtain results given 
in Tab. EZl 

For the antiparallel spin term, we have to distinguish the case when the antiparallel spins 
are those of interacting points, whose results are given in Tab. [TE[ from the case when they are 
those of the 1 and 3 particles, or those of the 2 and 3 particles. In this last case ui must be 
an exchange operator acting on the left hand side of the operator product. This produces the 
results given in Tab. [19l 

In the Tabs. [T71 [T8l \T9\ we have used the following values of the D terms: 

/ \ 

Dk„k2 = -4/3 -4/3 , (275) 
V -4/3 -4/3 / 
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Uikik2k3U2 

UlklU2k2k3 
UiU2kik2k3 



jkik2ki j^ks JUIU22 

Jfclfc2fc3^fc3/ni«22(l ^ 
jkik2k3 j^kz JU1U22 



Table 19: Spin traces for the case of antiparallel spins between the particles in 1 and 3 or 2 and 
3 in Eq. (pTIj) . 



which corresponds to the odd values of Epq defined in Tab. [13] of Appendix lAl 
All the isospin parts of the above operators can be written by using the xlljfi 
defined in Eq. (I257p . The contribution of the SOR to the interaction energy is: 



W, 



P2%M^r2) 



+ 
+ 
+ 

+ 



'9\fi^i^-2)C'^2{^2)+g'j^\r,,V2)&^^^^^^^^ 
)'^l'{^i^^2)C%{v,)+gif{v,,r2)Cl,,{^^^ 

)'afi^i^^2)C'^2{^2)+g'ti^i^^2)C^^^^^^^ 
'9'^l'{-i^-2)CU^i) + 9if{r,,V2)C%^^^^^^ 



^ jfc4fci k, jk2 fc3 fe5 j^k, A^' Va'eL (r 1 , r2 ) 

_l_ jkik-ik^ jk2kik^ jkz,k(i2 ^k4, ptit2 ^^^^ 



tlt2k2k4k5 ( \ 



where we have defined 

j^tit2kik2k3 
m,li,l2,l3 



ir^) 



j^tit2kik2kz 



ji^tit2k2k4,k^ 
d,ll,l2,l3,U 



ir^) 



^_o(r,) + (l+A3„jM^;o^;^V/2,/3(rO 

+ (l + A,„jM^yj^^,^^o(r,) 

+ (l + A,„jM^5^J^^,^(r,), 

= KTomIi.+i. (r.) + (1 + Dkzu, )<\ti+L^a (r^) 

+ (1 + I?.,„jM;^;™_,/3(rO + <\o7i+?.+/3,o(r«) ' 

= (1 + I).,„J(M£/-\^,^^,3^o(r,) + M*^/-V/2,Z3(rO) 

+ (l + Z).,„jM*£-+.^_,,^_o(rO 

+ (l + I),,.jM*^/-/^^,^_Jr,), 



(276) 



(277) 
(278) 

(279) 



Jijtlt2Ul / N 

^^^m,h,l2,hy'^^) 



A 



ui 



(1 '^«i,l)X/^o?20/3(^)^m,2«i-l,2ui-l(''«) 
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Direct 


Exchange 


uikik2ksU2 
uikik2U2k3 

kiUik2k3U2 
kiUlk2U2k3 
UiU2kik2k3 

kiUiU2k2k3 
kik2k3UiU2 

kik2UiU2k3 
kik2Uik3U2 


jkik2kQ jksk^kQ 
jk2k3k(i jkik^kQ 
jkskike jk2kske 
jkik2kg jkskske j^ke 
jkik2ke jkskskg j^ke 
jk2ksk(i jkikska 
jkik2ke jkskske j^ke 
jkik2kQ jk^k^k^ j^ka 
jkik2ke jkskske 


jkik2k(i jkQksk'Y jk4k5k7 
jkik2 fee jka k^k-j jk^k^ kg 
jk2 k:j ka jk^kikT jk-jk^ fcg 
jkskike jkekikr jk2k5k7 
jkikske jkik2k7 jkekrks j^ks 
jkikske jki k2kj jkgkrks ^fcs 
jk2k3kfi jkek4k7 jkikskr 
jki kske jkik2k7 jkekrks j^ka 
jk^k^k^ jkik2k-r jk^k-rks j^k^ 
jkikike jkek2k7 jkskskr 



Table 20: The traces obtained for the Wc{dd) term. 



^ Aiili20/3VV^»n,2wi,2u 
^ ^hll2ll3^^^rn,2ui,2u 



^_l(ri) + X/^o/2l«3(^)^m,2«i-l,2Mi' 



(280) 



j,Hlt2UlU2 ( \ 
^^^ej,h,l2,l3^^^) 



TUIU22 



^ MiOhlh^^^ ej,2ui-l,2u2y'-^^ ^ Mill2ll3y'')'-^ej,2ui,2u2\'-^) 



^«lO«20Z3 ,2ui-l,2u2-n^i) + X«il«20/3 ,2ui,2u2- 



lir^) 



41*2*3 



(281) 



In the above equations we have considered that i = 1,2, m = d,e , and fi = ti for i = 1 and 
H = t2 ioT i = 2. The expressions of the f^m,p<j(''«) ^l]*pgi^i) functions are given by Eqs. 
([2i5D . 

The structure of the Wc term is more involved. We calculate separately the various terms 
according to the direct or exchange nature of the correlations reaching the external points. 



Wc = Wc{dd) + Wc{ed) + Wc{de) + Wc{ee) + Wdcc) 
The operator structure that we must analyze in the dd direct case is: 



(282) 



^{ki,Ui}k2{k3,U2} + ^{kl,U2}k2{k3,Ui} + 



{{ui,U2},ki} + UikiU2 + U2kiUij ^2^3 + 
klk2({{ui,U2},k3} + Uik3U2 + U2k3Ui 



(283) 



where the various symbols have the same meaning as in Eq. (j274p . 

We calculate the value of this term by using the result of Eq. ()67p , and the results are given 
in Tab. [20j In this table, the factors ^'32'^*'^ present, since they will be included in the 

expressions of the nodal diagrams terms N which we shall define in analogy to Eq. (I7ip . The 
value of the spin traces are the same if we exchange ui and U2- For this reason, we have only 
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shown the results when ui is on the left hand side of U2- The isospin traces do not have this 
property and their values are given by the coefficients vllilili^i^i}) defined in Eq. (I258p . Since 
we have observed that the exchange contribution from antiparallel spins is much smaller than 
that from the parallel ones, we have neglected the contribution of the antiparallel spins in the 
SOC terms. 

By using the above definitions we can write: 



Wc{dd) 



+ A'=V2:eL(ri,r2: 



_l_ jkik2ke jk3k4kr jkrkske j^tit2k5 ^ ^ 

+ lk2k3kejk4k,krjkrk,ke^M2k5^^^^^^^^^^^^^ 

_|_ jk^kike jk^kikj jk2ksk-7 ^tit2k5 ^ ^ ^ 

_|_ jk^kske jkik2k7 jk(ik7k3 j^ks j\^tit2k5 ^ ^ ^ 

_j_ jk2k3k6 jkek4k7 jkikskT j^t^it2k5 ^ ^ ^^-^ 



where we have defined 



1 



dd,l\ ,l2,l3 ^ 

tlt2kz, 
I,dd,li,l2,l3 



^^£1, Jri, r2; 1) + M*;,*,^,':f,. Jri, r2; 2) 



.ri,r2j 



M2*it^./3(ri.r2) 

<.1,t^./3(ri'r2) 
<3i,^.,^3A(ri'r2) 



2 
3 

2 

^di;S,..+jri,r2) + M*:,;f^ 

3 

2 
3 

2 



dd.li ,l2,ls ^ 



dd,0,0,h+l2+l-J-^^'^'^^ 



ri,r2 + 



^dd,li,l2+l3fi^^^'^'^^ ^dd,0,h, 



(284) 



(285) 



(286) 
(287) 
(288) 

(289) 



and 



E E 



^/i,l,/2,l,/3 (^)-^?nn^2A;5,2Mi,2M2 ("^^ ' 



tl<2<3 
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+ ''^hJf2,0h^^')■^rnn,2k5-l,2ul,2u2-l^^^'^'^') 
+ ''^h,o/2^M^'''^'^mn,2k5-'i-,2ul-l,2u2^^^'^'^') 

+ (1 - 4,l)i\t40,^3(^)^^52l5-l,2«,-l,2«2-l(ri'r2) 



(290) 



The mn labels indicate that this last equation is valid not only for the dd case but also for the 
ed, de and ee ones. 

So far, we have calculated diagrams where only dynamical correlations reach the external 
points 1 and 2. The statistical correlations, labeled with e, are treated by using the spin-isospin 
exchange operator. In the ed part, the operator structure is: 



1 



(291) 



since ui is the spin-isospin exchange operator acting always on the left hand side. The various 
terms have already been calculated in the evaluation of the dd part (see Tab. [20j) . We obtain: 



Wcied) 



1 



dr,dr2Hfj^-'+'^'''''-'+''''''-'-'+''{n2)- 



(292) 



(5^^i^(ri,r2)C|(r2)+5^^/^(ri,r2)Cj22(r2))C'j22(ri 



jkik2 ke jks fcs ke j^e 



M*^55^^,^+,3,o(ri,r2; 1) + M:i;^«^,^^,^^,3(ri,r2; 1) 



41*2^5 



I Tkik2kQ jk^k^ke Tt^tit2k5 / r-c- ^\ 

-' ^'^ed,0,«i+Z2,i3'^ 1' 2> -LJ 

I rfc2fc3fc6 rfcifcsfce A/r*i*2fc5 j.^- ^\ 

^^^ed,0,h,l2+h^^'^' ^ 

The de term of Eq. (j282p has the same operator structure as the ed term when ui and U2 
are interchanged. This does not change the spin traces, therefore the result is: 



Wcide) 



dr,dr2H%'-'+^''''''''+''''''-'+''in2) 



9^1' (ri, r2)C^^2(ri) + 5*:,*^(ri, r2)C*^22(ri] 



(293) 



^^^!22(r2) 



jfci ^2 fce jkz fcs fee ^fee 



1 ^-tlt2fe5 
^'^(ie,0,/i+/2+«3,i 



o(ri,r2;2)+M*^5J^,^^,^^,^(ri,r2;2) 



I rfelfe2fc6 7-^3^5 ^6 A/ftlfefes /'„_, ,.„.9\ 

^'^de,0,«i+i2,«3'^ 1' 

_l_ Tk2k3kQ jkik^ka T\jtit2k5 / r^- 9'\ 

Since ui and U2 are spin-isospin exchange operators, in the ee term, the only possible ordering 
of operators is l/2{ui,U2}kik2k3. We obtain: 

W,{ee) = \ j dridr2//^^^-^+'^'''=^-^+'^'2'^-^+'^(n2) 

^^!22(ri)5dii*'(ri,r2)C*^22(r2)/'^'^'«/'^'^'M'=« 



M' 



tlt2k5 



{r,,r2-A)+MilXo:H_^l2+l,i^i,-2-,2) 



(294) 
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To calculate the Wc{cc) term of the eq. (j282p we found it useful to consider separately the 
situations where the p > 1 operators appear on the left, or on the right hand sides of the folding 
integrals. Specifically we define non-nodal diagrams as: 



2fiiri2)fp{ri2)+K^,%irur2) 



5*i(ri,r2) + 



(295) 



where the label Z can be L (for left) and R (for right). By using Eqs. (j219p and ()220p we define 
the left and right nodal diagrams as: 



q 



x!^c%a,3)C^l2'''C'/^,ri3)— X*3(3,2) + iV*3(3,2) -p*^(3,2) ) , 



A ^2 



(296) 



(297) 



The above equations (j295p . (|296p . (|297p form a set of hypernetted equations which can be 
solved iteratively. For example, one may start by setting the nodal diagrams equal to zero in 
Eq. (j295p . The (cc) nodal diagrams to be used in the evaluation of Wc{cc) are those where the 
left and right nodal diagrams are subtracted: 



iV-*/Hri,r2) = iV*^,p(ri,r2) -iV,^;^(ri,r2) -iVi;*;,(ri,r2) . 
The operator structure of the spin-tensor terms of the R diagrams is: 
1 



r{^4,iii} 



(298) 



(299) 



The terms related to the diagrams L are obtained by exchanging ui and U2. The various 
possibilities are given in Tab. [2TJ 





Traces 


kiUiU2kik2k^ 
kiUikiU2k2k3 
kiUikik2U2k3 
kiUikik2k-iU2 

Ml/c4M2fcl^2^3 

uikikiU2k2k^ 
uikikik2U2k^ 
uik4^kik2k^U2 


jk2k3k6 jkekikr jkikrks j^ks 
jk2k3k6 jkgkikr jkikrks 
jkik2ke,jk3k4kr jkrkak^ 
jk2k3 ke jkike k7 jk4kTks 
jfc2 fca kg jkikg kr jk^k^kc^ 
jk2k3 ke jk4,ki k7 jkgkrks 
jk\k2 kg jk^kekr jkrk3k5 
jk2k3ke jk4,kikT jk^k-jk^ j^k^ 



Table 21: Tensor-spin traces for the Eq. (|299l) . 
By putting together the various terms we obtain 
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+ 
+ 
+ 



jk2k3k(i jkgkikr j-kikTks j^t2k5 



jk2k3k6 jk(ik4kr jkikrks j^t2ks ^ ^ ^ 



cc,2k5 — l+l5 ^ 

> I 

CC,l4,,li,l2+l3 ^ 



Tkik2k(, jk-ikik-j jkrksks -iyjt2k5 / \ 

CC,l4,ll+l2,l[i^ 1' 2J 



I jk2k-ik(, rfcifcefcy jkik-jk^ M^'^^^ (y. yr,\ 

I Tk2ksk(i jkikik-j jkekTks Mhk^ (r r \ 

^ ^ ^^-'cc,0,i4+ii,«2+;3*^ 1' 2J 

I Tkik2k(i jkik^k-j jkjk^k^ A/f'a'^s 

^ ^ ^^-'cc,0,i4+ii+i2,i3*^'^^l' "^^j 

I Tk2k3k(i Tk4kik7 j-kekTks Aks -i\/ft2k5 / \ 



with 



M 



tiks 



CC,li,l2,l3 



ri,r2 



njRtiks . 



ri,r2;l)+M^^4';^^,3(ri,r2;2) 



(300) 



(301) 



In the above equation the R and L functions are defined as in Eq. (j290p by substituting the 
nodal diagrams N with the left and right nodal diagrams of Eqs. (I296p and (I297p . 



We give now the expression of the kinetic energy terms. For the t'^^ term we obtain: 



rr = -|^/*i/'T\(ri)C^Sn(ri) 



(302) 



where p^^(ri) has been defined in Eq. (1132^ and a sum on ti is understood. 
As indicated in Sect. 14. H we separate the remaining terms in three parts: 



^<t> — ^ <Pfi ~^ ^ <f>,s ~^ (t>,c n — 
In order to express the above quantity in a closed form, we define the function: 



JT^^ 3dd (ri, rsj - ()p,i()r,i 



(303) 



(304) 



where we have included the scalar dressing of the correlation operator acting on the external 
particles. We can see an analogy with the kinetic energy part in Eq. (j27ip . In addition, we have 
to consider a spin-isospin exchange operator acting on the external particles. By associating 
the indexes ki and li to p, /c3 and to r and and 1^ to the exchange operator, we can use 
the results obtained for the traces of the interaction energy in the case (? = 1, corresponding to 
k2 = I and I2 = 0. 



Consequently, for Tj"J we get: 



T 



(2) 



4>,o 



4m 



(305) 



Jfc4fclfc3^fc3ptlt2(r^^r2) + /^4fclfc5jfe5fc32^tlt2^ (j.^^ j.^) 



+ 
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26-, 



tlt2 



p^f (ri,r2)+/9^*^,(ri,r2] 



r(3) 



'2^ / *1*2( /l2fci-l+h,2fc3-l+i3V"^'"^^"" -^ii+ia+U 



(306) 



/'^'^'M'^VTf (ri, r2)iV*^ (ri, rs) + 

p^f (ri, r2)iVi^')*^ (r2, n) + ^^^,(^1, r2)ivg*^ (r2, ri)] (7*^22(^2) + 
p^^*^ {r^,r2)N(^>^ {r,,r,) + ^^^,(^1, r2)ivig*^ (r2, ri)^ 

^£22(^2)-!]}) , 

where the one-body densities pT2, PT2,j, Pts and prsj have been defined in Appendix [Dl A 
comparison with Eq. (j273p shows that we have only the exchange terms. We have substituted 
P2,exc(j) with PTn,{j) ^nd we have used I^^^^^ = Sk^k^- The term without operators that appears 
in the above equations is caused by the different vertex correction that must be used when no 
dynamical correlations reach the particle 2. 

Using the property mentioned above for the / matrices, we get the following expressions for 
the SOR contributions: 



.(2) 



4m 



j dridr2|^4^_i+,,2fe3-i+/3(ri'r2)A'=^ 
/^•4^'i^M'=Vt2' (ri, 1-2) + /''"^^'''/''''''VT2j-(ri, r2' 



(307) 



+ 



26-, 



p^2'(ri, 1-2) + p^r^l^ATl,Y2) C%2{'^l) 



Utsoci^l) Ct22i^2)-l + C'i,22(^2)Ulsoci^2 



■.(3) 



2m 



dvidr2 



^^2iti+h,2.3-i+/3(ri'r2)A'=^ 



(308) 



jfc4fc,fc3^fc3^^it2(r,,r2)iV*^(ri,r2) + 

jk,k,k, Ik,k,2pt,t2^ ^ ^^^^t2^ (^^ . 

Kt}!:it>i) + Kl.^%>^)] + 25M.Cj,,(ri) 

p^*^ (n , r2)iVi^ (ri , r2) + p^^^- (ri , r2)iV,% (ri , r2) 

C^j22(r2)f/^50c(r2) + f/feoc(ri; 

{v,,V2)N^'- (r2, ri) + p^f ^.(n, r2)ivg*^ (r2, r^: 



C. 



t2 

d,22 



r2 + 
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p^f (n, r2)iVi^')*^ (r2, ri) + /,^f^.(ri, r2)ivig*^ (r^, n) 
C%2(r2) - 1 



(3) 

In the evaluation of the T^^c terms we neglect T^^. This approximation is justified by the 

(3) (2) 

fact that the contribution of the T^q g terms are much smaller than those of the T^q g terms. 

(2) 

The structure of the ^ term is analogous to that of the exchange case of Wc{dd) when k2 = 1 
and I2 = 0. In this case, we use J^^^k? — A^'^Sj.^k^, and obtain the expression: 



.(2) 



48m 



j dri(ir2P^2''(''l'^2)/2fci-l+/i(n2)/2fe3-l+/3(^12) 



A'=^C*:,,(rOC*^22(r2)5^^i^(ri,r2) - 



I jkzkik^ jkekski 71^*1*2^:5 „\ , jkikikf, jkak^k-j, 71^*1*2^:5 „\ 



I Tk4k3k(i rkekiks Aks -n,rtit2k5 / \ 

I jk^k^ka jk^kik-i Aks T\jtit2ks / ^ 

+ ^ ^ ^ ^^^4.,dd,li,0,l3,U\^i^^2' 

+ ^'^'^'^•/'^'^'^^^iSzi,.3,o(ri,r2)|. 



kjkikfi jkiksks /\/f*i*2fc5 



(309) 



We present now the contribution of the spin-orbit terms. In our calculations we use correla- 
tions up to the tensor channels, therefore the presence of spin-orbit operators is only due to the 
interaction. We consider only the case of spin-orbit operators acting on the external particles. 
In this case we write: 



< >o 



9 J dridr2 



h+h {ri2W^^^ {ri2)h+iz (n2) 



/f(n2) 



P2%M^-M\ 



*1*2 



+^2+^3 



+ 



(310) 



Jl*2 



il*2 



*1*2 



where we have used the relation: 



C 



pf2i(L -8)12^1 =-mk,z5kzz . 



(311) 



Finally, we analyze the expressions used to calculate the contribution to the energy of the 
three-body potential, Eq. (ll4ip . The diagrams we consider in this calculation have been presented 
in Fig. [T9j We start by considering the term of the three-body force, Eq. (jl43p . which is a 
scalar function. The diagram (3.1) of Fig. [T9l is the leading term, and its contribution can be 
expressed as: 

< ^^f23 >3.i= \ j dridr2dr37;f23/'3*'*'(ri,r2,r3) , (312) 



109 



where a sum on the t indexes is understood, ps is the three-body density that can be written as: 

r2, rg) = pl]%t' i^i^^2,rs) + 6t,t,KtA,eJc (ri, r2, rg) 

mm',n7i',ZZ'=(i(i,(ie,ed 

+25t,t,5j,t35*^,(ri,r3)Kc'(r3)5*^,(r3,r2)V;*|(r2)5*^(r2,ri) (313) 

where l^'„(ri) has been defined in Eq. (jl98p . We have separated the direct and the exchange 
parts of the three-body density, and we have neglected the contribution of the Abe diagrams 
[95j . These Abe diagrams are simple non-nodal ones with three external points that play the 
same role as elementary diagrams in the TBDF. 
The expressions of the diagram (3.2) is: 

<^f23>3.2 = - /"dridr2dr3^^'=-^+'^^"^^^^^'=^-^+'^^^^^^ 



^'^5;.„;.,x[j+\/^35*'(ri,r2,r3)+ (314) 

A'=^/'=^'=^'=M'=^X^t\+Z3/^3:e^c^(ri'r2,r3)^ • 

where, as in all the previous equations, a sum on the operator channels, here indicated by the 
indexes k and I, is understood. As indicated in diagram (3.2) of Fig. UM the various operators 
act between the 1 and 2 particles. In this respect, we should point out that the scalar operator 
term, i. e. the case ki = k2 = 1 and h = h = 0, is not considered in Eq. (j314p . since this would 
produce a double counting with the contribution of the (3.1) diagram, Eq. (j313p . 

We consider now the contribution of the vf^s term of the three-body force, Eq. (jl42p . We 
define first two effective potentials: 



eff,mn 



^ r r 

(ri'^^) = E E4^2. / dr3[5^*nri,r3)C*|(r3)^/*f (r3,r2) 

fci,fc2=2 ts 

+ 5^^dnri,r3)C*^22(r3)5*f (r3,r2) 
+ 5^^*3(r,,r3)C*3,,(r3)5*f (r3,r2)^ 

C^3t'^X'=Hn3)^'^(r-32) (315) 
3 

,cc(ri>r2) = Yl E^^2. / *35cc(ri,r3)Cj22(r3)5?c(r3,r2) 

fci,fc2=2 is 

Cf3t'^X'=Hn3)^'^(r32) (316) 

with m,n = d, e and X^{r) defined in Eq. (jl45p . With the help of the above defined quantities 
we express the contribution of ^^23' related to the (2.1) diagram of Fig. [191 



k3,t 
eff,cc 



< VU3 >2.1 = ^ E / / 

tlt2 



f2ki-i+h (ru) f2k2-i+h (^12) 
fKri2) 



Tkik3k2 /|A;2,,iii2 

J ^ ^/l+/2 + l 



110 



+ g*^;^(ri,r2)C*^(ri)Cj22(r2) +5*;,*^(ri,r2)Cj22(ri)C|(r2) 



+ V, 



+ 



eff,de 
Jit2 I 



ri,r2)Q22(ri) 
+ ^e%*J(ri,r2)Cj,,(ri)[5*^,*^(ri,r2)C*^^ 

+ 9^l*^(ri,r2)Cj22(r2)^ 

+ ^e1fci^(ri,r2)9^;,*^(ri,r2)Cj22(ri)Cj22(r2) 

- ^ /'^''i / (i?r2/2fci-i+«i(ri2)/2fc2-i+/2(r-i2)5di'(ri,r2) 
• C7j22(ri)Cj22(r2)/'=^'=^^V'=^^^'=^^'=HxJ^t\+i + ^^1^+2)^'=^ 



ri,r2; 



iV*^(ri,r2)-p*Hri,r2: 



[iV*Kri,r2)-/9*^(ri,r2)] 
+ 2f5>%(ri, r2) k^ri, r2) - p*^(ri, r2) 



(317) 



In the above expression, we have considered only the contribution of the anticommutator term. 

In the (2.2) diagram of Fig. [T9t the operators act on the 13 and 23 particle pairs. In this 
case, the operator structure is: 



(318) 



where we have defined: 



(0?3>0^1)± = 0?30^1±02lO?3 



The contributions of Eq. (j318p . with p = Ik^ — l + li,q = 2ki, si = 2/c4 — I + /2 and S2 = 2A;2, 
are given in Tab. \T2\ where we used the matrix n^^i^^k^ki [ig£]2ed as 



R 



kik2k-j. ki 



T)k\ ■pk2 T)kz Tj,kA\ 
^ \^Vi -^23 -'^13 -'^23 ) 

A^Hkrk-A^^^k2kA(^ + Dk2k^ + 12P2(n3 • ^23) X 
1 - (1 - 5fc,3)(l - 433)] [1 - (1 - 423)(1 - <5fe43) 



(319) 



The plus signs in the table correspond to the anticommutator terms and the minus signs to the 
commutator ones. 



Ill 





Traces 


kikz{ki,k2)± 
{ki,k2)±ky,ki 
(ki, k2)±kikz 
ksiki, k2)±k4, 
ki{ki, k2)±k3 


A'^^dk^k, A^^5k,k,± i2'=4fc3fc2fci 
Rk,k2k3k,^ A^3<5fc3fci A^^5k,k, 
A^-^6k3k, A^^5k2k,± R^^^^^-^^^ 
A^Hk^k, A^^5k,k,± i2'=3fc2fcifc4 
^fc4fcifc2fc3± A^Hk^k, A^^5k2k, 



Table 22: Tensor-spin traces of the operator of Eq. pisp . 



When the isospin traces are included we obtain: 



Qili2i3± 

kik2k3k4,lil2 



R''^k,k,k2^A''Hk,k,A'''6k,k2 
Xili2iio(.'^) ^ Xili^oiii'^) =•= Xoi2?i+iio(^) 

±2x^l?o*il+i(2) 



(320) 



where we have defined the ^^k^k2k%4, symbols which allows us to express the contribution of 
diagram (2.2) of Fig. [19] as: 

A2n 



1-n 

< Vi23 >2.2 



dridr2dr; 



/2fc3-l+^l(n3)^fcl^^^^^^j 



/2fc4-l+Z2(''23) 



/i(n3) 



/l(r23) 



'^32j 



Pi 



1 



h ~^ ^^kik2k3k4,,lil2 



(321) 



where the sum over all the k, I and t indexes is understood. 

In the diagram (2.3) of Fig. \T9\ the operator dependent correlations act on the 1,2 and 1,3 
particle pairs. This implies that the only contribution comes from the commutator term of Eq. 
(jl42p . In this case, the operator structure is: 

i (uo'is, oii}[ol„o'2l] + hPls, o^l]{o?3> 



The spin and tensor traces of this set of operators are given in Tab. 
written a function Ci23^'^^ that is present in all the traces. 

After considering the isospin traces we obtain the expression: 



tit3t2 _ tit3t2 



l^jk^kiks _ jkzk^k^j^kr, 

+ioa,'ii(i) + ');fe!.u(2) - iSiS'ii(2) - 4'f.+iio(2) 



(322) 



where we have not 



»,*l*3t2 I 



+^0011(2) 



2^*1*3*2 



i/20io(2) + 2<tno(2) + 2^^01^,(1) 
11 



9„ili3*2 
'^%201/l- 



1^ 



(323) 



112 





Traces 


[ki, k2]k3ki 
[ki, k2]k4k3 
k3[ki,k2]k4 
k4[ki,k2]k3 


jfe3fcife5_ jkgkik^j^ks 
_jk'ikiks_^ jk^kik^j^ks 
jk2,kik^_ jkzk-^kz, J^kr, 
_jk3kik5 jk^kik^ j^ks 
_jk-ikik5_^ jkskiks j^k3 
jk3kik5_ jk^kik^ j^ks 



Table 23: Tensor-spin traces of the operators of Eq. p22p . 



Calling ^^^^ the above trace, we can write: 

<.123>2.3 - —Jdr,dr2dr3 ^^^^^^^ X (^3) ^^^^^^^ 

X'=^(r32)Cil3'^'=^pr^*^(ri, r2, rs) (324) 
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F The Euler procedure 

In this Appendix we present the method used to fix the correlation function. We have named 
this method Euler procedure. The starting point is the calculation of the hamiltonian mean 
value, considering only cluster terms up to the second order. In this case, the contribution to 
the W terms of Eq. (jl37p is given only by the Wq term, since the other terms are produced 
by clusters of higher order. In analogy, only the terms ' contribute to the r$ term of Eq. 
(jl28p . Therefore, the total energy in the two-body cluster approximation is given by: 

E2 = W2 + ri,^J + ri,'J , (325) 

where the expression of the various terms are: 

W2 = iy"dridr2if^^i-^+'i'''=^-^+'^'''^-^+'^(ri2) 

|p^^^^*Hri,ri)p^3^^*^(r2,r2)x^,(l)x+3(2)^f2^^i'l^i'|Xs2(l)x.4(^ 
X^(l)x*t(2)(Ti-T2)^^+'^+'^X*x(l)x*2(2) 

X+(l)Xit(2)(Ti • r2)'^+'^+'«x*2(l)x*.(2)| , (326) 
= "i^/ ^ri4r*^(ri,ri) J dr2pg^^^*^(r2, rs) 

(/2fci-l+h(n2)/2fc2-l+/2(n2) - '^2fci-l+/i,l'^2fc2-l+Z2,l) 

X+(l)xi(2)^5^f|Xs2(l)Xs4(2) 

X,t(l)Xit(2)(^i • r2)'^+'^Xt.(l)Xt2(2) , (327) 
= |^/dridr24--3-*^*^(ri,r2) 

(^/2A;i-l+«i (n2)/2fc2-l+/2(n2) - '^2fei-l+«i,l'^2fc2-l+/2,l) 

X+(l)Xit(2)(^i • r2)'^+'^Xt2(l)Xt,(2) . (328) 

As in the calculation of the energy expectation values, the k indexes can assume the values 1, 
2 and 3, and the I indexes the values and 1. The s = ±1/2 indexes indicate the spin third 
component and the t = ±1/2 the isospin third component. A sum on all the third components 
of spin and isospin and on the repeated indexes is understood in the above equations, and in 
the following ones. 

The other quantities of the above equations are the two-body correlation functions fp, Hj'p 
defined in Eq. (j27ip and the isospin independent operators, see Eq. (09]). We use the 
expressions of the one-body densities ()132p and ()133p where the spin dependence is explicit: 

p^\(ri) = ^p^\^^*^(ri,ri)x:,(l)Xs2(l), (329) 
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(ri,r2)V?pr'*^(r2,ri 



-Vip^i^^*i(ri,r2)-Vipg3^^*^(r2,ri) . (330) 
We obtain the optimal correlation functions fp by solving the Euler-Lagrange equation: 

<5(^2-C2) 



5fp 



, 



(331) 



where we have indicated with C2 the contributions of the constraints. The expression of C2 is 
analogous to that of W2 after substituting H^fp with fpXqfr, \ being the Lagrange multipliers. 
The values of these multipliers, are fixed by imposing the conditions that the various terms of 
the two-body correlation function assume their asymptotic values after a certain internucleon 
distance dp, called healing distance: 

fi{r >di) = 1 , (332) 



fp>i{r > dp) 



In addition we impose the condition: 



dfp 
dr 







(333) 



(334) 



Since we use correlation functions composed by six operator channels, the minimization 
procedure ()33ip should be applied to fix the values of six healing distances dp. This produces a 
system of six interconnected differential equations. It is possible to separate these equations by 
using a representation of the Euler equations in terms of the total spin and isospin S and T of 
the correlated nucleon pair. We use the projection operators: 

1 



nf2 = ^[25 + 1 + 



1 



[2r + i + (-i)^+Vi -Ts 



0-2] 



(335) 
(336) 



with 5, T = 0, 1. In the (T,S) representation we can write: 



Y^XpO{^= {xTS + Ss,lXTtSl2)Ilj2^f2 , (337) 

p=l 5,T=0,1 

where Xp can be the scalar part of the correlation, fp, that of the interaction v^, or the Langrange 
multiplier Xp. The relation between the expressions of these quantities in the two representations 
is: 



xsT = a;i + (4r-3)3;2 + (45-3)x3 + (4r-3)(45-3)x4 
xtT = X5 + (4r - 3)X6 . 



(338) 
(339) 



As discussed in Sect. 15.31 in our calculations we have used only two healing distances, one 
for the four central channels, dp=i,2,3,4 = dc, and the other one for the two tensor channels, 
(ip=5 g = dt. In terms of these quantities, we can rewrite the boundary conditions as: 



fsrir > dc) 
fTt{r > dt) 



1 



dfsT 



dr 



dr 



, 



. 



(340) 
(341) 



r=dt 



115 



After some algebra, and by using the properties of the Pauh matrices, we obtain for E2 the 
following expression: 



E2 



SlS2S3S4tlt2 



dridr2 

{Jt s{r 12)^'^ hsir 12) - (V/Ts(r-i2))^) 



4m 



+8(55,1 ( ( -VTi{ri2) - VTt{ri2) ) /rt(''i2) + m(n2)/ri(n2)/Tt(ri2) 



Am 

SlS2tl / 



/Tt(ri2) (vVTt(ri2) - ^/Tt(ri2)) - (V/Tt(ri2))' ) 



p^i^^^i (n, ri)p^3.4r2(r2, r2) - (-1)^ ^^^^^^^^ {^^,^2)9'^''' (ri, r2) 

/ / dv2{fh{ri2)+^Ss,iiUri2) - l) 
p^Y^*Hri,ri)p^^^^*^(r2,r2) - {-lf+^pf^^^'''''''\v^,V2) 



(342) 



In order to make the variation on £"2, we found convenient to rewrite the above expression 
as a function of two new quantities -Pt5 (ri,r2) and <5rs(ri, r2), defined as: 



^'T5(ri,r2) = 



2T+1 



<5Ti[/^(ri)pg(r2) + pS(ri)p^(r2) 

+(-l)^4(pg(ri,r2K(ri,r2)+pS(ri,r2K(ri,r2)) 
+i[pg(riK(r2) + p^(ri)pg(r2) 

-(-lf+^4(pg(ri,r2)p^(ri,r2)+pS(ri,r2)pg(ri,r2)) 

, i6(-ir r 
(2S + i)(2r + i) . 

(^Ti {fij (ri , r2)p[}j- (r2, ri) + pgj- (ri, r2)p^j- (r2, ri)) 
+ \ (Poj (r 1 , r2 )poj (r2 , ri ) + p[J^- (ri , r2 (r2 , ri )) 



(343) 



tit2=p,n 



(p^i(ri)p*^ (r2) - 4(-l)^+^ (p*^ (ri, r2)V?p*^ (r2, ri) 
-Vip*i(ri,r2)-Vip*^(r2,ri)))" 
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+ 



16(-1) 



{2S 



(344) 



-Viy9*Uri,r2)-Viy9g.(r2,ri) 



+ -PTs{T^l,r2)vTs{ri2) , 



where we have explicitly written the sum over the spin and the isospin. In order to obtain a 
quantity depending only on the relative distance ri2 between the particles 1 and 2 we integrate 
Eq. (j342p over ri and r2 by keeping fixed the value of ri2. We define the quantity: 



Pts = / dr2Prs'(ri,r2) 



27r 



r2dr2 



\r2+ri2\ 



\r2-ri2\ 



ru Jo 

and an analogous expression for QtSj therefore we write E2 as: 



ridriPTsiri,r2) 



(345) 



dvi2 QtsIts - -^PTs[fTs^ Its - hsY 



+85s,i I ( Qti - Privxt + -^^Pti ) frt + PxivrfTifTt 



(346) 



mr 



12 



--^Pti ( fTtV'fTt - {VfTtT ]]+Qts + ^Ptsvts 



where we must understand that all the functions and the operators act on the ri2 coordinate. 
The expression of the variation of E2 with respect to fxs ^ind fxt is: 

dvi2 5fTsPTs[ - -^^^TS + {VtS - \ts) FtS 

{vTt - \Tt) Frt] + ^5s,i5fTtPl'H - -V^Frt + 
) I m 



Vti — \ti — "^VTt + 2At'( + 



mr 



12 



FTt + 



{vti — Xti) -Frij 



, 



where we have defined 



Ets 
Frt 

Vts 



ItsPts > 



pTS 



2Qts + 



4m 



Pts 



(347) 

(348) 
(349) 

(350) 



The fact that Eq. (13470 has to be valid for every variation of fxs or fxt, implies that both 
expressions included in the braces has to be zero. By imposing this condition we obtain the 
Euler-Lagrange equations: 



m 



V'^Fts + ( Vts - ^ts ) Fts + Ss,i {vTt - ^Tt) Frt = , 
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V^Ftj + Vti - An - 2vTt + 2ATt + — ^ Frt 

m \ mrf2 ) 

+ (vTi - Ati)Fti = . (351) 

The expressions ()35ip represent a system of differential equations, the two equations corre- 
sponding to S = are not coupled, while the other equations are coupled. The solution of the 
above equations gives the optimal value for f^s and /rt and by using Eq. (j337p we obtain the 
correlation functions /p, in the representation used in the FHNC calculations. 

In the case of infinite systems, nuclear and neutron matter, these equations have been gen- 
eralized to include the spin-orbit channels in the correlation with v% [21] and f 14 [97] potentials. 



118 



G Acronyms 



AFDMC Auxiliary Field Diffusion Monte Carlo 

AV8' Argonne Vg nucleon-nucleon potential 

AV18 Argonne vis nucleon-nucleon potential 

CBF Correlated Basis Function 

FHNC Fermi HyperNetted Chain 

FHNC/SOC Fermi HyperNetted Chain/Single Operator Chain 

GFMC Green's Function Monte Carlo 

HF Hartree-Fock 

HNC HyperNetted Chain 

IPM Independent Particle Model 

MBCF Many-Body Correlation Function 

MF Mean Field 

NN Nucleon-Nucleon 

OBDF One-Body Distribution Function 

OBDM One-Body Density Matrix 

RFHNC Renormalized Fermi HyperNetted Chain 

QCD Quantum ChromoDynamics 

SOC Single Operator Chain 

SOR Single Operator Ring 

SRC Short Range Correlations 

TBCF Two-Body Correlation Function 

TBDF Two-Body Distribution Function 

TBDM Two-Body Density Matrix 

U14 Urbana vu nuclcon-miclcon potential 

UVH Urbana VH thrcc-nuclcon interaction 

UIX Urbana IX three-nucleon interaction 
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H Symbols 



Symbol Meaning Note 

^ ^ folding product Eq. ((MI) 

< I > Clebsch-Gordan coefficient 

A number of particles 

Og , parameters of the WS potential Eq. ([861) 

BP C-trace of two operators Tab. [10] 

C(ri2) sum of the composite diagrams Eq. (|20l) 

Clof^) C-trace Sect. [2XT] 

C^(l) vertex corrections Eq. (I203II204I] 

Cm.pi}(l) vertex corrections Eq. (j238ll239() 

cc sub- index for cyclic-cyclic Sect. 12.21 

Z)J,,^.(r) Eq. mil 

dd sub- index for dynamical-dynamical Sect. 12.21 

de sub- index for dynamical-exchange Sect. 12.21 

E{ri2) contribution of elementary diagrams 

Epq matrix used in SOC calculations Tab. [13] 

ee sub- index for exchange-exchage Sect. 12.21 

F{xi, ....,xa) many-body correlation function Eq. ([2]) 

^(1...^) operator dependent correlation function Eq. (|57l) 

/(r) two-body correlation function Eq. ^ 

fp{r) state dependent correlation function Eq. ([48]) 

g{xi,X2) two-body distribution function Eq. (jlOj) 

5*1*2(1^2) State dependent TBDF Eq. (l53l) 

Hij state dependent part of the TBCF Eq. ([561) 

Hf;{n2) Eq. dsn]) 
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h(r) 


two-body h-function 


Eq. (HD 


hc(r) 


used in SOC Eqs. 


Eq. dZZD 


hri(r) 


used in SOC Eqs. 


Eq. (SB 


jijk 


spin part of K^'''' 


Eq. (Ilinl) 


jijk 


spin part of L^'^'' 


Eq. (|114|) 




used in the product of two operators 


Tab. E] 


L{ri2) 


used in SOC Eqs. 


Eq. dZH]) 


jjpqr 


used in the product of four operators 


Tab. M 


l{x) 


Slater function 


Eq. (IMl) 






Eq. (12771I28TD 




pion mass 


Eq. §6!) 


N 


number of neutrons 






contribution of nodal diagrams 








Eq. (12T91I22D 


nljTTl V / 


used in the quasi-hole wave function 


Eq.(1172P 




momentum distribution 


Eq.(|l55D 


interaction and correlation operators 


Eq. (09]) 


Piicosd) 


Legendre polynomials 




p{ 


proton / neutron projector operator 


Eq. 


pk 

''J 


Q2k-1 


Eq. dsn) 


U ' IS ' ^ 


parameters of the WS potential 


Eqs. iSlETD 


nlj V / 


radial part of 


Eq. (USD 


J^kik2kzk4, 


spin C-traces 


Eq. (l:-}2()j) 




spatial coordinate of the particle i 




5(ri2) 




Eq. (|2()5|) 


Sij 


tensor operator 


Eq. jSOD 


gp,n 


one-body sum rule 


Eq. (|147|) 


gpp,pn,nn 


two-body sum rule 


Eqs. (|148D 


S2,a 


spin sum rule 


Eq. (|149|) 
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spectroscopic factor 


Eq. (jl8()j) 


s 


symmetrizer operator 






rp{i) rp rp 

' ''-c.m. 


kinetic energy terms 


Eqs. 




'7-fir2^3 (A\ 
t'lt2t3t4t5 / 


used in the isospin traces 


Eq. (|25lD 


used in the isospin traces 


Eq. (j246|) 




mean-field potential 


Eqs. 


(|85I86P 


X \ / 




Eqs. 








Eqs. 


(|24UII245]l 




volume of the system 






u ■ is 


parameter of the WS potential 




(lob]) 


z;^(ri2) 


scalar parts of the interaction 






,,27r 


part of the three-body force 


Eq. 


A-i A Ah 

(|i42p 


^123 


part of the three-body force 


Eq. 


/1 1 /I 1\ 

P4d|) 




contribution of non-nodal diagrams 


Fn 

i!jq. 








_cjq. 






used in the quasi-hole states 








generalized coordinate of particle i 








spherical harmonics 








spin spherical harmonics 


Eq. 




Wo Ws Wc Wcs 


interaction energy terms 


Eq. 




z 


number of protons 








p2fc-l 






Ap{l,...,p) 


sub-determinant 


Eq. 


(ini) 




used in the isospin traces 


Eq. 


(|258j) 


pp 


matrix used in SOC Eqs. 


Tab 


. m 


V 


spin-isospin degeneracy 


Eq. 




tW 
S123 


function used to calculate SOR 


Eq. 


mi 


^kiksksjhh 




Eq. (j32,S|) 


^k\k2k3 
''123 


spin part of ^^23 


Eq. 
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spin-isospin exchange operator 


Eq. 


(158]) 


p 


particle density 






PT1,..,T4 


kinetic energy densities 


Eqs. 


, (fT32lfT35D 


PT5,T6 


kinetic energy densities 


Eqs. 


, (12661 1267P 




OBDF 


Eq. 


Pi 


parallel spin OBDM 


Eqs. 


■ (|91|9.S|) 


Poj 


antiparallel spin OBDM 


Eqs. 


, (I92l9;-i|) 


pf^*^(l,2) 


operator dependent TBDF 


Eq. 


(iioip 






Eq. 


(I2U9I- 12111) 


p''''*''*(ri,r2) 


OBDM 


Eq. 


(|15ip 


«>(1,...,A) 


independent particle wave function 








single particle wave function 








natural orbit 


Eq. 


(|158j) 


Xs,t 


Pauli spinors 


Eq. 




An 


used in the isospin traces 


Eq. 


(|121D 




correlated many-body wave function 








quasi-hole wave function 


Eq. 


(|168|) 




the polar angles 6i and 






r\tit2t3± 

kik2k3k4,lil2 




Eq. 


(I320|) 
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